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1,  INTRODUCTION  AND  SUMMARY 


Year-round  daylight  saving  time  (YRDST)  has  recently  been 
observed  in  the  United  States.  The  observance  of  double  daylight 
saving  time  (DDST)  is  under  consideration.  One  of  the  principal 
expected  effects  of  the  adoption  of  these  time  systems  is  a change 
in  the  level  of  motor  vehicle  accidents.  For  this  reason,  Project 
LATEDAY  was  initiated  to  determine  the  Level  of  Accidents  under 
The  Effect  of  Daylight  saving  All  Year. 

The  two  principal  objectives  of  Project  LATEDAY  were: 

1.  The  determination  of  an  estimate  of  the  level  of  motor 
vehicle  accidents  under  the  effect  of  year-round  daylight 
saving  time. 

2.  The  determination  of  an  estimate  of  the  level  of  motor 
vehicle  accidents  under  the  effect  of  year-round  double 
daylight  saving  time. 

These  objectives  were  accomplished  through  a computer  analysis 
of  currently  available  motor  vehicle  accident  data.  An  experiment 
was  designed  which  made  use  of  these  data  to  model  the  accident 
situation  in  the  United  States  under  year-round  DST  and  under  DDST. 
A novel  technique,  the  Daylight  Saving  Equivalent  Day  approach,  was 
developed  to  allow  data  from  non-YRDST  and  non-DDST  years  to  be 
used  to  find  an  estimate  of  the  effect  of  YRDST  and  of  DDST. 

The  results  of  Project  LATEDAY  show  that  under  YRDST  there 
would  be  an  increase  in  morning  accidents  and  a somewhat  smaller 
increase  in  evening  accidents,  as  compared  to  a six-month  DST/six- 
month  standard  time  system.  Based  upon  these  results,  an  increase 
in  yearly  motor  vehicle  accidnets  of  the  order  of  one-half  to  one 
percent  is  estimated  if  YRDST  is  introduced. 

Under  DDST,  the  results  of  Project  LATEDAY  show  a large  in- 
crease in  morning  accidents  and  a moderate  increase  in  evening  ac- 
cidents compared  to  the  six-six  system.  Based  upon  these  results, 
an  increase  in  yearly  motor  vehicle  accidents  of  the  order  of  four 
percent  is  estimated  if  DDST  is  introduced. 
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It  should  be  noted  that  these  results  are  only  estimates,  based 
upon  data  that  were  not  necessarily  fully  representative.  Further- 
more, the  Equivalent  Day  methodology,  though  possibly  the  best 
available  technique,  was  based  upon  a set  of  assumptions  which  may 
not  be  entirely  satisfied.  Thus,  the  results  should  be  considered 
only  as  an  indication  that  there  would  be  a greater  number  of  motor 
vehicle  accidents  under  YRDST  or  under  DDST  than  there  would  be 
under  a policy  of  six  months  DST  and  six  months  standard  time,  as 
was  previously  in  effect. 
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2,  DAYLIGHT  SAVING  TIME  IN  THE  UNITED  STATES 


The  idea  of  advancing  clock  time  to  make  the  best  use  of  day- 
light seems  to  have  originated,  as  many  good  ideas  have,  with 
Benjamin  Franklin  in  1784.  Not  until  World  War  I,  however,  did 
pressures  for  fuel  conservation  and  increased  national  efficiency 
finally  produce  action.  The  Standard  Time  Act  of  1918  officially 
divided  the  country  into  time  zones  and  established  seven  months 
of  daylight  saving  time  (one  hour  advanced)  from  the  last  Sunday 
in  March  to  the  last  Sunday  in  October.  The  Interstate  Commerce 
Commission  was  empowered  to  fix  time  zone  boundaries  and  to  govern 
time  regulation. 

Following  the  war,  the  daylight  saving  time  provision  was  re- 
pealed, and  from  1920  to  1942  DST  was  at  local  option.  During 
World  War  II,  Year  Round  DST  (called  "War  Time")  was  introduced, 
primarily  for  power  conservation  reasons.  "War  Time"  was  repealed 
at  the  close  of  the  war  in  1945.  From  1946  to  1966,  DST  was  again 
at  local  option,  with  each  locality  free  to  choose  whether  to  have 
DST  and  when  the  DST  period  should  begin  and  end. 

In  1966,  the  Uniform  Time  Act  was  adopted  to  standardize 
the  use  of  DST  throughout  the  country.  The  Act  provided  that  all 
states  were  to  be  on  daylight  saving  time  from  the  last  Sunday  in 
April  until  the  last  Sunday  in  October.  Exemptions  were  allowed 
only  for  complete  states,  by  acts  of  their  legislatures.  Under  an 
amendment  to  the  Act,  state  legislatures  of  states  that  are  in  two 
time  zones  may  also  exempt  the  entire  area  of  the  state  lying  within 
one  of  the  time-zones. 

With  its  establishment  in  1967,  the  Department  of  Transporta- 
tion took  over  from  the  Interstate  Commerce  Commission  the  function 
of  fixing  time-zone  boundaries.  It  also  assumed  responsibility  for 
policing  the  observance  of  the  Uniform  Time  Act. 

In  December  1973,  under  the  impetus  of  the  Energy  Crisis,  the 
Emergency  Daylight  Saving  Time  Energy  Conservation  Act  of  1973  was 
passed.  It  provides  that  full-year  daylight  saving  time  shall  be 
in  effect  in  the  entire  country,  from  January  1974  until  April  (and 
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therefore  October)  1975.  Exemptions  are  allowed  only  for: 

1.  One  time  zone  of  a two-or-more  time  zone  state, 

2.  A state  in  one  time  zone  and  not  contiguous  to  any  other 
state,  and 

3.  Hardship  declared  by  proclamation  of  the  Governor  of  a 
state  prior  to  the  effective  date  of  the  Act. 


The  Act  directs  the  Secretary  of  Transportation 
ports  to  the  Congress  on  the  effects  of  this  Act 
mendations  for  legislative  or  other  action:  an 

due  June  1974,  and  a Final  Report,  due  June  1975 


to  submit  two  re- 
, including  recom- 
Interim  Report, 
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3,  EFFECTS  OF  LIGHT  ON  MOTOR  VEHICLE  ACCIDENTS 


It  is  generally  accepted  that  the  presence  of  daylight  de- 
creases the  motor  vehicle  fatality  rate  and  the  absence  of  daylight 
increases  the  motor  vehicle  fatality  rate.  A National  Safety 
Council  report  shown  in  Table  3-1  shows  a higher  fatality  rate  at 
night  than  during  the  day.  Other  factors  (e.g.,  differences  in 
driver  speeds  or  driver  intoxication)  may  be  responsible  for  some 
of  this  disparity,  but  the  light  effect  is  probably  at  least  a par- 
tial cause. 


TABLE  3-1.  DEATHS  AND  DEATH  RATES  BY  DAY  AND  NIGHT 


URBAN/RURAL 

DAY/ NIGHT 

PERCENT 
OF  DEATHS 

DEATH 

RATES* 

Day 

4 7% 

3 . 1 

Total 

Night 

53% 

8.3 

Urban 

Day 

4 7% 

1.9 

Night 

53% 

5.2 

Rural 

Day 

47% 

4.3 

Night 

53% 

11.6 

^Deaths  per  100,000,000  vehicle  miles 

Source:  "1972  Accident  Facts"  National  Safety  Council 


YRDST  causes  the  light  pattern  to  shift  one  hour  back  as  com- 
pared to  standard  time.  Thus  morning  periods  that  would  have  had 
light  under  standard  time  are  dark  under  YRDST.  and  evening  periods 
that  would  have  been  dark  under  standard  time  are  light  under  YRDST. 
With  this  change  of  light  pattern  should  come  a corresponding  change 
in  the  accident  rate,  so  that  the  affected  morning  hours  should  have 
a higher  accident  rate  and  the  affected  evening  hours  should  have  a 
lower  accident  rate.  Under  the  assumption  that  the  traffic  pattern 
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(vehicle -miles  per  hour)  remains  unaffected  by  YRDST,  the  new  ac- 
cident rates  can  be  applied  to  the  daily  traffic  pattern  to  find 
an  estimate  of  the  total  accidents  under  YRDST. 

However,  it  is  not  clear  that  the  daily  traffic  pattern  re- 
mains unchanged  under  YRDST.  There  are  three  classes  of  effects 
that  might  be  expected.  First,  the  shift  in  daylight  of  one  hour 
might  cause  a complete  shift  in  activity  of  one  hour  for  people  in 
"dayl ight - oriented"  activities.  This  might  include  persons  in  oc- 
cupations such  as  farming,  construction  work,  etc.,  where  light  is 
more  important  than  clock  time,  and  persons  not  working,  such  as 
retired  people  and  those  on  a day  off  from  work.  Second,  there 
may  be  a time  change  of  certain  trips.  With  an  extra  hour  of 
light  in  the  evening,  some  trips  such  as  midday  shopping  may  be 
deferred  to  the  evening.  Third,  new  trips  may  be  generated.  For 
example,  an  extra  hour  of  light  in  the  evening  may  induce  some 
additional  recreational  travel,  such  as  an  evening  swim,  baseball 
game,  or  tennis  match,  that  might  not  have  been  scheduled  if  dark- 
ness occurred  an  hour  earlier. 

For  this  reason,  an  approach  to  the  determination  of  the  level 
of  accidents  under  YRDST  (and  similarly  under  DDST)  should  not  dis- 
count a change  in  traffic  pattern.  As  will  be  seen,  the  approach 
used  in  this  study  finds  days  upon  which  the  sunrise/sunset  (and 
therefore  the  traffic  patterns  affected  by  them)  are  the  same  as 
they  would  be  under  YRDST  or  under  DDST,  and  uses  these  days  for 
comparison  purposes. 
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4,  AVAILABLE  DATA 


The  best  source  of  motor  vehicle  accident  data  available  to 
this  project  was  the  set  of  Standard  Accident  File  Extract  (SAFE) 
Computer  Tapes  which  were  prepared  under  contract  to  the  National 
Highway  Traffic  Safety  Administration  (NHTSA) . These  tapes  contain, 
in  a standardized  format,  information  concerning  each  reported 
motor  vehicle  accident  in  a state  for  a given  year.  The  original 
accident  data  were  taken  from  state  Accident  Data  Tapes  (ADTs) , 
which  were  supplied  to  NHTSA  by  the  individual  states.  The  con- 
tractor prepared  individual  computer  programs  to  convert  the  data 
from  each  ADT  into  the  SAFE  standard  format. 

The  SAFE  standard  format  specifies  the  form  in  which  the  data 
for  each  accident  are  placed  on  the  SAFE  tapes.  Each  group  of  re- 
cords of  the  tapes  contains  data  for  one  accident.  The  data  are 
broken  into  several  categories: 

1.  Modified  Accident  Severity  Index  (e.g.,  total  number  of 
people  severely  injured) 

2.  Accident  Related  Factors  (e.g.,  date,  time) 

3.  Vehicle  Identification  (e.g.,  make,  model) 

4.  Vehicle  Dynamics  (e.g.,  speed) 

5.  Roadway  Data  (e.g.,  number  of  lanes) 

6.  Occupant  Data  (e.g.,  driver,  age,  sobriety) 

7.  Non-occupant  Factors  (e.g.,  pedestrian  age,  sobriety) 

8.  Special  Data  (e.g.,  destination  for  buses). 

An  analysis  of  the  data  available  and  the  methodology  to  be 
used  by  this  project  determined  that  data  from  only  the  first  two 
categories  would  be  needed.  Since  tape  processing  on  a computer 
is  costly,  it  was  decided  to  have  excerpted  SAFE  tapes  produced 
specifically  for  Project  LATEDAY.  These  LATEDAY  tapes  would  con- 
tain, for  each  accident,  only  those  records  of  the  SAFE  tapes  in 
which  the  necessary  data  resided. 
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LATEDAY  tapes  were  produced,  under  contract  to  the  sponsor, 
for  all  state-year  pairs  for  which  the  corresponding  SAFE  tape 
data  were  available  and  usable.  Most  were  produced  from  SAFE 
tapes,  but  a few  were  produced  directly  from  state  ADTs. 

After  eliminating  all  LATEDAY  tapes  that  had  obvious  errors 
or  other  problems,  there  remained  a total  of  20  LATEDAY  tapes  which 
could  be  used  in  this  project,  as  listed  in  Table  4-1.  Note  that 
the  accident  data  are  all  for  either  1971  or  1972. 

The  LATEDAY  tape  format  for  each  accident  is  shown  in  Table 
4-2.  It  consists  of  one  36-byte  record.  The  first  half-record 
is  a copy  of  the  first  18  bytes  of  the  Modified  Accident  Severity 
Index.  The  second  half-record  is  a copy  of  the  first  18  bytes  of 
Accident  Related  Factors.  Only  those  bytes  used  by  the  LATEDAY 
programs  are  indicated  in  Table  4-2.  The  other  bytes  are  in- 
cluded to  adhere  to  the  SAFE  tape  format,  and  the  corresponding 
SAFE  tape  information  usually  does  appear  there.  However,  for 
those  LATEDAY  tapes  produced  directly  from  the  ADTs  these  bytes  may 
be  blank. 

It  should  be  noted  that  the  SAFE  tape  format  (and  thus  the 
LATEDAY  tape  format)  shows  the  time  of  each  accident  as  local  time 
without  any  indication  of  the  time  zone.  The  location  of  the 
accident  within  the  state  is  also  not  included  in  the  SAFE  informa- 
tion. Therefore,  for  a state  which  is  split  into  two  time  zones, 
it  is  impossible  to  determine  in  which  of  the  time  zones  an  ac- 
cident occurred.  Since  knowledge  of  the  time  of  the  accident  re- 
lative to  sunrise  or  sunset  is  necessary  to  determine  a DST  effect, 
the  SAFE  data  for  states  split  into  two  time  zones  may  not  be 
acceptable  for  use  in  this  project.  Of  the  eleven  such  states 
for  1971-1972,  only  a few  had  SAFE  data  available.  Of  these,  only 
two  were  used:  Kansas,  which  had  over  991  of  its  population  on 

one  side  of  the  time  zone  boundary,  and  Texas,  which  had  over  96% 
of  its  population  on  one  side  of  the  time  zone  boundary. 
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TABLE  4-1.  LATEDAY  TAPES  USED 


1971  (5) 

District  of  Columbia  (DC) 

Maine  (ME) 

South  Carolina  (SC) 

Utah  (UT) 

Vermont  (VT) 

1972  (15) 

New  Hampshire  (NH) 
New  Mexico  (NM) 
Oklahoma  (OK) 
Pennsylvania  (PA) 
South  Carolina  (SC) 
Texas  (TX) 
Washington  (WA) 


California  (CA) 
Connecticut  (CT) 

Delaware  (DE) 

District  of  Columbia  (DC) 
Georgia  (GA) 

Kansas  (KS) 

Michigan  (MI) 

Missouri  (MO) 


TABLE  4-2.  LATEDAY  TAPE  FORMAT 


Bytes 

Contents 

Code 

1-2 

Unused 

3 

Total  number  of  fatalities 

'O'  to  '9'  = 0 to  9 
'A'  to  'X'  = 10  to  28 
'Y'  = over  28 

4-7 

Unused 

8 

Total  number  of  non-fatal 

Same  as  Byte  3 

injuries . 

9-19 

Unused 

20-21 

Year 

For  example  '72'  = 1972 

22-23 

Month 

'01'  to  '12'  = January  to 

December 

24-25 

Day  of  month 

'01'  to  '31'  = 1st  to  31st 

26-28 

Unused 

29-32 

Military  clock  time 

'0000'  to  '2359'  = Midnight 

to 

11:59  p.m. 

33-36 

Unused 

For  all  Bytes: 

* = unknown/not  stated 
blank  = not  reportable 
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5,  DAYLIGHT  SAVING  EQUIVALENT  DAY  METHODOLOGY 


The  principal  technical  problem  faced  in  Project  LATEDAY  was 
that  of  using  data  which  came  from  years  when  neither  Year  Round 
Daylight  Saving  Time  nor  Double  Daylight  Saving  Time  were  in 
effect  to  determine  the  impact  of  YRDST  and  DDST  upon  motor  vehicle 
accidents.  A novel  approach,  the  Daylight  Saving  Equivalent  Day, 
was  developed. 

One  method  for  investigating  the  effect  of  YRDST  or  DDST  on 
motor  vehicle  accidents  would  be  to  compare  the  number  of  accidents 
in  two  similar  periods,  one  on  Standard  Time  and  the  other  on  Day- 
light Saving  Time  or  Double  Daylight  Saving  Time.  There  were  two 
obstacles  to  such  a comparison:  first,  as  mentioned,  data  were  not 

available  for  two  such  periods;  second,  even  if  such  data  were 
available,  changes  in  mot  or -vehicle  accidents  observed  between  the 
two  periods  may  be  caused  by  factors  other  than  Daylight  Saving 
Time,  e.g.,  weather  patterns. 

To  overcome  these  obstacles,  the  Daylight  Saving  Equivalent 
Day  Approach  was  devised.  It  enables  a comparison  to  be  made,  in 
effect,  between  a Standard  Time  period  and  a YRDST  or  DDST  period. 
In  addition,  it  allows  the  impact  of  YRDST  or  DDST  upon  motor 
vehicle  accidents  to  be  isolated  from  the  changes  due  to  other 
factors . 

5.1  THE  DAYLIGHT  SAVING  EQUIVALENT  DAY 

The  approach  used  in  this  project  is  based  upon  the  Daylight 
Saving  Equivalent  Day.  This  and  the  following  two  sections  will 
consider  this  concept  as  applied  to  the  study  of  Year  Round  Day- 
light Saving  Time.  The  application  to  Double  Daylight  Saving  Time 
will  be  discussed  in  Section  5.4. 

The  accident  data  available  to  this  project  were  for  years 
(1971  and  1972)  which  had  six  months  of  DST.  The  sunrise  and  sun- 
set curves  for  such  years,  for  a typical  location  in  the  United 
States,  are  shown  as  solid  lines  in  Figure  5-1.  If  these  years  had 
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SUXSET 


Figure  5-1.  1972  Sunrise  and  Sunset,  40°N  Latitude 
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Year  Round  DST,  the  sunrise  and  sunset  curves  would  follow  the 
broken  lines  in  Figure  5-1.  The  problem  is  therefore  to  use  acci- 
dent data  which  are  available  only  for  days  that  have  their  sunrise 
and  sunset  on  the  solid  line  to  estimate  the  number  of  accidents 
on  days  that  have  their  sunrise  and  sunset  on  the  broken  lines. 

The  Daylight  Saving  Equivalent  Day  Approach  was  designed  to  do 
this,  as  will  now  be  discussed. 

To  determine  the  impact  of  imposing  YRDST  upon  a Standard  Time 
day,  a Daylight  Saving  Equivalent  Day  corresponding  to  the  given 
day  is  found.  The  Daylight  Saving  Equivalent  Day  is  defined  as  a 
day  which  has  its  sunrise  one  hour  later  than  the  sunrise  of  the 
given  day  and  which  has  its  sunset  one  hour  later  than  the  sunset 
of  the  given  day.  Thus,  it  is  the  equivalent  of  the  Standard  Time 
day  having  DST. 

Using  the  definition  of  Equivalent  Day,  the  basic  idea  behind 
the  Daylight  Saving  Equivalent  Day  Approach  can  be  stated  as 
follows:  If,  for  a given  Standard  Time  day,  such  an  Equivalent  Day 

can  be  found  from  among  the  actual  days  for  which  data  are  available, 
then  the  number  of  accidents  that  were  recorded  on  the  Standard  Time 
day  can  be  compared  against  the  number  of  accidents  recorded  on  the 
Equivalent  Day.  All  other  things  being  equal,  the  only  difference 
between  the  two  days  would  be  the  one  hour  difference  in  sunrise 
and  sunset  time.  Thus,  the  comparison  would  show  the  effect 
of  having  DST  on  the  Standard  Time  day.  Finally,  if  such  a com- 
parison were  made  for  all  the  Standard  Time  days  of  the  year,  the 
overall  comparison  would  show  the  impact  of  having  Year  Round  Day- 
light Saving  Time. 

The  above  is  an  idealized  description  of  the  Equivalent  Day 
Approach.  Practically,  there  are  two  major  points  that  must  be 
considered.  First:  do  Equivalent  Days  exist  for  the  Standard  Time 

days  of  the  data  years  and  how  should  they  be  chosen?  This  will  be 
discussed  in  the  next  section.  Secondly:  what  should  be  done  if 

"all  other  things"  are  not  equal?  This  is  discussed  in 
Section  5.3. 
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5.2  SELECTION  OF  EQUIVALENT  DAYS 


In  this  section,  the  desirability  of  using  composite  Daylight 

Saving  Equivalent  Days  and  several  considerations  in  their  selection 
are  discussed. 

5.2.1  Composite  Equivalent  Days 

Examination  of  the  sunrise/sunset  curves  of  Figure  5-1  reveals 
that  there  are  very  few  Standard  Time  days  in  the  year  for  which  a 
day  exists  that  has  both  a one  hour  later  sunrise  and  a one  hour 
later  sunset.  This  is  due  to  the  approximately  inverse  shapes  of 
the  sunrise  and  sunset  curves  (i.e.,  a day  with  a later  sunrise 
than  a given  day  usually  has  an  earlier  sunset)  . Since  no  single 
Daylight  Saving  Equivalent  Day  exists  for  most  of  the  Standard  Time 
days  of  the  year,  it  was  decided  that  composite  Equivalent  Days 
would  be  used. 

A composite  Daylight  Saving  Equivalent  Day  is  constructed  out 
of  two  different  days,  one  with  its  sunrise  one  hour  later  than  the 
Standard  Time  day,  and  one  with  its  sunset  one  hour  later.  For 
example,  consider  March  14  (a  Standard  Time  day  for  the  data 
years  1971  and  1972)  in  a city  such  as  Denver  (latitude  40°)  . 

Sunset  occurs  at  6:06  p.m.  (see  Figure  5-2).  The  Daylight  Saving 
Equivalent  Day  for  this  day  would  therefore  have  a sunset  time  one 
hour  later,  7:06  p.m.  A search  for  such  a day  yields  September  17, 
which  is  selected  as  the  evening  portion  of  the  Daylight  Saving 
Equivalent  Day  (or  the  Evening  Equivalent  Day) , as  shown  in 
Figure  5-2.  Similarly,  sunrise  in  Denver  on  March  14  occurs  at 
6:14  a.m.  (see  Figure  5-3).  The  Daylight  Saving  Equivalent  Day 
thus  would  have  a sunrise  time  one  hour  later,  7:14  a.m.  A search 
for  such  a day  yields  October  19,  which  is  selected  as  the  morning 
portion  of  the  Daylight  Saving  Equivalent  Day  (or  the  Morning 
Equivalent  Day)  for  March  14  for  Denver,  as  shown  in  Figure  5-3. 

The  composite  Daylight  Saving  Equivalent  Day  is  thus  composed  as 
follows: 
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Figure  5-2.  Determination  of  the  Evening  Portion 
of  a Daylight  Saving  Equivalent  Day 
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SUNRISE 


MONTH 

r 

Figure  5-3.  Determination  of  the  Morning  Portion 
of  a Daylight  Saving  Equivalent  Day 
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Standard  Time 
Day 

Daylight  Saving 
Equivalent  Day 

Sunrise 

Sunset 

6:14  a. m. 

March  14 

6:06  p .m . 

7:14  a.m. , Oct.  19 
7:06  p ,m . , Sept . 17 

Thus  data  on  accidents  for  March  14  are  compared  with  data  for 
the  morning  portion  of  October  19  plus  data  for  the  evening  portion 
of  September  17. 

5.2.2  Existence  of  Morning  and  Evening  Equivalent  Days 

Using  the  sunset  curve  of  Figure  5-4  it  is  possible  to  find 
the  Evening  Equivalent  Day  corresponding  to  each  Standard  Time  day  of 
the  data  years.  These  are  shown  in  Figure  5-4.  It  can  be  seen 
that  an  Evening  Equivalent  Day  exists  for  each  Standard  Time  day; 
i.e.,  for  approximately  180  days  per  year.  For  several  Standard 
Time  days,  there  are  two  possible  Evening  Equivalent  Days.  The 
method  for  choosing  among  two  or  more  possible  Equivalent  Days  is 
discussed  in  Section  5.2.7. 

Similarly,  using  the  sunrise  curve  of  Figure  5-3,  the  Morning 
Equivalent  Day(s)  corresponding  to  each  Standard  Time  day  of  the 
data  years  are  shown  in  Figure  5-5.  A Morning  Equivalent  Day 
exists  for  about  one-third  of  the  Standard  Time  days  of  the  year; 
i.e.,  approximately  60  days  per  year.  For  most  Standard  Time  days, 
there  are  two  or  more  possible  Morning  Equivalent  Days. 

5.2.3  Latitude -Longitude  Considerations 

The  determination  of  Morning  and  Evening  Equivalent  Days  is 
based  upon  the  shape  of  the  sunrise  and  sunset  curves  for  the  region 
of  interest.  However,  a change  in  longitude  at  a given  latitude 
does  not  change  the  shape  of  the  curves  significantly;  it  just  in- 
creases or  decreases  all  of  the  values  by  approximately  the  same 
constant,  as  shown  in  Figure  5-6.  Therefore,  as  the  example  in 
the  figure  indicates,  two  regions  with  the  same  latitude  will  have 
the  same  set  of  Morning  and  Evening  Equivalent  Days. 

Latitude  changes,  however,  do  not  change  the  shape  of  the 
sunrise/sunset  curves,  as  illustrated  in  Figure  5-7.  This  figure 
shows  that  at  different  latitudes  the  same  day  can  have  significantly 
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POSSIBLE  EVENING  EQUIVALENT  DAYS 


DEC 


NOV. 

OCT. 

SEP. 

AUG. 

JUL. 

JUN. 

MAY 

APR. 

MAR. 

FEB. 

JAN. 


NOV.  DEC.  JAN.  FEB.  MAR.  APR.  MAY 
STANDARD  TIME  DAYS 


Figure  5-4. 


Evening  Equivalent  Days  Corresponding 
to  the  Standard  Time  Days  of  the  Year 
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POSSIBLE  MORNING  EQUIVALENT  DAYS 


JUN 


NOV.  DEC.  JAN.  FEB.  MAR.  APR  MAY 
STANDARD  TIME  DAYS 


Figure  5-5.  Morning  Equivalent  Days  Corresponding 
to  the  Standard  Time  Days  of  the  Year 
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SUNSET  LONG.  75° W SUNSET  LONG.  85°W 


SUNRISE 
LAT  40 °N 
LONG.  8 5° W 


Figure  5-6.  Equivalent  Days  for  Two  Different  Longitudes 
at  the  Same  Latitude 
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SUNSET  LAT.  30 


SUNSET 
LAT.  45°N 
LONG.  90°N 


Figure  5-7.  Equivalent  Day  for  Two  Different  Latitudes 
at  the  Same  Longitude 
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different  Equivalent  Days.  Therefore,  as  shown,  the  Daylight  Saving 
Equivalent  Day  is  a function  of  latitude.  Thus,  for  each  region  of 
interest,  the  central  latitude  is  used  to  determine  the  set  of 
Equivalent  Days,  and  longitude  is  not  considered. 

5.2.4  Base  Region  for  Accident  Data 

Since  the  data  are  available  by  state  and  year,  it  was  decided 
that  the  basic  region  for  use  in  accident  comparisons  would  be  the 
state.  The  geographical  center  of  each  state  was  found,  and  the 
states  were  tabulated  by  the  latitude  (to  the  nearest  whole  degree) 
of  their  geographical  centers  (see  Table  5-1).  This  latitude  was 
used  to  specify  the  sunr ise/sunset  curves  which  were  used  to  find 
the  Equivalent  Days  for  each  state. 

The  central  latitude  of  a state  will,  for  an  average  state,  be 
no  more  than  2°  of  latitude  from  the  latitude  extremes,  with  a max- 
imum of  about  5°  for  all  states  except  Alaska.  Data  on  the  loca- 
tion of  an  accident  within  a state  are  not  available  using  SAFE 
data.  If  such  data  became  available,  some  slight  gain  in  accuracy 
might  be  made  by  breaking  up  states  with  large  latitude  spreads 
into  two  or  more  regions. 

5.2.5  Same-Day-of-the-Week  Equivalent  Days 

In  comparing  the  number  of  accidents  on  a Standard  Time  day 
against  the  number  of  accidents  on  its  Equivalent  Day  there  may  be 
a problem  if  the  two  days  occur  on  different  days  of  the  week. 

This  is  illustrated  in  Table  5-2,  which  shows  the  average  number 
of  accidents  for  each  day  of  the  week  for  three  states.  It  can  be 
seen  that  there  can  be  large  differences  between  the  numbers  of 
accidents  on  two  different  days  of  the  week.  Therefore,  a dif- 
ference in  accidents  between  a day  and  its  Equivalent  Day  may  be 
due  to  the  difference  in  day  of  the  week  and  not  due  to  DST. 

To  eliminate  this  problem,  it  was  decided  that  an  Equivalent 
Day  would  only  be  chosen  to  be  the  same  day  of  the  week  as  its  cor- 
responding Standard  Time  day.  Thus,  if  the  Standard  Time  day  is  a 
Thursday,  then  both  its  Morning  Equivalent  Day  and  its  Evening 
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TABLE  5-1.  GEOGRAPHIC  CENTER  LATITUDES 


LATITUDE 

STATE 

47° 

Montana,  North  Dakota,  Washington 

46° 

Minnesota 

45° 

Maine 

44° 

Idaho,  Michigan,  South  Dakota,  Vermont,  Wisconsin 

43° 

New  Hampshire,  New  York,  Wyoming 

42° 

Iowa,  Massachusetts 

41° 

Connecticut,  Nebraska,  Pennsylvania,  Rhode  Island 

O 

o 

Illinois,  Indiana,  New  Jersey,  Ohio 

39° 

Colorado,  Delaware,  District  of  Columbia,  Maryland, 

Nevada,  Utah 

04 

oo 

o 

Kansas,  Missouri,  West  Virginia 

37° 

California,  Kentucky,  Virginia 

36° 

35° 

Arkansas,  North  Carolina,  Oklahoma,  Tennessee 

34° 

Arizona,  New  Mexico,  South  Carolina 

33° 

32° 

Alabama,  Georgia,  Mississippi 

31° 

Louisiana,  Texas 

04 

O 

o 

29° 

Florida 
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TABLE  5-2.  AVERAGE  NUMBER  OF  1972  INJURY-CAUSING  ACCIDENTS 
BY  DAY  OF  THE  WEEK  FOR  THREE  SELECTED  STATES 


Day  of  the  Week 

California 

Kansas 

Vermont 

Monday 

448.7 

52.4 

10.9 

Tuesday 

439.5 

54.1 

10.2 

Wednesday 

446.5 

54.3 

10.4 

Thursday 

464.0 

55.2 

13 . 1 

Friday 

634.9 

76 . 2 

19.0 

Saturday 

582 . 2 

72.2 

19.1 

Sunday 

442.5 

42.3 

11.8 

Equivalent  Day  are  chosen  to  be  Thursdays.  This  introduces  a slight 
inaccuracy,  since  the  Equivalent  Day  may  not  have  a sunrise  (sunset) 
exactly  one  hour  later  than  that  of  the  Standard  Time  day.  However, 
whatever  day  of  the  week  the  Equivalent  Day  would  normally  have 
been,  changing  it  to  another  specified  day  of  the  week  requires  a 
shift  of  no  more  than  three  days.  Since  sunrise  (sunset)  times 
typically  change  between  zero  and  two  minutes  per  day,  the  new 
Equivalent  Day  will  usually  be  no  more  than  three  to  five  minutes 
off  from  the  desired  sixty  minute  difference  between  Standard  Time 
day  and  Equivalent  Day.  To  eliminate  any  large  deviations  from  the 
desired  sixty  minute  difference,  we  can  consider  that  no  Equivalent 
Day  exists  whenever  the  best  same -day-of -the -week  Equivalent  Day 
has  a deviation  more  than  a specified  amount  (e.g.,  five  minutes). 

5.2.6  Holidays 

It  is  well  known  that  the  traffic  pattern  on  a holiday  often 
differs  greatly  from  the  "normal"  traffic  pattern  that  would  have 
occurred  on  the  same  day  if  it  were  not  a holiday.  Thus,  if  a 
Standard  Time  day  and  its  Equivalent  Day  are  being  compared  and  if 
either  or  both  are  holidays,  a difference  in  accidents  may  well  be 
due  in  large  part  to  the  unusual  holiday  traffic  patterns  rather  than 
being  due  to  a DST  effect.  Since  there  are  typically  just  eight  to 
ten  legal  holidays  per  year  in  a state  it  was  decided  that  it 
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would  be  reasonable  to  eliminate  the  holidays  from  consideration 
rather  than  to  try  to  adjust  for  holiday  traffic  behavior.  Since 
the  traffic  on  the  evening  before  a holiday  is  often  affected  by 
the  holiday,  these  were  also  eliminated. 

Each  state  was  contacted  to  determine  the  legal  holidays  for 
the  data  years.  All  of  these  holidays  and  all  the  evenings  before 
the  holidays  were  then  eliminated  from  both  the  list  of  Standard 
Time  days  and  the  list  of  possible  Equivalent  Days. 

5.2.7  Choosing  Between  Multiple  Equivalent  Days 

As  mentioned  in  Section  5.2.2,  and  shown  in  Figures  5-4  and  5-5, 
more  than  one  Morning  or  Evening  Equivalent  Day  often  corresponds  to 
a given  Standard  Time  day.  An  example  is  shown  in  Figure  5-8. 

Since  it  was  convenient  in  the  processing  to  have  a unique  Morning 
Equivalent  Day  and  a unique  Evening  Equivalent  Day  corresponding  to 
each  Standard  Time  day,  a method  for  choosing  between  multiple 
Equivalent  Days  was  developed. 

One  method  of  choosing  between  multiple  Equivalent  Days  would 
be  to  choose  the  Equivalent  Day  nearest  to  the  Standard  Time  day. 
However,  the  major  non-DST  differences  between  two  days  that  can 
affect  accidents  are  differences  related  to  weather.  Therefore, 
it  is  preferable  to  choose  the  Equivalent  Day  closest  to  the  Standard 
Time  day  in  weather. 

As  a gross  approximation  to  the  yearly  weather  pattern,  the 
following  function  was  used: 


As  shown  in  Figure  5-9,  this  function  will  be  +1  for  February  1 
(Day  number  32  of  the  year)  which  might  be  considered  the  middle  of 
winter.  Also,  WF  will  be  -1  on  August  2,  approximately  the  middle 
of  summer.  It  will  be  near  zero  for  spring  and  fall.  Thus,  when 
there  are  multiple  Equivalent  Days,  the  Equivalent  Day  which  has  its 
Weather  Function  closest  to  that  of  the  Standard  Time  Day  will  be 
chosen . 


WEATHER  FUNCTION 
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SUNRISE 
LAT.  40°N 


MONTH 

Figure  5-8.  Multiple  Equivalent  Days 
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JFMAMJJASONDJ 


MONTH 


Figure  5-9.  The  Weather  Function 
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As  an  example  of  this  process,  consider  the  multiple  Morning 
Equivalent  Days  shown  in  Figure  5-8.  This  figure  shows  that  for 
April  1,  at  40°N  latitude,  there  are  three  possible  Equivalent 
Mornings  (ignoring  same-day-of - the-week  Equivalent  Days  for  this 
example):  February  21,  September  20,  and  November  15- 

The  WF  value  for  April  1,  the  Standard  Time  day,  is  found  to 
be  0.5148.  The  WF  values  for  the  three  possible  Morning  Equivalent 
Days  are  found,  yielding  the  following: 


DAY 

WF (DAY) 

| WF (DAY) -WF ( APRI L 1) | 

February  21 

. 9416 

0.4268 

September  20 

-.6666 

1.1814 

November  15 

. 2297 

0.2851 

Thus  November  15  has  the  closest  WF  to  that  of  April  1 and  is 
chosen  as  the  Morning  Equivalent  Day.  Note  that  though  February  21 
is  only  40  days  from  April  1,  whereas  November  15  is  98  days  away, 
November  15  is  chosen  since  it  should  be  somewhat  closer  in  weather 
to  April  1. 

The  use  of  same-day-of -the-week  Equivalent  Days  makes  the 
choice  of  a method  for  finding  a unique  Morning  or  Evening  Equivalent 
Day  less  critical.  It  is  not  usually  possible  to  find  two  Equivalent 
Days  both  equally  close  to  the  desired  sunrise  (sunset)  time  and 
both  of  the  same-day-of- the-week  as  the  Standard  Time  day.  However, 
when  such  a situation  does  occur,  the  above  procedure  is  used. 

5.2.8  Influenced  Periods 

A change  from  Standard  Time  to  Daylight  Saving  Time  will  have 
its  greatest  effect  on  motor-vehicle  accidents  in  certain  parts  of 
the  day  (those  near  sunrise  and  sunset)  and  may  have  essentially  no 
effect  on  accidents  in  certain  other  parts  of  the  day  (those  near 
noon  and  midnight).  Therefore,  the  morning  and  evening  periods  for 
which  data  are  collected  and  for  which  the  comparisons  of  Standard 
Time  days  and  Equivalent  Days  are  made  can  be  defined  to  exclude 
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those  parts  of  the  day  in  which  the  change  will  not  influence 
motor  vehicle  accidents.  This  divides  the  day  into  those  parts 
(morning  and  evening)  in  which  accidents  may  be  influenced  by  DST 
and  those  parts  (midday  and  night)  in  which  accidents  are  considered 
not  to  be  influenced  by  DST. 

For  convenience  of  computation,  it  is  desirable  to  make  stand- 
ard definitions  of  the  morning  and  evening  periods  that  will  apply 
for  all  state-years.  Also,  since  the  effect  of  the  change  to  DST 
will  be  in  the  hours  near  sunrise  and  sunset,  the  morning  and 
evening  periods  must  be  chosen  so  that  all  sunrises  and  sunsets 
(both  Standard  Time  and  DST)  for  all  parts  of  the  country  for  all 
parts  of  the  year  are  included  in  the  defined  periods.  These  con- 
siderations led  to  the  choice  of  4 a.m.  to  10  a.m.  as  the  Morning 
period,  and  4 p.m.  to  10  p.m.  as  the  Evening  period.  This  splits 
each  twenty-four  hours  into  four  quarter-days: 


4 a.m.  to  10  a.m. : 

10  a.m.  to  4p.m.: 

4 p.m.  to  10  p.m.: 

10  p.m.  to  4 a.m.: 


Morning  (DST- Inf luenced) 
Midday  (DST-Uninfluenced) 
Evening  (DST- Inf luenced) 
Night  (DST-Uninfluenced) 


Therefore,  when  the  Equivalent  Day  approach  requires  a comparison 
of  Morning  or  Evening  accidents,  the  number  of  accidents  for 
the  particular  six  hour  period,  as  just  defined,  is  used. 


5.3  ISOLATING  THE  EFFECTS  OF  DAYLIGHT  SAVING  TIME 

To  find  the  effect  of  DST  on  motor  vehicle  accidents,  the  Equi- 
valent Day  approach  compares  the  number  of  accidents  on  several 
pairs  of  days,  each  pair  consisting  of  a Standard  Time  day  and 
its  corresponding  Daylight  Saving  Equivalent  Day.  The  discovered 
difference  in  accidents  will  be  due  partly  to  Daylight  Saving  Time, 
but  will  also  be  due  to  any  differences  in  other  acc ident -affect ing 
factors  (such  as  weather)  that  might  exist  between  the  set  of 
Standard  Time  days  and  the  set  of  Equivalent  Days.  Thus,  to  find 
the  impact  of  DST,  its  effect  must  be  isolated  from  the  effects 
of  other  accident -related  factors.  There  are  two  methods  of  treat- 
ing this  situation:  one  is  to  make  an  assumption  of  randomization; 

the  other  is  to  use  a normalization  technique. 
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5.3.1  Assumption  of  Randomization 


The  simplest  way  to  treat  the  problem  of  isolating  the  effects 
of  DST  is  to  use  a sufficiently  large  set  of  data  so  that  the  as- 
sumption can  be  made  that  the  random  effects  of  non-DST  accident 
factors  are  balanced  out.  Accident  data  were  available  for  20  state- 
years.  For  each  state-year  there  were  about  60  Standard  Time  days 
with  Equivalent  Mornings,  for  a total  of  1119  Morning -Equ ivalent 
Morning  pairs.  Also,  for  each  state-year  there  were  about  160  days 
with  Equivalent  Evenings,  for  a total  of  3192  Evening -Equivalent 
Evening  pairs.  Let  us  consider  the  assumption  that  any  effect  of 
non-DST  accident  factors  on  the  1119  Standard  Time  mornings  and 
the  3192  Standard  Time  evenings  was  approximately  balanced  by  a 
similar  effect  upon  the  1119  Equivalent  Mornings  and  the  3192 
Equivalent  Evenings.  For  accident  factors  that  are  truly  random, 
this  is  a good  assumption.  However,  since  the  Equivalent  Days  and 
Standard  Time  days  often  come  from  different  parts  of  the  year 
(see  Figures  5-4  and  5-5),  there  may  be  a general  bias  which  does 
not  balance  out  over  the  large  number  of  comparisons.  If  so,  this 
will  produce  some  degree  of  error  in  the  results  under  the  Random- 
ization Assumption.  This  possibility  of  error  leads  to  the  use  of 
the  Normalization  technique  discussed  in  the  next  section. 

The  Randomization  Assumption,  however,  does  produce  a direct 
comparison  between  Standard  Time  days  and  Equivalent  Days,  and 
thus  does  not  introduce  any  additional  assumptions  in  the  attempt 
to  isolate  the  DST  effect.  Therefore,  there  is  some  advantage  to 
using  this  assumption.  For  this  reason,  and  since  it  required 
little  added  computation,  results  based  upon  the  Randomization  As- 
sumption method  were  computed  in  addition  to  those  based  upon  the 
Normalization  technique. 

5.3.2  Normalization 

Normalization  is  a procedure  designed  to  isolate  the  effects 
of  Daylight  Saving  Time  from  other  effects,  such  as  weather,  which 
may  influence  the  number  of  accidents.  This  procedure  is  applied  to 
the  pairs  of  Standard  Time  and  Daylight  Saving  Equivalent  Days. 

The  normalization  procedure  separates  each  day  into  two  parts:  the 


29 


part  of  the  day  potentially  influenced  by  Daylight  Saving  Time 
(i.e.,  morning  and  evening),  and  the  part  of  the  day  uninfluenced 
by  Daylight  Saving  Time  (i.e,  midday  and  night).  Then  for  the  two 
days  in  each  pair  the  parts  of  the  days  influenced  by  Daylight 
Saving  Time  are  compared,  and  the  parts  of  the  days  uninfluenced 
by  Daylight  Saving  Time  are  compared.  Any  change  in  the  number  of 
accidents  in  the  hours  influenced  by  Daylight  Saving  Time  over  and 
above  the  change  in  the  uninfluenced  hours  is  assumed  to  represent 
the  effect  of  Daylight  Saving  Time. 

A step  by  step  description  of  the  normalization  methodology  is 
presented  below,  followed  by  an  example. 


Step  1:  Assumption  of  Isolated  Day-Parts 

Assume  that  for  each  day  two  parts  can  be  isolated: 

Part  I = DST - Inf luenced  parts  of  the  day 
(i.e.,  4-10  a.m.  and  4-10  p.m.) 

Part  U = DST-Uninf luenced  parts  of  the  day 

(i.e.,  10  a.m.  - 4 p.m.  and  10  p.m.  - 4 a.m.) 


Step  2:  Values  of  Impact  Parameters 

Define  four  accident  totals  for  each  state-year: 

U-cj-j,  = Total  number  of  accidents  for  the  uninfluenced  hours 
of  all  Standard  Time  days  that  have  Equivalent  Days. 


= Total  number  of  accidents  for  the  uninfluenced  hours 
of  the  corresponding  Equivalent  Days. 


I 


ST 


I 


EQ 


= Total  number  of  accidents  for  the  influenced  hours 
of  all  Standard  Time  days  that  have  Equivalent  Days. 

= Total  number  of  accidents  for  the  influenced  hours 
of  the  corresponding  Equivalent  Days. 


The  morning  and  evening  periods  are  considered  separately. 

For  each  morning  (4  a.m.  to  10  a.m.),  the  corresponding  DST-uninflu- 
enced  hours  are  chosen  to  be  the  night  preceding  (10  p.m.  to  4 a.m.) 
and  the  midday  following  (10  a.m.  to  4 p.m.).  For  each  evening 
(4  p.m.  to  10  p.m.),  the  corresponding  DST-uninf luenced  hours  are 
chosen  to  be  the  midday  preceding  (10  a.m.  to  4 p.m.)  and  the 
night  following  (10  p.m.  to  4 a.m.). 


30 


Step  5:  Change  Factors 


Define  two  change  factors: 

Change  factor  for  uninfluenced  hours  (this  is  the  comparison 
between  parts  of  the  days  uninfluenced  by  DST)  : 

c = UeQ 
U 

Change  factor  for  influenced  hours  (the  comparison  between 
parts  of  the  days  influenced  by  DST) : 

r - :EQ 

I T 

1 ST 

Change  factors  are  similar  to  percent  changes.  For  example, 
a change  factor  of  1.03  corresponds  to  a 3%  increase  in  accidents 
when  going  from  the  Standard  Time  day  to  the  Equivalent  Day;  a 
change  factor  of  0.96  corresponds  to  a 4%  decrease. 

Step  4:  Break  Up  Cj  into  a DST  and  a Non-DST  Part 

Change  factor  for  DST  influenced  hours  = 

(Change  factor  due  to  DST)  x (Change  factor  due  to  all  other  causes) 
or 


CI  C l/DSTxC I / OTHER 

Therefore 


'I 


1 1 /DST 


I/OTHER 


Step  5:  Assumption  of  Same  Non-DST  Change  in  Both  Day-Parts 


Assume  that 
Saving  Time 
hours,  that 


the  change  due  to  all  factors  other  than  Daylight 
is  the  same  in  both  the  influenced  and  uninfluenced 
is  : 


CU  C 1/ OTHER 
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Substituting  this  into  the  previous  equation  yields: 


I/DST  ui7other  ^ 


Step  6:  Change  in  the  DST- Inf luenced  Hours  Due  to  DST 

Substituting  the  change  factor  definitions  of  Step  3 into  the 
equation  of  Step  5 yields: 


'I/DST 


is 

tst 

% 

bT 


and  the  percent  change  of  accidents  in  the  influenced  hours  due 
to  DST  can  be  found  directly  from  Cp^^g^. 

Using  the  case  of  morning  accidents  in  Michigan,  1972,  as  an 
example  of  this  Equivalent  Day  Normalization  procedure,  we 
calculate : 

Ic;T=2468  = Sum  of  morning  (4  a.m.  to  10  a.m.)  injury -caus ing  accidents 


on  all  Standard  Time  days  that  have  Morning  Equivalent 
Days . 


I =:2620  = Sum  of  morning  (4  a.m.  to  10  a.m.)  injury-causing  accidents 
on  all  the  corresponding  Morning  Equivalent  Days, 

Ugrp=10277  = Sum  of  all  injury-causing  accidents  on  the  night  pre- 
ceding (10  p.m.  to  4 a.m.)  and  the  midday  following 
(10  a.m.  to  4 p.m.)  each  of  the  above  Standard  Time  days. 

UgQ=:10725  = Sum  of  all  injury-causing  accidents  on  the  night  pre- 
ceding (10  p.m.  to  4 a.m.)  and  the  midday  following 
(10  a.m.  to  4 p.m.)  each  of  the  above  Morning  Equivalent 
Days . 


= = HIS-  = 1.062  (i.e.,  + 6.2%) 


ST 


2468 
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+ 4.4%) 


r _ UEQ  10725 

U UgT  10277 


1 . 044  (i . e . , 


r = CI  _ 1.062 

I/DST  Cy  1.044 


1.017  (i.e.,  + 1.7%) 


Note  that  the  Normalization  technique  is  based  upon  the  two 
assumptions  stated  in  Steps  1 and  5.  The  Step  1 assumption  seems 
reasonable  and  has  been  discussed  in  Section  5.2.8.  The  Step  5 
assumption  is  that  the  acc ident -related  factors  which  differ  be- 
tween the  Standard  Time  days  and  the  Equivalent  Days  differ  equally 
in  the  influenced  hours  and  in  the  uninfluenced  hours.  This  will 
be  true,  in  general,  for  phenomena  that  occur  randomly  in  dif- 
ferent parts  of  the  day,  for  phenomena  that  occur  for  more  than 
half  a day  at  a time  (i.e.,  spanning  both  influenced  and  uninflu- 
enced hours)  , and  for  phenomena  that  occur  possibly  in  certain 
parts  of  the  day,  but  with  equal  likelihood  throughout  the  year 
(i.e.,  spanning  both  Standard  Time  and  Equivalent  Days).  This 
list  seems  to  include  most  acc ident - related  factors.  However,  if 
there  are  accident -related  phenomena  which  occur  only  in  certain 
parts  of  the  year  (e.g.,  more  in  the  Standard  Time  days  than  in  the 
Equivalent  Days)  and  occur  more  in  certain  times  of  the  day  than 
in  others  (e.g.,  more  in  the  influenced  hours  than  in  the  uninflu- 
enced hours) , then  there  will  be  inaccuracies  introduced  in  the 
Normalization. 

5.4  APPLICATION  TO  DOUBLE  DAYLIGHT  SAVING  TIME  (DDST) 

The  Daylight  Saving  Equivalent  Day  Methodology  can  be  applied 
to  the  estimation  of  the  number  of  accidents  under  Year  Round 
Double  Daylight  Saving  Time.  The  method  is  the  same  as  that  used 
for  Year  Round  Daylight  Saving  Time  (as  discussed  in  Sections  5.2 
and  5.3),  with  some  exceptions  that  will  be  discussed  in  the  fol- 
lowing paragraphs. 


33 


The  adoption  of  DDST  will,  of  course,  change  the  sunrise  and 
sunset  times  of  all  the  days  of  the  year,  not  just  the  Standard 
Time  days.  The  previously  Standard  Time  days  will  have  their  sun- 
rise and  sunset  times  delayed  by  two  hours,  and  the  previously  DST 
days  will  have  their  sunrise  and  sunset  times  delayed  by  one  hour. 
Therefore,  the  Equivalent  Day  methodology  for  DDST  ideally  must 
compare  the  number  of  accidents  on  all  the  days  of  the  year  against 
the  number  of  accidents  on  all  of  their  corresponding  Equivalent 
Days . 

Composite  Double  DST  Equivalent  Days  are  found  in  a similar 
manner  to  that  of  composite  Daylight  Saving  Equivalent  Days.  An 
example  is  shown  in  Figure  5-10.  The  Morning  DDST  Equivalent 
Days  corresponding  to  each  day  of  the  year  are  shown  in  Figure  5-11. 
A Morning  DDST  Equivalent  Day  exists  for  (depending  on  latitude) 
about  one-third  of  the  days  of  the  year,  approximately  130  days. 

The  Evening  DDST  Equivalent  Days  corresponding  to  each  day  of  the 
year  are  shown  in  Figure  5-12.  An  Evening  DDST  Equivalent  Day 
exists  for  (depending  on  latitude)  a little  over  half  of  the  days 
of  the  year,  approximately  200  days.  Note  that  the  days  which 
have  DDST  Equivalent  Days  are  not  equally  distributed  throughout 
the  year.  For  the  morning,  they  are  generally  April  through 
September  days;  for  the  evening,  they  they  are  generally  September 
through  April  days.  The  Equivalent  Day  approach,  therefore,  found 
the  change  due  to  having  DDST  in  these  periods  of  the  year. 
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SUNRISE 
LAT.  40 


Figure  5-10.  Example  of  Finding  a Morning  Double 
DST  Equivalent  Day 
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POSSIBLE  MORNING  DOUBLE  DST  EQUIVALENT  DAYS 


JAN,  r- 


DEC,  - 


NOV, 


OCT, 


SEP. 


AUG. 


JUL. 


JUN. 


MAY 


APR. 


MAR. 

FEB, 

JAN. 


JAN.  FEB.  MAR.  APR.  MAY  JUN.  JUL.  AUG.  SEP.  OCT.  NOV, 

DAYS  OF  THE  YEAR 


Figure  5-11.  Morning  DDST  Equivalent  Days  Corresponding 
to  the  Days  of  the  Year 


I 

DEC  . 
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POSSIBLE  EVENING  DOUBLE  DST  EQU I VALENT  DAYS 


JAN. 


I I I I I I I I I I I 

JAN.  FEB.  MAR.  APR,  MAY  JUN.  JUL.  AUG.  SEP.  OCT . NOV.  DEC. 

DAYS  OF  THE  YEAR 


Figure  5-12.  Evening  DDST  Equivalent  Days  Corresponding 
to  the  Days  of  the  Year 
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6,  DESCRIPTION  OF  THE  LATEDAY  PROGRAMS 


There  are  five  major  LATEDAY  Programs.  Three  of  these 
programs  are  used  in  the  determination  of  the  level  of  accidents 
under  YRDST: 

1.  The  EQDAY  program  (discussed  in  Section  6.1)  computes  the 
Daylight  Saving  Equivalent  Days  for  a given  state  for  a 
given  year.  This  program  accepts  as  input  a sunrise/sunset 
table  for  the  latitude  of  the  state's  geographic  center 
and  a list  of  the  holidays  for  the  state  for  the  given 
year.  It  produces  a list  of  the  Morning  and  Evening 
Equivalent  Days  for  each  non-D.S.T.  day  of  the  year.  " 

2.  The  QDY  program  (Section  6.2)  reads  in  a LATEDAY  Tape, 
extracts  pertinent  accident  data,  and  produces  an  array 
showing  the  number  of  accidents  involving  at  least  one 
fatality  or  injury,  for  each  quarter-day  of  the  year.  This 
program  was  optimized  to  minimize  the  time  of  execution, 
since  it  was  used  to  process  data  on  several  million 
accidents . 

3.  The  COMPAR  program  (Section  6.3)  reads  in  both  the  computed 
table  of  Daylight  Saving  Equivalent  Days,  and  the  accident 
array  discussed  above,  and  then  computes,  tabulates,  and 
compares  the  accidents  for  each  non-D.S.T.  day  and  its 
Equivalent  Day. 

To  determine  the  level  of  accidents  under  Year-Round  Double  DST, 
three  programs  (including  one  of  the  above)  are  used: 

1.  The  DBLDAY  program  (Section  6.4)  computes  the  Double  Day- 
light Saving  Equivalent  Days  for  a given  state  for  a given 
year.  DBLDAY,  as  EQDAY  above,  accepts  a sunrise/sunset 
table  and  a list  of  holidays  as  input.  It  produces  a list 
of  Morning  and  Evening  DDST  Equivalent  Days  for  each  day 
of  the  year. 
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2. 


The  QDY  program,  mentioned  above,  finds  the  number  of 
accidents  for  each  quarter-day  of  the  year. 

3.  The  DBLCOMP  program  (Section  6.5)  reads  in  both  the 

computed  table  of  DDST  Equivalent  Days  and  the  accident 
array,  and  then  computes,  tabulates,  and  compares  the 
accidents  for  each  day  of  the  year  and  its  DDST  Equivalent 
Day. 

The  five  LATEDAY  programs  were  written  in  FORTRAN  IV  on  an 
IBM  370/155  computer.  Listings  of  the  programs  appear  in 
Appendix  A. 

6.1  EQUIVALENT  DAY  (EQDAY)  PROGRAM 

The  Equivalent  Day  (EQDAY)  Program  computes  the  Daylight 
Saving  Equivalent  Days  for  a state-year.  Based  upon  the  latitude 
of  the  geographic  center  of  the  state,  a set  of  sunrise  and  sunset 
times  are  found  in  the  Tables  of  Sunrise,  Sunset,  and  Twilight 
(Supplement  to  the  American  Ephemeris,  1946),  published  by  the 
Naval  Almanac  Office,  U.S.  Naval  Observatory,  Washington,  D.C. 

These  times  are  keypunched  and  then  entered  into  a computer  file 
for  use  as  input  to  EQDAY.  In  addition,  state  governments  are 
contacted  to  determine  the  legal  holidays  for  each  state  for  the 
data  year.  These  are  coded  (January  1=1,  December  31  = 366), 
keypunched,  and  entered  into  the  input  computer  file  for  EQDAY. 

Using  the  above  inputs,  the  EQDAY  program  (see  listing, 
Appendix  A.l)  is  run.  An  outline  of  the  operation  of  the  EQDAY 
program  follows: 

1.  The  state  name,  data  year  name,  list  of  holidays  and 
lists  of  sunrise/sunset  times  are  read  in. 

2.  An  hour  is  added  to  the  appropriate  sunrise/ sunset  times 
to  account  for  DST  in  the  data  year. 

3.  February  29  is  deleted  from  the  sunrise/sunset  tables  for 
non-leap  years. 
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4. 


For  each  Standard  Time  day  of  the  year,  the  best  same-day- 
of-the-week  Equivalent  Day  is  found.  Only  days  which  have 
a sunrise  (sunset)  60  ^+5  minutes  later  than  the  sunrise 
(sunset)  of  the  Standard  Time  day  are  considered.  Also, 
no  day  which  has  already  been  a Morning  (Evening)  Equi- 
valent Day  4 times  is  considered.  The  day  closest  to  a 
60  minute  difference  is  used.  If  2 or  more  days  are 
equally  close  to  a 60  minute  difference,  the  Weather 
Function  described  in  Section  5.2.7  is  used. 

5.  A list  of  the  Morning  and  Evening  Equivalent  Days  cor- 
responding to  each  Standard  Time  day  is  stored  in  a com- 
puter file  for  use  as  input  to  the  COMPAR  program. 


6.2  QUARTER  DAY  (QDY)  PROGRAM 

The  Quarter  Day  (QDY)  Program  (see  listing.  Appendix  A. 2)  ex- 
tracts accident  data  from  a LATEDAY  tape  for  each  quarter  day  of 
the  year. 

The  input  to  QDY  is  the  LATEDAY  tape  which  has  the  data  on 
accidents  for  the  state-year. 


Using  this  input,  the  QDY  program  operates  as  follows: 

1.  Each  accident  entry  is  read  from  the  LATEDAY  tape  using 
the  format  shown  in  Table  4-2. 


2.  The  data  are  converted  into  numerical  form  from  the  byte 
codes  of  Table  4-2. 

3.  The  quarter  day  of  the  accident  is  found  based  on  the 
time  of  the  accident: 


000 

to 

359  = Quarter 

day 

4 

400 

to 

1059  = Quarter 

day 

1 

1100 

to 

1559  = Quarter 

day 

2 

1600 

to 

2159  = Quarter 

day 

3 

2200 

to 

2359  = Quarter 

day 

4 

of  previous  day 


4.  Two  matrices  are  filled,  a fatality  matrix  and  an  injury- 
causing  accidents  matrix.  These  matrices  have  366  x 4 
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entries,  one  for  each  quarter  day  of  the  year.  If  the  ac- 
cident involved  one  or  more  fatalities,  the  corresponding 
quarter  day's  entry  in  the  fatality  matrix  is  increased 
by  one.  Similarly,  if  the  accident  involved  one  or  more 
injuries  or  fatalities  (i.e.,  injuries  of  severity  K,  A, 

B,  or  C) , the  corresponding  quarter  day's  entry  in  the 
in j ury- caus ing  accident  matrix  is  increased  by  one. 

5.  A total  count  is  kept  of  all  fatalities  and  all  non-in- 
jury accidents. 

6.  A matrix  of  the  hourly  inj ury- causing  accidents  for  the 
Evening  quarter  day  is  also  found. 

7.  The  matrices  and  totals  found  are  printed  out  and  are  stored 
in  a computer  file  for  use  as  input  to  the  COMPAR  and 
DBLCOMP  programs.  The  fatality  matrix  and  the  fatality 

and  non-injury  accident  totals  are  found  and  printed  for 
comparison  with  the  injury-causing  accident  matrix,  but 
are  not  now  used  in  the  COMPAR  or  DBLCOMP  programs. 

6.3  THE  COMPARISON  (COMPAR)  PROGRAM 

The  Comparison  (COMPAR)  Program  (see  listing,  Appendix  A. 3) 
computes  and  tabulates  a comparison  between  the  number  of  accidents 
on  Standard  Time  days  and  the  number  of  accidents  on  their  cor- 
responding Morning  and  Evening  Equivalent  Days.  Based  upon  this, 
the  program  produces  an  estimate  of  the  level  of  accidents  under  Year 
Round  Daylight  Saving  Time. 

The  inputs  to  COMPAR  are  the  data  files  produced  by  the  EQDAY 
and  QDY  programs. 

Using  these  inputs,  the  COMPAR  program  operates  as  follows  (a 
sample  printout  of  this  program's  output  is  shown  in  Appendix  B.l): 

1.  The  accident  matrices  produced  by  QDY  are  read  in.  The 
fatality  and  injury-causing  accident  matrices  are  printed 
out . 

2.  For  information  and  comparison  purposes,  four  tables  are 

computed  and  printed  out:  the  accident  totals  by  day-of- 
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the-week  and  quarter  day,  the  accident  totals  by  month 
and  quarter  day,  the  average  accidents  by  day-of- the -week 
and  quarter  day,  and  the  average  accidents  by  month  and 
quarter  day. 

3.  The  Equivalent  Days  found  by  EQDAY  are  read  in  and  printed 
out . 

4.  The  accident  totals  in  the  Morning  quarter  day  and  in  the 
Evening  quarter  day  are  found  for  the  DST  and  the  non-DST 
parts  of  the  year. 

5.  The  actual  accidents  for  each  Standard  Time  day  of  the 
year  are  compared  with  the  accidents  of  their  Equivalent 
Days.  A large  table  is  printed  out.  For  each  Standard 
Time  day,  the  table  contains: 

a.  Day- of - the -week  (1  = Monday,..., 7 = Sunday) 

b.  Date 

c.  Number  of  Morning  accidents  (preceded  by  a minus  sign 
if  no  Equivalent  Morning  exists) 

d.  Number  of  Evening  accidents  (preceded  by  a minus  sign 
if  no  Equivalent  Evening  exists) 

e.  Number  of  Night-of-Previous-Day  accidents  + Number  of 
Midday  accidents  (Normalization  for  the  Morning) 

f.  Number  of  Midday  accidents  + Number  of  Night  accidents 
(Normalization  for  the  Evening) 

For  the  corresponding  Equivalent  Day,  the  table  con- 
tains : 

a.  Day-of- the-week  of  the  Equivalent  Morning 

b.  Date  of  the  Equivalent  Morning 

c.  Number  of  Equivalent  Morning  accidents  (0  = no  Equi- 
valent Morning) 

d.  Day-of- the-week  of  the  Equivalent  Evening 

e.  Date  of  the  Equivalent  Evening 
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£.  Number  of  Equivalent  Evening  accidents  (0  = No  Equi- 
valent Evening) 

g.  Number  of  Night  - of - Previous  - Day  accidents  + Number  of 
Midday  accidents  for  the  Equivalent  Morning  (Equiva- 
lent Morning  Normalization) 

h.  Number  of  Midday  accidents  + Number  of  Night  accidents 
for  the  Equivalent  Evening  (Equivalent  Evening  Normal- 
ization) 

The  overall  totals  of  all  of  the  above  accident  entries 
are  found  and  printed  out  at  the  bottom  of  the  table. 

For  the  Standard  Time  days,  only  the  entries  for  which  a 
corresponding  Equivalent  Morning  or  Evening  exists  (i.e., 
those  not  preceded  by  minus  signs)  are  totaled. 

6.  The  Morning  Apparent  Percent  Change  and  the  Evening  Ap- 
parent Percent  Change  are  found  and  printed.  Each  of  these 
is  the  percent  difference  between  the  actual  total  number 
of  accidents  on  the  Standard  Time  days  and  the  correspond- 
ing total  number  of  accidents  on  the  Equivalent  Days.  If 
the  Randomization  Assumption  is  accepted,  these  Apparent 
Percent  Changes  indicate  the  change  in  the  level  of  Morn- 
ing and  Evening  accidents  under  the  effect  of  YRDST. 

7.  The  percent  change  for  the  Night  + Midday  accidents  (Morn- 
ing Normalization)  and  for  the  Midday  + Night  accidents 
(Evening  Normalization)  are  found. 

8.  Using  these,  the  Morning  Net  Percent  Change  and  the  Even- 
ing Net  Percent  Change  are  found  and  printed.  These  are 
the  results  of  the  Normalization  process  and  indicate  the 
change  in  the  level  of  Morning  and  Evening  accidents  under 
the  effect  of  YRDST. 

9.  The  Apparent  and  Net  Percent  changes  refer  to  the  change 
in  the  accident  totals  for  the  Morning  or  Evening  quarter 
day  for  the  six  months  which  had  Standard  Time  in  the  data 
years.  To  get  an  estimate  of  the  effect  of  YRDST  on  the 
yearly  accident  totals,  an  accident  total  for  a year  with 
YRDST  was  computed  as  follows: 
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Yearly  accidents  under  YRDST  = 


(Year's  DST  Accidents) 

+(100%  + Year's  Morning  % Change)  x (Year's  Non-DST 
Morning  Accidents) 

+ (Year's  Non-DST  Midday  Accidents) 

+ ( 100%  + Year's  Evening  % Change)  x (Year's  Non-DST 
Evening  accidents) 

+(Year's  Non-DST  Night  Accidents) 

This  is  computed  in  two  different  ways:  first,  by  using 

the  Apparent  Changes,  and  second,  by  using  the  Net  Percent 
Changes.  These  are  then  compared  with  the  actual  data 
year's  total  accidents. 

6.4  DOUBLE  DST  EQUIVALENT  DAY  (DBLDAY)  PROGRAM 

The  Double  DST  Equivalent  Day  (DBLDAY)  Program  computes  the 
Double  Daylight  Saving  Equivalent  Days  for  a state-year.  DBLDAY 
is  similar  to  EQDAY  and  accepts  the  same  inputs,  a set  of  sunrise/ 
sunset  times  and  a list  of  state  holidays.  Using  the  above  inputs, 
the  DBLDAY  program  (see  listing,  Appendix  A. 4)  is  run.  An  outline 
of  the  operation  of  the  DBLDAY  program  follows: 

1.  The  state  name,  the  data  year  name,  the  list  of  holidays, 
and  the  lists  of  sunrise/sunset  times  are  read  in. 

2.  An  hour  is  added  to  the  appropriate  sunrise/sunset  times 
to  account  for  DST  in  the  data  year. 

3.  February  29  is  deleted  from  the  sunrise/sunset  tables  for 
non-leap  years. 

4.  For  each  day  of  the  year,  the  best  same-day-of-the-week 
Equivalent  Day  is  found.  Only  days  which  have  a sunrise 
(sunset)  120  ±5  minutes  later  than  the  sunrise  (sunset) 
of  a Standard  Time  day  or  60  ±5  minutes  later  than  the 
sunrise  (sunset)  on  a DST  day  are  considered.  Also,  no 
day  which  has  already  been  a Morning  (Evening)  Equivalent 
Day  four  times  is  considered.  The  day  closest  to  the  de- 
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sired  time  difference  is  used.  If  two  or  more  days  are 
equally  close  to  the  desired  time  difference,  the  Weather 
Function  described  in  Section  5.2.7  is  used. 

5.  A list  of  the  Morning  and  Evening  Equivalent  Days  cor- 
responding to  each  Standard  Time  day  is  output  to  a com- 
puter file  for  use  as  input  to  the  DBLCOMP  program. 

6.5  THE  DOUBLE  DST  COMPARISON  (DBLCOMP)  PROGRAM 

The  Double  DST  Comparison  (DBLCOMP)  Program  (see  listing,  Ap- 
pendix A. 5)  computes  and  tabulates  a comparison  between  the  number 
of  accidents  on  Standard  Time  days  and  the  number  of  accidents  on 
their  corresponding  Double  DST  Morning  and  Evening  Equivalent  Days. 
Based  upon  this,  the  program  produces  an  estimate  of  the  level  of 
accidents  under  Double  Daylight  Saving  Time.  The  inputs  to  DBLCOMP 
are  the  data  files  produced  by  the  DBLDAY  and  QDY  programs. 

Using  these  inputs,  the  DBLCOMP  program,  in  a similar  manner 
to  the  COMPAR  program,  operates  as  follows  (a  sample  printout  of 
this  program's  output  is  shown  in  Appendix  B.2): 

1.  The  accident  matrices  produced  by  QDY  are  read  in.  The 
fatality  and  inj ury-caus ing  accidents  matrices  are  printed 
out . 

2.  For  information  and  comparison  purposes,  four  tables  are 

computed  and  printed  out:  the  accident  totals  by  day-of- 

the-week  and  quarter  day,  the  accident  totals  by  month 
and  quarter  day,  the  average  accidents  by  day-of-the-week 
and  quarter  day,  and  the  average  accidents  by  month  and 
quarter  day. 

3.  The  Double  DST  Equivalent  Days  found  by  DBLDAY  are  read 
in  and  printed  out. 

4.  The  accident  totals  in  the  Morning  quarter  day  and  in  the 
Evening  quarter  day  are  found  for  the  DST  and  the  non- DST 
parts  of  the  year. 

5.  The  actual  accidents  for  each  day  of  the  year  are  com- 
pared with  the  accidents  of  their  DDST  Equivalent  Days. 
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A large  table  is  printed  out.  For  each  day  of  the  year, 

the  table  contains: 

a.  Day- of - the -week  (1  = Monday,..., 7 = Sunday) 

b.  Date 

c.  Number  of  Morning  accidents  (preceded  by  a minus  sign 
if  no  DDST  Equivalent  Morning  exists) 

d.  Number  of  Evening  accidents  (preceded  by  a minus  sign 
if  no  DDST  Equivalent  Evening  exists) 

e.  Number  of  Night-of - Previous -Day  accidents  + Number  of 
Midday  accidents  (Normalization  for  the  Morning) 

f.  Number  of  Midday  accidents  + Number  of  Night  accidents 
(Normalization  for  the  Evening) 

For  the  corresponding  DDST  Equivalent  Day,  the  table 

contains : 

a.  Day-of - the-week  of  the  DDST  Equivalent  Morning 

b.  Date  of  the  DDST  Equivalent  Morning 

c.  Number  of  DDST  Equivalent  Morning  accidents  (0  = no 
DDST  Equivalent  Morning) 

d.  Day-of -the-week  of  the  DDST  Equivalent  Evening 

e.  Date  of  the  DDST  Equivalent  Evening 

f.  Number  of  DDST  Equivalent  Evening  accidents  (0  = no 
DDST  Equivalent  Evening) 

g.  Number  of  Night  - of - Previous  - Day  accidents  + Number  of 
Midday  accidents  for  the  DDST  Equivalent  Morning  (DDST 
Equivalent  Morning  Normalization) 

h.  Number  of  Midday  accidents  + Number  of  Night  accidents 
for  the  DDST  Equivalent  Evening  (DDST  Equivalent  Even- 
ing Normalization) 

The  overall  totals  of  all  of  the  above  accident  entries 

are  found  and  printed  out  at  the  bottom  of  the  table. 

For  the  actual  days  of  the  year,  only  the  entries  for 
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which  a corresponding  DDST  Equivalent  Morning  or  Evening 
exists  (i.e.,  those  not  preceded  by  minus  signs)  are  to- 
taled. 

6.  The  Morning  DDST  Apparent  Percent  Change  and  the  Evening 
DDST  Apparent  Percent  Change  are  found  and  printed.  Each 
of  these  is  the  percent  difference  between  the  actual  to- 
tal number  of  accidents  on  the  days  of  the  year  and  the 
corresponding  total  number  of  accidents  on  the  DDST  Equi- 
valent Days.  If  the  Randomization  Assumption  is  accepted, 
these  Apparent  Percent  Changes  indicate  the  change  in  the 
level  of  Morning  and  Evening  accidents  under  the  effect 

of  DDST. 

7.  The  percent  change  for  the  Night  + Midday  accidents  (Morn- 
ing Normalization)  and  the  percent  change  for  the  Midday 

+ Night  accidents  (Evening  Normalization)  are  found. 

8.  Using  these,  the  Morning  DDST  Net  Percent  Change  and  the 
Evening  DDST  Net  Percent  Change  are  found  and  printed. 
These  are  the  results  of  the  Normalization  process  and 
indicate  the  change  in  the  level  of  Morning  and  Evening 
accidents  under  the  effect  of  DDST. 

9.  All  totals  for  Morning,  Evening,  Night  + Midday,  and  Mid- 

day + Night  accidents  for  both  actual  days  of  the  year  and 
DDST  Equivalent  Days  are  broken  up  into  two  parts:  the 

total  corresponding  to  the  DST  period  of  the  data  year 
and  the  total  corresponding  to  the  non-DST  period  of  the 
year.  These  are  used  to  find  and  print  out  Apparent  and 
Net  Percent  Changes  for  the  DST  and  non-DST  periods  of 
the  year. 

10.  The  Apparent  and  Net  Percent  changes  refer  to  the  change 
in  the  accident  totals  for  the  Morning  or  Evening  quarter 
days  for  the  days  which  had  DDST  Equivalent  Days.  To  get 
an  estimate  of  the  effect  of  DDST  on  the  yearly  accident 
totals,  an  accident  total  for  a year  with  DDST  was  com- 
puted as  follows: 
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Yearly  accidents  under  DDST  = 


(1001  + Year's  Morning  % Change)  x (Year's  Morning 
Accidents ) 

+(Year's  Midday  Accidents) 

+(100%  + Year's  Evening  I Change)  x (Year's  Evening 
Accidents ) 

+(Year's  Night  Accidents) 

This  is  found  using  the  Apparent  Percent  Changes  and  using 
the  Net  Percent  Changes,  and  these  are  compared  with  the 
actual  data  year's  total  accidents. 

11.  The  Yearly  comparisons  above  are  also  made  using  the  DST 
and  non-DST  period  data.  In  this  case,  the  accident  to- 
tal for  a year  with  DDST  was  computed  as  follows: 

Yearly  accidents  under  DDST  = 

(100%  + Year's  Morning  DST  % Change)  x (Year's  DST 
Morning  Accidents) 

+(100%  + Year's  Morning  Non-DST  % Change)  x (Year's 
Non-DST  Morning  Accidents) 

+(Year's  Midday  Accidents) 

+(100%  + Year's  Evening  DST  % Change)  x (Year's 
DST  Evening  Accidents) 

+(100%  + Year's  Evening  Non-DST  % Change)  x (Year's 
Non-DST  Evening  Accidents) 

+(Year's  Night  Accidents) 

This  is  computed  in  two  different  ways:  first,  by  using 

the  Apparent  Changes,  and  second,  by  using  the  Net  Per- 
cent changes.  These  are  then  compared  with  the  actual 
data  year's  total  accidents. 
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7.  RESULTS 


7.1  YRDST  RESULTS 

The  LATEDAY  YRDST  Programs  (EQDAY , QDAY  AND  COMPAR)  were  run 
for  each  of  the  20  state-years  for  which  data  were  available. 

The  results  are  shown  in  Tables  7-1  to  7-6. 

Tables  7-1  and  7-2  show,  for  the  non-DST  period  of  the  year, 
the  estimated  change  in  in j ury- caus ing  accidents  that  would  have 
occurred  in  the  morning  (4  a.m.  to  10  a.m.)  and  the  evening  (4  p.m. 
to  10  p.m.)  parts  of  the  day  if  there  had  been  YRDST.  For  each 
state  and  year  for  which  data  were  available.  Table  7-1  contains: 

1.  The  number  of  non-DST  days  that  had  Equivalent  mornings. 

2.  The  total  number  of  actual  morning  accidents  on  the  days 
that  had  Equivalent  mornings. 

3.  The  corresponding  total  number  of  accidents  on  the  Equi- 
valent mornings. 

4.  The  difference  between  these  two  totals. 

5.  The  apparent  percent  change  between  these  two  totals. 

6.  The  morning  normal i zat ion ^ i.e.,  the  change  in  accidents 
in  the  "Night  + Midday"  period. 

7.  The  net  percent  change  (found  as  shown  in  Section  5.3.2). 

The  columns  are  totaled  and  the  overall  morning  percent  changes 
are  shown. 

Table  7-2  shows  the  evening  results  in  a form  similar  to  the 
morning  results.  For  each  state  and  year,  Table  7-2  contains: 

1.  The  number  of  non-DST  days  that  had  Equivalent  evenings. 

2.  The  total  number  of  actual  evening  accidents  on  the  days 
that  had  Equivalent  evenings. 

3.  The  corresponding  total  number  of  accidents  on  the  Equi- 
valent evenings. 

4.  The  difference  between  these  two  totals. 
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TABLE  7-1.  MORNING  (4  A.M.  TO  10  A.M.)  INJURY - CAUS ING  ACCIDENTS 
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TABLE  7-2.  EVENING  (4  P.M.  TO  10  P.M.)  INJURY- CAUS ING  ACCIDENTS 
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TX72 

163 

14  036 

15  019 

333 

7 . ON 

1 . ON 

=i  _ QV 

IJT  71 

153 

2150 

2243 

9 3 

4 . 3N 

1 . 2N 

3 . ON 

VT7 1 

139 

57 1 

63  0 

59 

1 0 . 3". 

8 .4N 

1.0  \ 

WR72 

148 

5754 

cr  -1  cr 

-519 

-9  . ON 

-3 . 3N 

-0 .7N 

TOTRl 

3 1 92 

1 02372 

1 054  04 

3 032 

2 . 96 N 

0 . 04 N 

2 . 92 N 
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TABLE  7-3.  NIGHT  AND  MIDDAY  INJURY- CAUSING  ACCIDENTS 
(MORNING  NORMALIZATION) 


State  ^ Year 

# Eqvl . Days 

Actual 
Acc.  Totals 

Eqvl.  Day 
Acc.  Totals 

Change 

^Change 

CA  72 

53 

13032 

12627 

-405 

-3.1 

CT  72 

59 

1732 

1997 

265 

15.3 

DC  71 

51 

399 

456 

57 

14.3 

DC  72 

55 

466 

461 

-5 

-1.1 

DE  72 

54 

280 

303 

23 

8.2 

GA  72 

53 

1488 

1520 

32 

2 . 2 

KS  72 

52 

1386 

1424 

38 

2.7 

ME  71 

60 

523 

651 

128 

24.5 

MI  72 

69 

10277 

10725 

448 

4.4 

MO  72 

52 

531 

549 

18 

3.4 

NH  72 

65 

289 

389 

100 

34.6 

NM  72 

51 

603 

703 

100 

16.6 

OK  72 

52 

1031 

1185 

154 

14.9 

PA  72 

58 

4478 

5294 

816 

18.2 

SC  71 

48 

882 

790 

-32 

-3.9 

SC  72 

51 

934 

1011 

77 

8.2 

TX  72 

50 

5611 

5823 

212 

3.8 

UT  71 

51 

689 

883 

194 

25.8 

VT  71 

60 

345 

434 

89 

28.2 

WA  72 

75 

4226 

4286 

60 

1.4 

TOTAL 

1119 

49142 

51511 

2369 

4.82% 
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TABLE  7-4.  MIDDAY  AND  NIGHT  INJURY  - CAUSING  ACCIDENTS 
(EVENING  NORMALIZATION) 


State 

§ Year 

# Eqvl.  Days 

Actual 
Acc.  Totals 

Eqvl.  Day 
Acc.  Totals 

Change 

^Change 

CA 

72 

162 

38623 

39304 

681 

1 . 8 

CT 

72 

154 

4760 

4625 

-135 

-2.8 

DC 

71 

158 

1294 

1367 

73 

5.6 

DC 

72 

169 

1436 

1419 

-17 

-1.2 

DE 

72 

163 

804 

897 

93 

11.6 

GA 

72 

168 

4531 

4661 

130 

2.9 

KS 

72 

167 

4787 

4398 

-389 

8.1 

ME 

71 

150 

1569 

1435 

-134 

-8.5 

MI 

72 

161 

26153 

25618 

-535 

-2.0 

MO 

72 

164 

1605 

1762 

157 

9.8 

NH 

72 

158 

944 

1002 

58 

6.1 

NM 

72 

169 

2027 

2119 

92 

4.5 

OK 

72 

166 

3269 

3485 

216 

6 . 6 

PA 

72 

152 

12240 

12314 

74 

0.6 

SC 

71 

161 

2532 

2697 

165 

6. 5 

SC 

72 

167 

3013 

2943 

-70 

-2 . 3 

TX 

72 

163 

18124 

18312 

188 

1.0 

UT 

71 

153 

2409 

2439 

30 

00 

VT 

71 

139 

909 

985 

76 

1 . 2 

WA 

72 

148 

8350 

7653 

-697 

-8.3 

TOTAL 

3192 

139379 

139435 

56 

0.04% 
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TABLE  7-5.  IMPACT  OF  YRDST  ON  YEARLY  ACCIDENTS 
(USING  APPARENT  PERCENT  CHANGE) 


STATE 

.} 

YEAR 

YEAR  WITH 
SIX-MONTH 
D . 3 . T . 

■ACTUAL > 

YEAR  WITH 
FULL -YEAR 

D . 3 . T . 

< CALCULATED  > 

CHANGE 

PERCENT 

CHANGE 

CATE' 

1 8 0849 

133353 

£5  04 

1 . 3 3 

CT7£ 

£3511 

£37  09 

1 93 

0 .34’: 

DC?  1 

6745 

697  0 

££5 

3 .34"; 

DC7£ 

68 1 1 

69  33 

1 33 

1 .79?. 

DE7£ 

3967 

3913 

-54 

- 1 . 36X 

GA73 

£3139 

c!  >’•  6 c!  7 

433 

l .39*: 

K.  o r i_ 

£ 1 £6 1 

£1399 

133 

0 .65’: 

ME?  1 

7710 

7347 

137 

l .73*: 

MI7£ 

116711 

1 1 6664 

-47 

-0 . 04*. 

MQ7£ 

8419 

3 734 

3 05 

3 . 63"- 

MH7£ 

4837 

5 1 5 0 

313 

6 .47’: 

liM7£ 

1 0595 

1 1 04  0 

445 

4 .£0’: 

OK  73 

1644£ 

16301 

359 

£ .13*. 

PA7£ 

65548 

66549 

1 0 0 1 

l .53’: 

SC  71 

1 3£64 

13463 

£04 

i .54’-: 

SC?£ 

14£01 

14474 

£73 

1 .93'. 

TX73 

86701 

38314 

1613 

l .36*: 

UT7 1 

13538 

136£6 

38 

0 .£37; 

VT7 1 

493  0 

5107 

177 

3 . 59". 

UJA7£ 

41635 

4131  0 

-3  £5 

-0 .73". 

TOTAL 

670914 

678973 

3064 

i . £ O': 

TABLE  7-6.  IMPACT  OF  YRDST  ON  YEARLY  ACCIDENTS 
(USING  NET  PERCENT  CHANGE) 


STATE 

•> 

Cx 

YEAR 

YEAR  WITH 
SIX-MONTH 

D .3  .T  . 

< ACTUAL) 

YEAR  WITH 
FULL— YEAR 

d . : . t . 

< CALCULATED  ) 

CHANGE 

PERCENT 

CHANGE 

CA7£ 

130349 

1 333 03 

£354 

1 .30’. 

CT7£ 

£351  1 

£3613 

1 01 

0.4  3*. 

DC?  1 

6745 

63  3 0 

35 

1 .£6’. 

DC7£ 

63 1 1 

6954 

143 

£ . 1 O’. 

DE7£ 

3967 

1 | 

-140 

7*  ^ -.v 

GA7£ 

£3139 

£3463 

1 . 13’. 

KS7£ 

£ 1 £6 1 

£1659 

393 

1 .37’. 

ME?  1 

7710 

7 P,  IS 

135 

1 .63". 

MI7E 

1 16711 

1 16745 

34 

0 . 03’. 

M07£ 

3419 

3579 

160 

1 .90’. 

NH7£ 

4337 

4 354 

1 17 

£.43’. 

NM7£ 

1 0595 

1 0303 

£13 

£.01’. 

0K73 

16443 

16430 

-13 

-0.07*. 

PA7£ 

65548 

65733 

174 

0 .£?’. 

SC71 

1 3 £6  4 

1 3363 

98 

0.74’. 

SC-73 

14301 

1443  0 

£19 

1 .54’. 

TX7£ 

367  0 1 

37888 

1 137 

1 .37*. 

IJT71 

1 3583 

13373 

-£  1 0 

-1  .55*. 

VT71 

4930 

4953 

■> 

i—  L. 

0.45’. 

WR7£ 

41635 

4 1367 

£3£ 

0.56*. 

TOTAL 

67  09 1 4 

67 6437 

5573 

0.33’. 

53 
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The  apparent  percent  change  between  these  two  totals. 

6.  The  evening  normalization,  i.e.,  the  change  in  accidents 
in  the  "Midday  + Night"  period. 

7.  The  net  percent  change  (found  as  shown  in  Section  5.3.2). 

The  columns  are  totaled  and  the  overall  evening  percent  changes 
are  shown. 

The  morning  and  evening  normalizations  used  in  Tables  7-1 
and  7-2  were  obtained  as  shown  in  Tables  7-3  and  7-4.  Table  7-3 
shows  the  morning  normalization,  using  accident  totals  from  the 
"Night  + Midday"  period  (the  10  p.m.  to  4 a.m.  period  preceding  the 
specified  morning  and  the  10  a.m.  to  4 p.m.  period  following  it). 
Table  7-4  shows  the  evening  normalization,  using  accident  totals 
from  the  "Midday  + Night"  period  (the  preceding  10  a.m.  to  4 p.m. 
period  and  the  10  p.m.  to  4 a.m.  period  following). 

The  overall  results  from  Tables  7-1  and  7-2  show  a Net  Percent 
Change  of  +4.71%  in  morning  injury-causing  accidents  for  the  non- 
DST  six  months  of  the  year  (based  upon  an  Apparent  Percent  Change 
of  +9.76%),  and  a Net  Percent  Change  of  +2.92%  in  evening  injury- 
causing  accidents  for  the  non-DST  six  months  of  the  year  (based 
upon  an  Apparent  Percent  Change  of  +2.96%).  These  results  indi- 
cate an  increase  in  accidents  under  YRDST  in  the  morning  and  a 
somewhat  smaller  increase  in  the  evening. 

Tables  7-5  and  7-6  show  the  impact  of  YRDST  on  the  yearly  in- 
jury-causing accident  totals  (under  the  assumption  that  the  acci- 
dents in  the  non-DST  morning  and  evening  periods  are  changed  by 
the  percentages  shown  in  Tables  7-1  and  7-2,  and  that  the  acci- 
dents in  all  other  parts  of  the  year  remain  unchanged).  For  each 
state-year,  the  actual  year's  accidents  are  shown  and  compared 
with  the  calculated  YRDST  accidents.  Table  7-5  uses  the  Apparent 
Percent  Change  (from  Tables  7-1  and  7-2),  and  Table  7-6  uses  the 
Net  Percent  Change  (From  Tables  7-1  and  7-2). 

The  yearly  in j ury- causing  accident  totals  show  a normalized 
Net  Percent  Change  of  +0.83%  (based  upon  an  Apparent  Percent  Change 
of  +1.20%).  This  result  indicates  an  increase  in  yearly  accidents 
of  the  order  of  one-half  to  one  percent  if  YRDST  is  introduced. 
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It  should  be  noted  that  these  results  are,  of  course,  only 
estimates  based  upon  data  from  a not  necessarily  representative 
group  or  states.  Furthermore,  the  Equivalent  Day  Methodology, 
though  possible  the  best  available  technique,  still  includes  var- 
ious inaccuracies,  as  was  noted  in  the  discussion  of  the  method- 
ology. 

Thus,  the  results  should  be  considered  only  as  an  indication 
of  a small  increase  in  morning  and  evening  accidents  yielding  a 
net  change  in  yearly  injury-causing  accidents  of  between  one-half 
and  one  percent,  and  the  calculated  percent  changes  should  not  be 
taken  as  exact  figures. 

7.2  DDST  RESULTS 

The  LATEDAY  DDST  Programs  (DBLDAY , QDY , and  DBLCOMP)  were 
run  for  each  of  the  20  state-years  for  which  data  were  avail- 
able. The  results  are  shown  in  Tables  7-7  through  7-12. 

Table  7-7  shows  the  estimated  change  in  in j ury- caus ing  acci- 
dents that  would  have  occurred  in  the  morning  (4  a.m.  to  10  a.m.) 
and  the  evening  (4  p.m.  to  10  p.m.)  parts  of  the  day  if  there 
had  been  Year  Round  Double  DST.  For  each  state  and  year  for  which 
data  were  available,  Table  7-7  contains: 

1.  The  number  of  days  of  the  year  that  had  DDST  Equivalent 
mornings . 

2.  The  total  number  of  actual  morning  accidents  on  the  days 
that  had  Equivalent  mornings. 

3.  The  corresponding  total  number  of  accidents  on  the  DDST 
Equivalent  mornings. 

4.  The  difference  between  these  two  totals. 

5.  The  apparent  percent  change  between  these  two  totals. 

6.  The  morning  normalization,  i.e.,  the  change  in  accidents 
in  the  "Night  + Midday"  period. 

7.  The  net  percent  change  (found  as  shown  in  Section  5. 5. 2). 
The  columns  are  totaled  and  the  overall  DDST  morning  percent 
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TABLE  7-7.  MORNING  (4  A.M.  TO  10  A.M.)  DDST 
INJURY-  CAUSING  ACCIDENTS 


ST  RTF. 
YERP 

MO  . 
EQVL 
PR  V S 

RCTURL 
RCC I PENT 
TOTRL 3 

EQVL  DRY 
RCC I PENT 
totrl: 

LHhMhE 

RFPRPEMT 

PERCENT 

CHRNGE 

NIGHT  + 
MIPPRY 
CHRNGE 

MET 

PERCENT 
CHRN hE 

CRTS 

1 £ 0 

7 356 

3615 

1 £59 

17  . r. 

-8.17 

£7.47 

C T 7 S' 

146 

945 

1137 

1 9£ 

£0 . 3’. 

-0.47 

£0.37 

PC  71 

1 35 

34  0 

4 1 £ 

i iZ. 

£1  .£'. 

-4 .17 

CD  . 3'. 

DC  7 £ 

1 34 

34  0 

4£4 

34 

£4 

c.  • r • 

£1  .47 

DE7£ 

1 33 

1 73 

1 7£ 

0 

0 . ON 

-9 .57 

1 0 .57 

GR7£ 

94 

761 

q 2 3 

1 6£ 

£1  . 3'-. 

-4.77 

C l'  ml’ 

KS7£ 

1 3 0 

9 36 

1 04  0 

1 04 

11  . 1 ‘ . 

—y . 6*. 

£1  .67 

ME71 

167 

4 03 

48  1 

7?  Z; 

1 7 . 9'. 

5 . 07 

1£  . 37 

M I 7c! 

165 

575  3 

65  08 

165 

1 3 .in 

-0.57 

1 Z>  7 V 

1 • 1 B 

MO  7 £ 

1 3 0 

4 34 

5£  0 

y6 

1 9 . 3\ 

-11  .£7 

34  B qv 

NH7£ 

159 

1 6 0 

£79 

1 19 

74.47 

3.47 

6 y . 6’. 

MM7£ 

1 07 

3 8 6 

396 

1 0 

8 .6". 

— O 7V 

•_*  • 1 e 

1£  .47 

□K7£ 

1 1 1 

7 0£ 

859 

157 

££.47 

“6 . 7". 

31  .£7 

F'R7£ 

144 

3£  07 

344  0 

iZ.  Z*  Z' 

"7  Z* 

-3.1 7 

10.77 

SC  7 1 

1 09 

559 

63  7 

6*  8 

1£  .£7 

-14.97 

31  .37 

SC  72 

1 06 

593 

655 

6*8 

1 0 .57 

-5.97 

17.37 

TX7£ 

Ou 

£418 

d 9 8 5 

567 

£3.47 

-5.  £7 

3 0.37 

UT7 1 

1 3 3 

631 

735 

154 

£4  .47 

f ■ iZ.  • 

34 .17 

VT7 1 

164 

iz!ii!4 

£97 

"7‘Z* 

32  .6*. 

Z*  ■ iZ 

£3.57 

WR7£ 

1 7y 

£45  0 

26*87 

iZ  _•  r 

9.77 

0 .57 

9 . 1 7 

TOTRL 

3647 

iz' 775 

3 3£4£ 

4467 

15  . 5£7 

-4 . 097 

£0  .457 

TABLE  7-8.  EVENING  (4  P.M.  TO  10  P.M.)  DDST 
INJURY-CAUSING  ACCIDENTS 


STRTE 

Z. 

YERP 

MO. 

EQVL 

PRYS 

RlTURL 
RCC I PENT 
TOTRLS 

EQVL  PRY 
RCC I DENT 
TOTRLS 

CHRNGE 

RPPRPEMT 

PERCENT 

CHRNGE 

MIPPRY 

+MIGHT 

CHRNGE 

NET 

PERCENT 

CHRNGE 

CR7£ 

193 

338 1 1 

37165 

3354 

9 . 9 "a. 

4 . 5'. 

5 . £N 

CT7£ 

£07 

4997 

5£  1 5 

£ 1 3 

4 .4*: 

C • j 

£ . ON 

DC  7 1 

£ 09 

1 38  0 

1575 

195 

1 4 . r . 

-J  B ?v 

1 0 . ON 

DC7  £ 

£11 

1465 

1671 

£06 

14.1*. 

6 . 5'. 

7 . 1 N 

PE7£ 

£07 

3 1 1 

9£  0 

1 09 

13  .4". 

14  .8N 

- 1 . £N 

GR7£ 

178 

4177 

4713 

j 36 

1 £ . 8 ’. 

9 -5*. 

3 . ON 

K S7£ 

£07 

4 £ £ 7 

4337 

1 1 0 

8 .6’. 

1 . 3‘. 

1 . 3m\ 

ME  71 

£1  0 

1 4 06 

15  07 

1 01 

l'  ■ tZ  B 

- 0 . 4‘. 

7.6N 

M 1 7£ 

££  1 

£4908 

£549  0 

53  £ 

z*  z>  ■.-■ 

1 ■ _•  B 

0.  IN 

(Z  ■ iZ  'o 

MD7£ 

£04 

1434 

1776 

£9£ 

1 9 . 7’. 

15.4N 

•z*  "7  v 

— ' a 1 o 

NH7£ 

£17 

1 038 

1 1 08 

£0 

1 . 3‘. 

-£  . 5N 

4 . 5N 

NM7£ 

133 

££3£ 

o 3 9 8 

166 

7.4'. 

1 4 . 5N 

6*  • c.  a 

□K7£ 

19£ 

£997 

3156 

159 

5 . 3'. 

1 8 . 3". 

-11.  ON 

PR7£ 

£07 

13467 

14745 

1273 

9 . 5’i 

5 . 6N 

z- . 6"-. 

SC  7 1 

13£ 

£394 

£60  0 

£06 

8 . 6‘> 

1 £ . 4N 

-3.4N 

SC7£ 

134 

8616 

£6£9 

13 

0 . 5 

£ . 9N 

-£  . 4N 

TK7£ 

163 

14710 

15595 

335 

6 . O'. 

3 . ON 

£ . 9N 

UT7 1 

£06 

3054 

3 £68 

£14 

7 . 0"i 

d n \ '• 

4 . £N 

VT7 1 

865 

993 

1011 

18 

1 .8". 

-4.3N 

6 .4N 

WR7£ 

Z*  Z*  ~7 

9 3£8 

9££3 

- 1 0 0 

-1  .IN 

-4  .IN 

z*  . kz! 

TOTRL 

4 033 

131545 

140107 

356£ 

6 .51*: 

3.38N 

3 . 03N 

56 


TABLE  7-9.  NIGHT  AND  MIDDAY  DDST  IN JURY - CAUS ING  ACCIDENTS 
(DDST  MORNING  NORMALIZATION) 


STATE 

ACTUAL 

EQUIVALENT 

CHANGE 

% CHANGE 

CA  72 

30870 

28368 

-2502 

-8.1 

CT  72 

5020 

5002 

-18 

-0 . 4 

DC  71 

1234 

1184 

-50 

-4.1 

DC  72 

1176 

1208 

32 

2 . 7 

DE  72 

821 

743 

-78 

-9.5 

GA  72 

3103 

2956 

-147 

-4.7 

KS  7 2 

4002 

3657 

-345 

-8.6 

ME  71 

1898 

1992 

94 

5.0 

MI  7 2 

27669 

27519 

-150 

-0.5 

MO  72 

1743 

1548 

-195 

-11.2 

NH  72 

1045 

1081 

36 

3.4 

NM  72 

1507 

1376 

-131 

-8.7 

OK  72 

2620 

2444 

-176 

-6.7 

PA  72 

13262 

12846 

-416 

-31.5 

SC  71 

2023 

1722 

-301 

-14.9 

SC  72 

2063 

1942 

-121 

-5.9 

TX  72 

9674 

9167 

-507 

-5.2 

UT  71 

2410 

2236 

-174 

-7.2 

VT  71 

1172 

1209 

37 

3.2 

WA  72 

10344 

10400 

56 

0.5 

TOTAL 

123656 

118600 

- 5056 

-4.09% 

57 


72 

72 

71 

72 

72 

72 

72 

71 

72 

72 

72 

72 

72 

72 

71 

72 

72 

71 

71 

72 


7-10.  MIDDAY  AND  NIGHT  DDST  INJURY-CAUSING  ACCIDENTS 
(DDST  EVENING  NORMALIZATION) 


ACTUAL 

EQUIVALENT 

CHANGE 

% CHANGE 

46975 

49094 

2119 

4.5 

6787 

6945 

158 

2 . 3 

1815 

1883 

68 

3.7 

1855 

1975 

120 

6.5 

1083 

1243 

160 

14.8 

4807 

5266 

459 

9.5 

6041 

6121 

80 

1.3 

2261 

2253 

-8 

^1" 

o 

1 

36459 

36499 

40 

0.1 

2260 

2609 

349 

15.4 

1391 

1356 

-35 

-2.5 

2251 

2578 

327 

14.5 

3929 

4648 

719 

18.3 

17659 

18655 

996 

5.6 

2805 

3154 

349 

12.4 

3370 

3469 

99 

2.9 

18952 

19519 

567 

3.0 

3407 

3498 

91 

2.7 

1573 

1505 

-68 

-4.3 

13327 

12781 

-546 

-4.1 

179007 

185051 

6044 

3.381 
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TABLE  7-11.  IMPACT  OF  DDST  ON  YEARLY  ACCIDENTS 
(USING  APPARENT  PERCENT  CHANGE) 


STATE 

& 

YEAR 

YEAR  WITH 
S I X- MONTH 
D . D . S . T . 
(ACTUAL) 

YEAR  WITH 
FULL-  (EAR 
D . D . S . T . 
(.CALCULATED) 

CHANGE 

PERCENT 

CHANGE 

CATE 

1 3 0349 

191496 

1 0647 

5 . 89*. 

CT7£ 

£3511 

£44£6 

915 

3 b 39'-. 

DC?  1 

3745 

7 3 1 £ 

567 

3.41 

DC  7 c 

6811 

74E5 

614 

9 . 0 1 N 

DE7  £ 

3967 

4165 

193 

4 . 99‘-. 

bA7c 

£3189 

£5  06  0 

1371 

3 . 07  N 
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TABLE  7-12.  IMPACT  OF  DDST  ON  YEARLY  ACCIDENTS 
ruSING  NET  PERCENT  CHANGE) 
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D . D . S . T . 
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changes  are  shown. 

Table  7-8  shows  the  evening  results  in  a form  similar  to  the 
morning  results.  Tables  7-9  and  7-10  show  how  the  morning  normal- 
ization and  the  evening  normalization  were  obtained. 

The  overall  results  from  Tables  7-7  and  7-8  show  a Net  Per- 
cent Change  of  +20.45%  in  morning  in j ury- caus ing  accidents  for  the 
year  (based  upon  an  Apparent  Percent  Change  of  +15.52%),  and  a Net 
Percent  Change  of  +3.03%  in  evening  injury-causing  accidents  for 
the  year  (based  upon  an  Apparent  Percent  Change  of  +6.51%).  These 
results  indicate  a large  increase  in  accidents  under  DDST  in  the 
morning  and  a small  increase  in  the  evening. 

Tables  7-11  and  7-12  show  the  impact  of  DDST  on  the  yearly  in- 
jury-causing accident  totals  under  the  assumption  that  the  acci- 
dents in  the  morning  and  evening  periods  are  changed  by  the  per- 
centages shown  in  Tables  7-7  and  7-8  and  that  the  accidents  in  the 
other  parts  of  the  day  remain  unchanged  for  the  year.  For  each 
state-year,  the  actual  year's  accidents  are  shown  and  compared  with 
the  calculated  DDST  yearly  accidents.  Table  7-11  uses  the  Apparent 
Percent  Change  (from  Tables  7-7  and  7-8),  and  Table  7-12  uses  the 
Net  Percent  Change  (from  Tables  7-7  and  7-8). 

The  yearly  in j ury- caus ing  accident  totals  show  a normalized 
Net  Percent  Change  of  +4.01%  (based  upon  an  Apparent  Percent  Change 
of  +4.64%).  This  result  indicates  an  increase  in  yearly  accidents 
of  the  order  of  4%  if  DDST  is  introduced. 

It  should  be  noted  that  these  results  are,  as  for  YRDST,  only 
estimates  based  upon  data  from  a not  necessarily  representative 
group  of  states,  and  include  the  various  inaccuracies  of  the  Equi- 
valent Day  Methodology.  It  should  especially  be  noted  that  the 
morning  increase  is  based  upon  April  through  September  days,  while 
the  evening  increase  is  based  upon  September  through  April  days. 
Thus,  the  results  should  be  considered  only  as  an  indication  of 
a large  increase  in  April  through  September  morning  accidents  and 
a small  increase  in  September  through  April  evening  accidents, 
yielding  a net  change  in  yearly  injury-causing  accidents  of  about 
4%,  and  the  calculated  percent  changes  should  not  be  taken  as  ex- 
act numbers. 
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APPENDIX  A - THE  PROGRAMS 

A,  1 THE  EQUIVALENT  DAY 
PROGRAM/  EQDAY 
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C EQDAT 

c 

C F INUS  THE  EQUIVALENT  DAY  FOR  EACH  NON-D.S.T.  DAY 
C 

C THE  EQUIVALENT  DAY  IS  CHOSEN  TO  BE  ON  THE  SAME  DAY  OF  THE 

C WEEK  AS  THE  GIVEN  DAY. 

C HOLIDAYS  ANC  I HE  NIGHTS  BEFORE  HOLIDAYS  ARE  NOT  USED  AS  DAYS 
C UR  AS  EQUIVALENT  DAYS. 

C NO  EQUIVALENT  MOKMNG/EVENING  IS  USED  MORE  THAN  FOUR  TIMES. 

C 

C 

IMPLICIT  INTEGER  t A - /i  > 

UIMENSION  MTHBGN(12) . MTHEND  ( 12 ) 

DATA  FITHtSGN/1. 32. 61. 92. 122. 153. 183. 1:14,245.275. 

1 306,336/ 

DATA  MTHENL/31. 60. 91. 121. 152. 182. 213. 244. 274, 305. 335, 

1366/ 

UIMENSION  STATE (5) .HOL (50 ) .M0SED(366 ) ,EDSED<366> 

UIMENSION  RISEI366).  SETC366).  RUSAGE(366>.  SUSAGE(366) 
LOGICAL  LEAP 

DATA  MDSLD.EUSED.RUSAGE.SUSAGE/366*0.366» 0.366* 0.366*0/ 

READ  (6.601)  STATE. YEAR. NHOL 

601  FORMAT  (AAA, A2. IX. 14, IX. 12) 

WRITE  <6.691 ) STATE. YEAR .NHOL 

691  FORMA  I (1H1.AAA.A2.1X.IA.1X.I2) 

READ  ( 6 , b02 ) (HOL(I) .1=1. NHOL) 

602  FORMAT  (20(13. IX)) 

HR  1 TE ( 6 . 692 ) (HOL ( I ). 1 = 1, NHOL) 

692  FORMAT  ( 1 X , 2 0 ( 1 3 . 1 X ) ) 

UO  AU1  M=I,12 
M1=MTHBGN(M) 

M2  = M I HENU  ( M, ) 

REAL  (6,603)  (RISE(I),  I=M1,M2) 

603  FORMAT  (16(14, IX)) 

WK 1 I L ( 6 « 693 ) (RISE(I),  I=M1,M2) 

693  FORMAT  (16(1X,I4) ) 

AO  1 CONTINUE 

DO  AU2  M=1 , 12 
M1=M I HEGN ( M ) 

M2=M(HLNU(M) 

READ  (6.603)  ( SET ( I ) . I=M1,M2) 

WR1TL(6,693)  ( SET ( I ) , I =M1 , M2 ) 

A02  CONTINUE 
C CHECK  TEAR 

IF  ( TEAR.EG.1972IG0  TO  1972 
IF  (YEAR. EC. 1971) GO  TO  1971 
1972  LEAP=,TRUE, 

USTBGN=121 
USTENU=302 
GO  TU  701 

1971  LEAP=. FALSE. 

DST bGN=116 
DST  ENU  = 3UA 
GU  TU  701 

C ADD  ONE  HOUR  FOR  DST  DAYS 
701  DO  A99  I=DSTBGN,CSTEND 
R1SE11)=RISE(I)+100 
SET  II  )=SET ( 1 1+100 
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499  continue: 

IF  (LEAP)  GO  TO  703 
C NON-LEAP  YEAR 
C DELETE  F Lb  29 
RISE(b0)=0 
SET ( bU )=0 
GO  TO  703 
C LOAF  TLAK 

C REMOVE  HOLIUATS  AND  NIGHT  BEFORE  HOLIDATS 
703  UO  403  I=l.NHOL 
HOLDAT  = HOL( I ) 

RISE  <H0L0AT)=0 
SET (HOLDAY  >=0 
IF  (HOLDAT. NE.l)  GOTO  702 
SET (366)=0 
GOTO  403 

702  SET (HULUAT-1 >=0 
403  CONTINUE 
C FIND  EUUiVALENT  DAYS 
DO  409  1=1,366 

IF  (I.GE.DSTBGN.AND.I.LE.DSTEND)  GO  TO  409 


MDSEU ( 1 ) =EGU1V (RISE. I .RUSAGE .LEAP ) 
EOSEU ( 1 )=EGU1V (SET, I . SUSAGE . LEAP ) 

409 

CONTINUE 

C PRINT  OUIPUT 

WRITE! 9.617 (TEAR. ST  ATE 

617 

F ORMAT ( IX , 14 , * EQUIVALENT  UAYS  FOR 
WR1TE(6.607)YEAR. STATE 

• .4A4.A2//) 

607 

FOkMAI  (1H1.I4.21H  EUUIVALENT  DAYS 
WKITEI9.604) 

WRITE(6,604) 

FOR  .4A4.A2//I 

604 

FORMAT  (2X.* DAY*. 7X«* MORNING' «2X.' 
00  411  1=1,366 

EVENINGV13X,  ,EQDAY’.4X.'EQDAY’ 

IF  (l.GE.USTbGN.AND.I.LE.OSTENO)  GOTO  411 
IF  ( 1.EQ.60.AND. .NOT .LEAP )G0T0  411 
M=MDSEU ( I ) 

E=EDSED ( 1 ) 

WRlTE(9.bl4 ) 1 • H , E 
WH1TE(6. 614)1. M.E 

614  FORMAT (2X,13»*)',8X,13,6X,I3) 

411  CONTINUE 

WRITE<6,607> 

WR1TE(6,604) 

DO  410  1=1.366 

IF  (I.GE.DSTBGN.AND.I.LE.DSTEND)  GO  TO  410 
IF  ( 1 .EQ.60 .AND. .NOT .LEAP )GO  TO  410 
J=MDSEU ( I ) 

K=EDSED ( 1 ) 

M1  = KI'JT  HNU  ( I ) 

ni=datno(I) 

M2=MNT  HNU ( J ) 

N2  = DAT  NO ( J ) 

M3=MNT  hNU ( K ) 

N3  = DA  Y NO ( K ) 

WRITE!6.605)  M1.N1.M2.N2.M3.N3 
605  FORMAT  (1H  .I2,1H/.I2,7X.I2.1H/,12«4X«I2.1H/.I2) 
410  CONTINUE 
END 
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C FINU  MORNING/EVENING  EQUIVALENT  DAT 

INTEGER  FUNCTION  EQUI V ( R ISET . I . USAGE . LEAP ) 
IMPLICIT  INTLGEK(A-Z) 

LOGICAL  LEAP 

DIMENSION  RISETI366) .USAGE(366) 

IF  (K1SLT (I ) .EQ.O)  GOTO  706 
UYOF  WK=  I - ( 1 / 7 ) * 7 

IF  (.NOT.  LEAP  .AND.  I . GE . 60 ) OYOF WK=UTOF WK-1 
IF  (UYOFWK.LL.O  ) DY0FWK=DY0FWK+7 

dstime=riset(i>+iou 

EUU I V = U 
MINU1F=5 

UO  All  J=DY0FWK.366.7 
EQ  = J 

IF  (.NOT. LEAP  .AND.  J . GE . 60  > EQ  = J+ 1 
IF  (EU.G1 ,366)GGT0A11 
IF  (USAGL(EQ).GI .3)  GOTO  All 
0IFF=GETU1F (UST IME.RISET (£Q ) ) 

IF  (UlFF.GT.MNUIF)  GO  TO  All 
IF  (UIFF.LT.MINOIF)  GO  TO  705 
IF  ( APART ( 1 ,EU) .GT. APART ( I .EQUIV ) ) GO  TO  All 

705  EGUI V = EQ 
MIND1F=D1FF 

All  CON11NUE 

IF  (LUUIV.GT.O)  USAGE(EQUIV)=USAGE( EQUIV)  + 1 

RETURN 

706  EGUIV=0 
RETURN 
ENU 


C GET  TIME  DIFFERENCE  IN  MINUTES  BETWEEN  TWO  GIVEN  TIMES 
INTEGER  FUNCTION  GET01F  ( OST IME »UT IME ) 

INTEGER  1 1KE1 .TIMES. HOUR 1.H0UR2.0STI ME 
IF (USTIME. LT.JTIMEIGO  TO  712 
T IME1=UT IME 
T IME2=UST IME 
GO  TO  711 

712  TIMEi=USTIME 
T I M E 2 = u T IME 
711  H 0 U R 1 = T 1M El/100 

H0UK2= I IML2/100 
GET01F=(H0UR2-H0UR1 )»60 
MI  NUT  1-T  1 ML  1 -HOUR  1*  1(J0 

MINU1 2=T lME2-HGUR2*lu0 

GEIUiF=6ETC,IF  + MINUT2-MINUTl 
IFIGETU1F.GT.720)  GETDIF=1AAO-GETOIF 

return 

ENU 


C FIND  HOW  FAR  APART  TWO  DAYS  ARE  (BASED  ON  WEATHER) 
INTEGER  FUNCTION  APART  (0AY1.DAY2) 

INTEGER  DAY1.DAY2 
DATA  P1/3.1A159/ 

UY1=UAY1 

UY2=UAY2 

WE T HR1=C0S ( 2. 0*PI*(DTl-32. 01/366.0) 

— WETHR2=COS(2.0*PI*(DY2-32.0)/366.0) 

APAR1 =ABS(NETHR2-WETHR1 ) «10U0 .0 
IF  (UAT2.LO.U)APART=5000 
RE1 URN 
ENU 
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C FIND  MUN1H  NDMEER  FROM  DAYCOOE 
FUNCTION  MNTHNOtCAYCOL ) 

INTEGER  DAY  CDF 

iftuaycdl.e:q.o)rnthno=o 

IF (DAYCDL.GT.U.AND.DAYCOE.LE.3J )MNTFiNO  = 1 
IF ( DA YCDL.GT. 31. AND. OA YCOE.LE. 60 )RNTHN0=2 
IF (UAYCUE. GT.60. AND. QAYCUL.LE. 91 >RNTFiNO  = 3 
_ IF(UAYC0E.GT.91.AND.0AYC0E.LE.121)MNTHN0  = ‘F 
IF (UAYCUt.GT.121. AND.DAYCUE.LE.lb2 )MNTHN0=5 

IF  (UAYCUE. GT. 152. AND. UAYCUE. LE.182)FNTFiN0  = 6 

1FIUAYCDL.GT.1B2.AN0.UAYCDL.LE.  213)  F,  NT  UNO  = 7 
IF  (UAYCUL.GT.213.AND.DAYCDE.LE.2A‘4)FiNTFTN0  = 8 
IFlUATCUL.GT. 244. AND. UAYCDL. LE.274)FNTFiN  0 = 9 
IF(DAYCDL.GT.274.AND.DAYCDE.LE.305)MNTF1NO=10 
IF  (UAYCUL.GT.305.  AND  . UAYCDL  . LE  . 335  )KNTHNO=ll 
IF(UAYCUL.GT.535)MNTFIN0  = 12 
RETURN 
LNU 


C FIND  OAT  NUMBER  FROM  DAYCOOE 

INTEGER  FUNCTION  DAYNO ( DAYCUE ) 

INTEGER  UAYCUE 

IF (UAYCUE. LE.31)CAYN0=DAYCDE 

IF (UAYCDL.  GT.31.AND.UA YCOE.LE.60  S OA YNO=DA YCOE-31 
IF ( DAYCUE • GT .60.AND.OAYCDE.LE.91)OAYNQ=DAYCDE-60 
IF (UAYCUE.GT.91. AND . D A YCDE . LE . 12 1 ) OA YNO=0 AYCDE-91 
IF (OAT  COL.  GT.121.AND.0A YCOE.LE. 152) U A Y NO=UA YCUE-121 
IF (UAYCUE. GT .152. AND. DAYCUE. LE.182)UAYM0=UAYCUE-152 
IF (UAYCUE.GT.ia2.AND.DAYCUE.LE.213)UAYN0=0AYCUE-182 
IF (OAYCOL.GT.213.AND.OAYCOE.LE.24‘4  )0AYN0  = UAYCDE-213 
IF ( UAYCUE. GT .244. AND . 0 A YCOE . LE . 2 74 ) UAYN0=UAYCDE-244 
IF  (OATCUL.GT  .274  . ANO.  DA  YCOE.LE.  305  )U  AT  N0=0A YCUE-274 
IF ( UAYCUL.GT. 305. ANO. U A YCOE.LE. 335)0. A rsO-UAYCUE-305 
IF ( UAYCUE. GT. 335) DA YN0=UAYCUE-335 
RETURN 
LNU 
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A. 2 THE  QUARTER  DAY  PROGRAM,  QDY 
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C QUARTER-DAY  STATISTICS 
C 

C FINCS  THE  NUMBER  OF  ACCIDENTS  FOR  EACH  QUARTER-DAY 
C OF  THE  YEAR 

• i 

C 

C 

C (INCLUDES  FEB.  29  ENTRY  IN  ARRAYS) 

C 

C 

IMPLICIT  INTEGER  (A -7) 

DIMENSION  DYINMH ( 12 ) .DYBFMHI12) 

DIMENSION  NEVE  AC ( 366 . 6 ) 

DIMENSION  T0T46 ( 366 ) 

DIMENSION  TOT610(366> 

INTEGERS  NUMBER  ( 256  ) 

INTEGER*2  ACSVIX<256) 

INTEGER*?  KI(2)/2*0/.NNF(2)/2*0/.M10(2)/2*0/.M1(2)/2*0/. 
1D10(2)/2*0/,D1(2>/2*I)/,TMHR10(2)/2*0/.TMHR1(21/2*0/. 

2TMMN10 (2  )/2*0/,TKMNl (2 )/2*0/ 

INTEGER*2  L(12) 

data  blanks  /zocoooono/ 

DATA  NEVEAC  , T0T46 , TO T61 0/2 196*0 . 366*0 . 366*0/ 

DATA  DYBFMH/O  .31 ,60 ,91 ,121 ,152, 182, 2 13,244, 279, 305. 335/ 

DATA  DYINMH /31 , 29, 31. 30, 31, 30, 31, 31, 30, 31, 30. 31/ 

EQUIVALENCE  (KI(1).BYTS4A),(NNF(1),BYTS4B).(M10(1).RYTS4C) 
EGUI VALENCF  (Ml ( 1 ) , 0 YTS40 ) .(010(1). BYTS4E ) , (Dl(l) ,BYTS4F! 
EQUIVALENCE  ( TMHR10 ( 1 ),RYTS4G).(TMHR1(1).BYTS4H) 

EQUIVALENCE  (TKMN10 ( 1 ),BYTS4I),(TMMN1(1),BYTS4J) 

DATA  ACSVIX/192*-20,10,ll,12,13,14,15,16,17,8*-20, 

118, 19, 20. 21  . -20. -20. 22, -20, 23. 8*-20. 24, 25, 

2-20. 26, 27, 28. 29, 7*-20, 0.1, 2, 3, 4, 5, 6, 7. 8,9, 

37*-20/ 

DATA  NUMBER/239*-20,0,l,2,3.4,5.6,7.8,9,7*-20/ 

DIMENSION  NFATAL ( 366 .4 ) .NCASUL  ( 366 ,4  ) 

DATA  NFATAL. NCASUL/1464*0. 1464*0/ 

DATA  NZEKO.NERROR.NFACC.NINUAC.NFTLTY.NINJRO/O, 0,0. 0.0,0/ 

C 

CALL  ERRSET (218 ,250 ,250 , 1 ) 

C 

C LOOP  TO  READ  AND  FILL  ARRAYS 
C 

1111  RE AD (9, 16  Cl ,END  = 150  0 ) K I (2) ,NNF(2) ,M10(2) , Ml ( 2 ) , D 1 0 < 2 > , 01 < 2 ) , 
1 TMHR10 (2)  ,TMHR1 (2)  »TMMN10(2)  ,TMMN1(2) 

1601  FORMAT (2X,A1,4X,A1,10X,3X,4A1,3X.4A1.4X) 

INJK=ACSVIX ( (BYTS4A-RLANKS 1/256) 

NNFINJ=ACSVIX( (BYTS4B-BLANKS 1/256) 

MNTH10  = NUMPER( ( BYTS4F-BL ANKS  ) /256 ) 

MNTH1=NUMBER ( (BYT S4D -BLANK S)/2 56 ) 

- - DAY10=NUM3ER ( ( BYTS4E-BL ANKS ) /256 ) 

DAY1=NUMUER((BYTS4F-RLANKS>/256) 

HR10=NUMdER( ( BYTS4G-RL ANKS 1/256) 

HR1=NUMBER ( (BYTS4H-BL ANKS ) /256 ) 

MIN10=NUM8f  R ( ( H YT  S4  I -BL  ANK  S ) /256  ) 

MIN1=NUM8ER ( ( BYTS4 J-RL ANKS ) /256 ) 

M0NTH=10*MNTH10  + MNTH1 
CYCFMH  = 10*r.AY10  + D A Y 1 
HCUR  = 1 00  0 *).R  1 0 + 1 0 0 *HR  1 
MIN=10*M1N10  ♦ MINI 
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time=hour  + min 

IF  (INJK.LT.O  .OR.  NNFINJ.LT. 0 .OR.  MONTH. LT . 1 
1 .OR.  MONTH. GT.12)  GO  TO  1405 

IF  (HOUR. GT. 2300  . OR  . HOUR  . L T . 0 > GO  TO  l<f05 
IF  (DYOFMH.LT. 1 .OR.  DYOFMH.GT .DYINMH(MONTH) ) GOTO  1405 
IF  ( MIN.GT. 59  .OR.  MTN.LT.O)  GO  TO  1405 
DAY=DYOFMH  + DYBFMH ( MONTH ) 

IF  (TIME.GF.400.AND.TIMF.LT.1000)QDAY=1 
IF  ( TIME.  GF.  1000  . AND.T IMF.LT. 1600 )QDAY  = 2 
IF  ( TIME.  GF.160  0.AND.TIME.LT.  2200  )(5DAY  = 3 
IF  (TIME. GF. 2200. AND. TIME. LE.2359)0DAY=4 
IF (TIME.GE.O. AND.TIMF.LT. 4 00) GOTO  1701 
GOTO  1702 

1701  ODAY=4 
DAY=DAY-1 

IF  (DAY. EG. 0 )DAY  = 366 

1702  IF (INJK.EQ.O)GOTC  1703 

NFATAL (DAY.QDAY ) =NF AT AL ( D A Y , QD A Y ) +1 

NFACC=NFACC+1 

NFTLTY  = NFTL  TY  + INJK 

GOTO  1704 

1703  IF(NNFINJ.EQ.O)  GOTO  1798 

1704  NCASUL(DAY.QDAY)=NCASUL(DAY,QDAY)+1 
NINJAC=NINJAC+1 
NINJRD=NINJPD+NNFINJ+INJK 

IF (GDAY.NE.3)G0T0  1111 
EVEHR=10*HR10  + HR1  - 15 
NEVE AC (DAY .EVEHR )=NEVEAC ( D A Y . E VEHR > +1 
GOTO  1111 

1405  NERR0R=NERR0R+1 
GOTO  1111 

1798  NZERO=NZERO+l 
GOTO  1111 
C 

C WRITE  OUT  ARRAYS 
C 

1500  WRITE (6.1608) 

1609  FORMAT ( 1H1 ) 

WRITE(6.1603) 

WRITE(8.1603) 

1603  FORMAT!////’  NO.  OF  ACCIDENTS  WITH  FATALITIES  (K)’///) 

WRITE (6,1609) 

_ WRITE(8,1609) 

1609  FORMAT ( IX . ’DAY’ , 9 X . • 4 AM- • . 6X . ’ 1 0 AM  • , 6X . • 4PM- • ,6X . • 1 0PM • / 

1 IX,  ’NO. ’ ,9X, ’10AM’ ,6X, ’-4PM* ,6X. *10PM’ ,6X, ’-4AM’/) 

WRITE  (6.1602HI . < NFATAL  ( I . J > . J = 1 ♦ 4 ) , 1 = 1 . 366 ) 

WRITE (8, 1602) (I. ( NFATAL ( I . J ) ,J  = 1.4) ,1  = 1.366) 

1602  FORMAT (( IX . 13. • : ',4(1X,I9))> 

WRITE  <6.161U)NFACC.NFTLTY 
WRITE  (8,1610)NFACC,NFTLTY 

1610  FORMAT!//’  NUMBER  OF  FATAL  ACCIDENTS  =’.17/ 

1 ’ NUMBER  OF  FATALITIES  =’.I7) 

WRITE(6,16U8) 

WRITE(6.1604) 

WRITE(8,1604) 

1604  FORMAT!////’  NO.  OF  ACCIDENTS  WITH  ANY  INJURY  (K.A.B.C) 

WRITE  (6,1609) 

WR ITE ( 8 , 1609 ) 


•///) 


69 


WR  ITU  6,1 602)  ( I . (NCASUL  ( I . J)  . J=1.4)  . 1 = 1.366) 

WRITE(8,1602) (I. (NCASUL ( I . J ) . J=1 . 4 ) , 1=1.366) 

WRITE(6«1611)NINJAC,NINJR0 

WRITE  (8,1611  ) N IN J AC  « N IN  JRD 

1611  FORMAT!//*  NUMBER  OF  INJURY-CAUSING  ACCIDENTS  =*,17/ 

1 • NUMBER  OF  INJURED  (INCL.  FATALITIES)  =*.I7) 

WRITE(6,1607)N2ERO.NFRROR 

WRITE(8.1607)N2ERO.NERROP 

1607  FORMAT ( ////23H  NC-INJURY  ACCIDENTS  = .I7//10H  ERRORS  = .17) 

DO  14  01  1 = 1,366 

T 0 T 4 6 ( I ) = N F V E A C ( I ,1 ) + N E V E A C ( I .2) 

TOT610(I)=NEVEAC(I.3)+NEVEAC(I,4)4NEVEAC(I.5)+NEVEAC(I.6) 

1401  CONTINUE 

WR I TE ( 6 , 1608 ) 

WRITE (6.1612) (I . (NEVFAC(I.J) , J=1.6) .T0T4  6 ( I ) « TOT610 ( I ) .1  = 1.366 ) 
WRITE (8. 161 2) ( I . (NEVFAC( I . J) . J=1 .6) .T0T46(I ) , T0T61 0 ( I ) .1=1.366 ) 

1612  FORMAT (//12X, * HOURLY  EVENING  INJURY-CAUSING  ACCIDENT  TOTALS*/// 

1 IX,  *DAY* , 1 X » 4 X » • 4-5 • ,4X. *5-6*  «4X, *6-7* .4X. *7-8* , 4x . * 8-9 • , 3X . 

2 *9-10* .6X.3X. *4-6* .3X. *6-10*// 

3 ( IX, I3,6( IX, 16) ,5X .2< IX.  16)  > ) 

END 
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A. 3 THE  COMPARISON  PROGRAM/  COMPAR 


C COMPAR 
C 

C ROUTINE  TO  COMPARE  THE  ACCIDENT  TOTALS  OF 
C NON-D.S.T.  DAYS  WITH  THOSE  OF  THEIR  EQUIVALENT  DAYS 
C 
C 

C ACTUAL  EQUIVALENT  DAY  ACCIDENT  TOTALS  ARE  USED. 

C 

C 

IMPLICIT  INTEGER  ( A-7  ) 

REAL  MPCTCG,EPCTCG,YPCTCG 
REAL  PCEVAC  .NCGEVA 

REAL  MNMPCG ,ENMPCG,MNETCG,ENETCG,YNTPCG 

REAL  RMPCTC  ,REPCTC 

REAL  RMNMPC .RENMPC 

REAL  RMNETC .RENETC ,R YPCTC .RYNTPC 

REAL  QAVDWK , AVDWK , DMNRM , OENRM , MNRM . ENRM , Q AVMNH, AVMNH 
LOGICAL  LEAP 

DIMENSION  NFATALI 366.4) . NC ASUL ( 366 . 4 ) 

DIMENSION  CLEVAC(C) . PCEV A C ( 8 ) .NCGEVA(6) 

DIMENSION  MNHALL ( 12  ) 

DIMENSION  EVAC<8) .EQFVAC(fi) tSEVAC ( 8 ) , SEQE  VA ( 8 ) 

DIMENSION  NEVE  AC < 366.6 ) ,T0T46( 366) .T0T61 0(366) 

DIMENSION  DST BGN ( 2 ) . DSTEND ( 2 ) 

DIMENSION  ST  AT  E ( 5 ) 

DIMENSION  CWKT0T(7,4>,DWKN0<7)»QAVDWK(7.4) . AVDWK ( 7 ) 
DIMENSION  MiNHTOK  12.4)  ,MNHN0(12)  . Q AVMNH  ( 12 . 4 ) .AVMNH  (12) 
DIMENSION  DYEFM2 ( 13 ) 

DIMENSION  E QDA  Y (366.3) 

DATA  LEAP/. FALSE./ 

DATA  DSTUGN/116,121/ 

DATA  DSTEND/304 ,302/ 

OAT A SMNON.SENON.SMDST .SEDST .SNOON.SMIDN/O, 0.0, 0,0.0/ 

DATA  DWK TOT .DWKNO .MNHTOT ,MNHN0/28*0 .7*0 .48*0,12*0/ 

DATA  MNHALL/12*0/ 

DATA  SEVAC,SEQEVA/8*0,8*0/ 

DATA  DYEFM2/0, 0.3 1,60, 91, 121, 152, 182, 213, 244, 274, 305, 335/ 
DATA  NERROR/O/ 

DATA  SMOYNM,SEDYNM,SMEQNM,SEEQNM/0. 0,0,0/ 

DATA  NMEQDS ,NEEODS/0 . 0/ 

C 

c 

c 

C GET  STATE  AND  YEAR 
C 

RE  ADI  10.2630 ) YEAR .STATE 
2630  FORMAT (1X,I4,21X,4A4.A2////) 

IF  (YEAR. EQ. 1972)  LE  AP= .TRUE  . 

C 

C READ  QUARTER-DAY  ACCIDENT  TOTALS 
C 

READ (8 ,2601  ) ( ( NFATAL ( I . J) ,U  = 1 .4 ) . 1 = 1 ,366  ) 

2601  FORMAT ( /////////// (7X, 4 (IX, 19) ) ) 

READ (8,2641 ) NF ACC , NFTL T Y 

2641  FORMAT! //2RX , I7/28X , 17  ) 

RE  AD l 8 , 2b0 1 ) ( ( NCASUL (I,J),J  = 1.4). 1 = 1*366) 
PEA0(8,2642)NINJAC»NINJRD 

264 2 FORMAT (//39X,  I7/39X.  17  ) 
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READ  (8,2652)  ( (NEVE AC  ( I »J)  » J=1 ,6 ) »T0T46  (I  > ,TOT610  (X  ) ,1  = 1 ,366  ) 
2652  FORM AT (////////////// (4X, 6 (IX, 16 ) ,5X,2( IX, 16) ) ) 

IF  (LEAP)  GOTO  2720 

NFATAL(59,4)=NFATAL(59,4)+NFATAL(60.4) 
NCASUL(59.4)=NCASUL(59.4>+NCASUL<60,4)  ' 

2720  WRITE(6. 2631) YEAR, STATE 
2631  FORMAT ( 'l* ,10X,I4,1X,4A4,A2) 

WRITE (6,2610)  (I, (NFATAL(I.J) ,J=1.4)*I=1  ,366) 

2610  FORMAT! • NUMBER  OF  ACCIDENTS  WITH  F AT AL IT IES * /// ( 1 X , 13 , 

1 3X,4(1X,I9) ) ) 

WRITE! 6. 2643) NF ACC. NFTLTY 

2643  FORMAT!//*  NUMBER  OF  FATAL  ACCIDENTS  =*,17/ 

1 ’ NUMBER  OF  FATALITIES  =*,I7) 

WRITE (6,2631 ) YEAR, ST ATE 

WRITE  (6.2611) (I, (nCASUL(I,J),J=1,4).I=1.366) 

2611  FORMAT!*  NUMBER  OF  ACCIDENTS  WITH  ANY  INJURY*/// 

1 (1x«I3,3x,4(1x»t9  ) ) ) 

WRITE (6, 2644 ) NINJAC * NINJRD 

2644  FORMAT!//*  NUMBER  OF  INJURY-CAUSING  ACCIDENTS  =’.17/ 

1 • NUMBER  OF  INJURED  (INCL.  FATALITIES)  =*.I7) 

C 

C GET  DAY-OF-WEEK  ANO  MONTHLY  TOTALS  AND  AVERAGES 
C 

DO  2420  1=1,366 

IF  (.NOT. LEAP  .AND.  I.EO.60)GOTO  2420 
DYOFwK  = WEEKDY (YEAR,  I ) 

MHOFYR=MNTHNO ( I ) 

DC  2421  J= 1 , 4 

DWKTOT (OYOFWK ,J)=DWKT0T (DYOFWK.J)  + NCASUL(I.J) 

MNHTOTIMHOFYR, J)=MNHTOT(MHOFYR. J)+NCASUL( I. J) 
MNHALL(MHOFYR)=MNHALl  < MHOF YR ) +NC ASUL ( I . J ) 

2421  CONTINUE 

DWKNO( OYOFWK ) =DWKN0 ( DYOFWK > +1 
MNHNO (MHOFYR )=MNHNO(MHOFYR)+l 
2420  CONTINUE 

WRITE! 6 ,2631  ) YEAR  .STATE 

WRITE  (6,262  0 ) (D,  (DWKTOT  (D.  J)  , J = l,4:  , ( D > ,D  = 1.7) 

2620  FORMAT!*  ACCIDENT  TOTALS  FOR  EACH  DAY  OF  THE  WEEK*// 

1 (3X.I2. * : * .4 (IX. 19) »3X. 'N0,DAYS=* , 19) ) 

WRITE (6,2631 )YFAR .STATE  . 

WRITE (6,2621 ) (M, ( MNHTOT ( M , J ) , J=1 , 4 ) , MNH ALL ( M ) , MNHNO ( M ) , M=1 , 12 ) 

2621  FORMAT! * ACCIDENT  TOTALS  FOR  EACH  MONTH*// 

1 (3X.I2, • : ’ ,4( IX. 19) ,3X. *TOT=* . 19 , 3X . • NO . DAYS= * . 1 9 ) ) 

DO  2426  D=1.7 
DO  2428  J= 1 , 4 

QAVDWK(D, J )=<!.* DWKTOT (D, J) )/DWKNO(D> 

2428  CONTINUE 

SIGDWK=DWKT0T(0.1 ) +DWKTOT ( D . 2 ) ♦DWKTOT < 0 . 3 ) +OWKTOT ( D . 4 ) 
AVDWK(D)=(1.*SIGDWK)/0WKN0(D) 

2426  CONTINUE 

DO  2427  M = 1 , 12 
DO  2429  J= 1 , 4 

GAVMWH ( M , J ) = ( 1 .*MNHTOT ( M . J ) ) /MNHNO (M) 

2429  CONTINUE 

SIGMNh  = MNHTOT ( M , 1 ) +MNHTOT ( M . 2 ) +MNHTOT ( M , 3 ) +MNHTOT ( M , 4 ) 

AVMNH ( M ) = ( 1 . *S IGMNH ) /MNHNO ( M ) 

2427  CONTINUE 

WRITE (6, 2631) YEAR, ST ATE 
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5 ’DAY  DATE.  ACCDNTS*  «9X»  ’MIDDAY  + NIGHTV//) 

SMACC=0 
SDMACC=0 

SEACC=0  ✓ 

SDE ACC=0 
SDYTOT=0 
SEOTOT  = 0 

DO  2404  I = 1 « NDA  YS 
DA Y=EODAY ( 1,1) 

ME0DAY=ECDAY(I,2) 

EEGDA Y=EUDAY ( I ,3) 

MHDAY=MNTHNO(DAY  ) 

DYDAY=DAY-DYBFM2(MHDAY+1) 

DWKDAYrWEEKDY (YEAR, DAY) 

MHMEG=MNTHN0< MEQDAY ) 

DYMEQ=MEQDAY-DYBFM2<  MHMEQ+1) 

DWKMEQ=WEEKDY ( YE AR . MFQDAY  ) 

MHEEQ=MNTHN0<EEQCAY) 

DYEED=EEQDA Y-DYBFM2 ( MHEEQ+1 ) 
DWKEEQ=WEEkDY(YEAR.EFGDAY) 

IF  (DAY.LT. 1 .OR.  DAY. GT. 366 

1 .OR.  MEQDAY. LT.O  .OR.  ME0DAY.GT.366 

2 .OR.  EEQDAY.LT. 0 .OR.  EEOD A Y . GT . 366 ) GOTO  2709 
IF  ( • NOT . L t AP  .AND.  (DAY.EGl.60  .OR.  MEODAY.E0.60 

1 .OR.  EEQDAY.EQ.60 ) ) GOTO  2709 

DYBEFR=DAY-1 
IF  ( DAY .Eq. 1 )0YBEFR=366 

IF  (DAY. tO. 61  .AND.  .NOT. LEAP)  DYBEFR=59 

MBEF0R=MEQDAY-1 

IF  (MEGDAY.EQ.l )MBEFoR=366 

IF  (MEQDAY.EG.61  .AND.  .NOT. LEAP)  MBEF0R=59 
C 

DMACC=NCASUL(0AY,1) 

DEACC=NCASUL (DAY. 3) 

C 

IF  (MEQDAY.EQ.O)GOTO  2705 
C MORNING  EQUIVALENT  DAY 
NMEQDS=NMEUDS+1 
MACC=NCASUL ( MEQDA Y , 1 ) 

SMACC=SMACC+MACC 

SDMACC=SDMACC+DMACC 

MEGNMD  = NCASUL(MBEF0R.4)-»NCASUL(MEQDAY,2) 
MDYNMD=NCASUL (DYBEFR .4 )+NCASUL (DAY ,2) 

GOTO  2706 

C NO  MORNING  EQUIVALENT  DAY 

2705  MACC=0 
MOYNMD=0 
MEQNMD=0 
DMACC=-DMACC 

C 

2706  IF  (EEQDAY.EQ.O)GOTO  2707 

C EQUIVALENT  EVENING 

NEEQDS=NEEGDS+1 
EACC  = NCASUL (EEqOAY  *3 ) 

SEACC=SEACC+EACC 

SDEACC=SDEACC+DEACC 

ECYNMD  = NCASUL (DAY  ,2  )+NCASUL  (DAY  ,4  ) 

EEQNMD  = NCASUL (EEQDAY .2  )+NCASUL (EEQOAY ,4  ) 
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GOTO  2708 

C NO  EQUIVALENT  EVENING 

2707  EACC=0 
ECYNMD=0 

EEQNMC=0  ' 

DEACC=-DEACC 
C TOTALS 

2708  SMDYNM=SMDYNM+MOYNMD 
SMEQNM=SMEQNM+MEQNMD 
SEDYNM=SEDYNM+EDYNMD 
SEEQNM=SEEGNM+EEQNMD 

WRITE ( 6.2602 ) DWKC A Y . WHO A Y . DYD A Y . DM ACC . DE ACC . MDYNMO . EDYNMD » 

1 DWKMEO  .MHMEQ  • DYMEQ.MACC. 

2 DWKEEG.MHEEQ.  DYEEQ . E ACC , MEQNMO . EEGNMD 

2602  FORMAT ( 1 X . 1 X . I1.1X.IP,'/',I2,1h:,I7,3X.I7,8X.'<'.I6.I6.2X, 

1 •>  • . 12X  « 1H (.IltlX.I2,«/*,I2.2H:).I7.3X.lH(.Il»lX»I2.,/,« 

2 I2,':)',I7,6X,'<',I6,I6,2X,'>') 

GOTO  2404 

2709  NERROR=NERROR+l 
2404  CONTINUE 

WRITEI6.2604) 

26  04  FORMAT  (///IX.  • . • , 

1 . ., 

2 . . ) 

WR ITE ( 6. 2603 )SDM ACC. SDE ACC. SMOYNM.SEDYNM.SMACC.se ACC. 

1 SMEQNM « SEE  QNM 

2603  F0RMAT(//1X,9X,I7.3X.I7.8X. •<* .16, I6.2X. •>• .12X.10X, T7.3X. 

1 10X, 17. 8X. *<• .16.16. 2X. »>• ) 

WRITE (6. 2646 ) NMEQDS . NEEQDS 

2646  FORMATI///////////////1X.40X, 'NO.  OF  MORNING  EQUIVALENT  DAYS 

1 I7//1X.40X. ‘NO.  OF  EVENING  EQUIVALENT  DAYS  ='.I7) 

MCHNG=SMACC-SDMACC 
ECHN6=SEACC-SDFACC 
MNMCHG=SME&NM-SMDYNM 
ENMCHG=SEEQNM-SEDYNM 

WRITE (6,2650 ) MCHNG . ECHNG . MNMCHG . ENMCHG 
2650  FORMAT(////1X.40X, 'CHANGE  IN  MORNING  ACCIDENTS  =».I7/ 

1 IX, 40X  . 'CHANGE  IN  EVENING  ACCIDENTS  =*.I7// 

2 IX,  40X,  'CHANGE  IN  N I GHT-M I OD A Y ACCIDENTS  =*.17/ 

3 IX, 40X, 'CHANGE  IN  MIDDAY-NIGHT  ACCIDENTS  ='.I7) 

WRITE (6,2638) YEAR, ST  ATE 

2638  FORMAT ( 'l'/////40X, • CHANGE  IN  THE  LEVFL  OF  ACCIDENTS'/ 

1 4 0 X . ' UNDER  THE  EFFECT  OF'/ 

2 40X , ' DAYLIGHT-SAVING  ALL  YEAR'/ 

3 40X . ’ FOR  '.14,*  FOR  '.4A4.A2///) 

MPCTCG=(MCHNG*100. l/SDMACC 

EPCTCG= (ECHNG* 100. )/RDE ACC 
RMPCTC=MPCTCG+0.00 
REPCTC=EPCTCG+0.00 
WRITE (6, 2637 JRMPCTC.REPCTC 

2637  FORMAT ( 3X , 40X . • MORNING  ACCIDENTS:  ',F6.1.'X'// 

1 3X.40X.  ' EVENING  ACCIDENTS:  '.F6.1.*%') 

MNMPCG=(MNMCHG*100. ) /SMD YNM 
ENMPCG=(ENMCHG*100. ) /SEDYNM 
RMNMPC=MNMPCG+0 . 00 
RENMPC=ENMPCG+0.00 

MNETCG=( (100.*SMACC)/SDMACC)/( (l.*SMEQNM)/SMDYNM)  - 100. 
ENETCG=( (100.*SEACC)/SCEACC)/( ( 1 . *SEE GNM > /SED YNM ) - 100. 
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RMNETC=MNETCG+0.00 

RENETC=ENETCG+0.00 

WRITE (6,2649) YEAR .STATE .RMPCTC .RMNMPC .RMNETC .REPCTC .RENMPC .RENETC 
2649  FORMAT(*1*/////40X.*  NET  CHANGE  IN  *.14^,*  ACCIDENTS  FOR  • ,4A4  , A2 

1 /////IX .40X.F6.1 . •%  CHANGE  IN  MORNING  ACCIDENTS*/ 

2 1X.40X.F6.1. •%  CHANGE  IN  NI GH  T + MiIDD  A Y ACCIDENTS*// 

3 1X,40X,F6.1, •%  NET  CHANGE  IN  MORNING  ACCIDENTS*//// 

4 1X.40X.F6.1, •%  CHANGE  IN  EVENING  ACCIDENTS*/ 

5 IX ,40X ,F6.1 , '%  CHANGE  IN  MIDDAY+NIGHT  ACCIDENTS*// 

6 1X,40X,F6.1,*X  NET  CHANGE  IN  EVENING  ACCIDENTS*) 

WRITE (6. 2639) YEAR, ST ATE 

2639  FORMAT) *1*/////1X,30X, • COMPARISON  OF  *.I4.*  ACCIDENT  TOTALS*, 

1 • FOR  • ,4A4,A2) 

WRITE  (6,2651  ) 

2651  FORMAT (30X,*  USING  ACTUAL  PERCENT  ACCIDENT  CHANGES*) 
WRITL(6,2646) 

2648  FORMAT(///50X. '6-MONTH  O.S.T.  FULL-YEAR  D.S.T.*/) 
M6MTH=SMN0N+SMDST 
E6MTH=SEN0N+SEDST 

MFYEAR= (MPCTCG/100  . + l.)*SMNON  + SMDST 
EFYEAR=(EPCTCG/100.  + l.)*SENON  + SEDST 
T0T6M=M6MTH4SN00N+E6MTH+SMIDN 
TOTFYR  = MFYEAf<  + SNOON  + FFYEAR+Sr*IDN 

WRITE (6,2640 ) M6MTH , MF YE AR , SNOON , SNOON , E6MTH , EFYEAR , SMIDN , 

1 SMIDN, TOTbN.TOTFYR 

2640  FORMAT ( 30X , * MORNING  ACC  I DENTS : • , 4X , 17 , 1 2X , 17/ 

1 30 X , * MID-DAY  ACCIDENTS: • , 4X . 1 7 , 12X , 1 7/ 

2 30X  , • EVENING  ACC  I DENTS : * . 4X , 1 7 , 12X , 17/ 

3 30X  , ♦ NIGHT  ACCIDENTS:*, 4X, 17, IPX, 17/// 

4 30X , ' TOTAL* *17x, I7.12X.I7) 

YCHNGE=T0TFYR-T0T6M 

YPCTCG=(100.*YCHNGE)/TOT6M 

RYPCTC=YPCTCG+0.000 

WRITE (6 ,2645) YCHNGE .RYPCTC 

2645  FORMAT (////30X, • CHANGE  IN  YEARLY  ACCIDENT  TOTAL  =*,I7/// 

1 30X , • PERCFNT  CHANGE  IN  YEARLY  ACCIDENTS  =*,F7.2,*X*) 

WRITE (6. 2639) YEAR, ST ATE 
WPITL(6,2647) 

2647  F ORMAT ( 30X , * USING  NET  PERCENT  ACCIDENT  CHANGES*) 

WR ITE  < 6 ,2648  ) 

MMTFYR= (MNETCG/100.  + l.)*SMNON  + SMDST 
ENTFYR=(ENETCG/100.  + l.)*SENON  + SEDST 
TNTFYR=NNTFYR+SNOON+FNTFYR+SMIDN 

WRITE (6 ,264  0 )M6MTH,MNTFYR ,SN00N,SN00N,E6MTH,ENTFYR,SMIDN, 

1 SMIDN,T0T6M,TNTFYR 
YNTCHG=TNTFYR-T0T6M 
YNTPCG= ( 100.*YNTCHG ) /T0T6M 
RYNTPC=YNTPCG+0.000 
WRITE (6.2645) YNTCHG . R YNTPC 
C 

C HOURLY  ACCIDENTS  FOR  EVENING  QUARTER-DAY 
C 

WR ITE ( 6 , 2655  ) YEAR, STATE 

2655  FORMAT ( 1H1//1X. 14, * HOURLY  EVENING  ACCIDENTS  FOR  • , 


1 4A4.A2) 

WR ITE  1 6 . 2656 ) 
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c 

2436 

C 

C 

2437 

2438 
C 

2737 

2439 
2730 

2653 

2739 

2435 

2654 
2658 


2 2x«  *6-10 ‘.Sx tlX. •FOOAY* ,lX.3x.'4-5' ,3X. '5-6* ,3X. * 

3 3x , •7-8’ t3Xi *8-9*  , 2 x » ’9-10*  , 5x  < 3X « *4-6' » 2X , *6-10'/ 

DO  2435  I=1,MDAYS 
0AY=EGDAY ( 1,1) 

EEQDAY=EQDAY(I,3) 

MHDAY=MNTHNO (DAY  ) 

DYDAY=DAY-DYBFM2(MHDAY+1 ) 

DWKDAY=WEEKDY (YEAR, DAY) 

MHEEQ=MNTHNO(EEQCAY) 

DTEEG=EEQDAY-DYBFM2(MHEE0+1) 

DWKEEQ=WEEKDY(YEAR,EFQ0AY) 

IF  (OAY.LT.l  .or.  day.  GT. 366 
2 .OR.  EEQDAY.LT. 0 .OR.  EEQDA Y .GT .366 ) GOTO  2739 

IF  (.NOT.LFAP  .AND.  (DAY  .EQ .60 
1 .OR.  EEQDAY.Efl.60 ) ) GOTO  2739 

DO  2436  J=l,6 

EVACI J)=NEVEAC(DAY.J) 

CONTINUE 

EVAC(7)=T0T‘+6(DAY) 

EVAC(8)=T0T610(DAY) 


IF  (EEGDAY.EQ.O (GOTO  2737 
EQUIVALENT  EVENING 
DO  2437  J=l,6 

EQEVAC( J)=NEVEAC(EEQOAY. J) 

CONTINUE 

EGEVAC (7)=T0T46(EEQDAY) 

EQEVAC ( 8 )=TOT610 (EEQDAY  > 

00  2438  J=l,8 

SEVAC ( J ) =SF  V AC ( J ) +E V AC (J) 

SEGEVA ( J)=SEGEVA ( J ) +FQE V AC ( J ) 

CONTINUE 
GOTO  2738 

NO  EQUIVALENT  EVENING 
DO  2439  J = 1 , 8 
EVAC ( J)=-EVAC ( J) 

EQEvAC ( J )=0 
CONTINUE 

WRITE  <6 ,2653  )MHD AY. DYO AY. E VAC. MHEEO.DYEEQ, EQEVAC 
FORMAT (IX, 12, */• ,12, • I', IX, 616, 5X, 216, 

1 7X  ,1X  ,12,  •/•  ,12, * ! • ,1X  »6I6,5X ,216 ) 

GOTO  2435 
NERROR=NERROR+1 
CONTINUE 
WRITE (6 ,2604  ) 

WRITE (6 ,2654 ) SEVAC ,SFOEV A 

FORMAT (1X,6X,1X,6I6,5X,2I6,7X,1X,6X,1X, 


L EI6.5X.2I6) 

WRITE (6, 2658) YEAR, ST ATE 
F0RMAT(1H1//1X,I‘».  • EVENING 
L 'SUMMARY  FOR  *,4A4,A2/ 

> //1X.10X.5X. • ACTUAL  ' , 5X. 'EGUI V 

’ , 8 X . ' NET  •/ 

* ,5X, 'ACCIDENTS  ' . 5X . • ACC IDENTS  • 
• 8X , ' PERCENT  '.5X.*  +NIGHT  *,0X 
TOTALS  ’ , 5 X . • TOTALS 
' , 5 X , • CHANGE  *.8X,' 


INJURY-CAUSING  ACCIDENT 
DAY  '.8X.10X.8X 


5X 

IX 


MIDDAY 
HOURS 

• change 

1X.10X.5X, • 

' CHANGE 


i 8X 

, » 


♦ ,8X.10X.8X< 
CHANGE’//) 


6-7  • . . 

> 


,'  APPARENT 
, 

PERCENT'/ 
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DO  2440  1=1.8 
IF  <I.GE.7)WRITE<6.2659) 

2659  FORMAT!/) 

OLEVAC ( I ) = S E Q E V A ( I ) - S E V A C ( I ) ^ 

PCEVAC ( I ) = (OLEVAC ( I ) * 1 0 0 . ) /SE V AC (I) 

NCGEVA ( I )=( ( 100.*SECFVA( I > )/SEVAC(I> )/( ( 1 . *SEEONK ) /SEO YNM ) - 100. 

TIMl=I+3 

TIM2=I+4 

IF ( I .EQ.7)TIM1=4 
IF ( I . EG.7)TIM2"6 
IF  ( I ,EQ.8)TIM1=6 
IF  (I. EQ. 8)11  M 2 = 10 

WRITE  <6 ,2657 )TIM1 .TIM2.SEVAC  (I ) .SEOEVA  (I  ) , OLEVAC (I  ) .PCEVAC  (I  ) . 

1 ENMPCG.MCGEVA ( I ) 

2657  FORMAT ( IX. Il . 'PM-' . 12. *Pm:  ’,5X.l7.3X.5X.l7,3X,8X.I7.3X,eX, 

1 F7 .2, ,2x.5X»F7.2. •*' .2X.8X.F7.2, * % * ) 

2440  CONTI NUF 

WRITL(6,2660) 

2660  FORMAK////) 

TIM1=4 

WRITE (6. 2657) TIM1 .TIM2, SDEACC.SE ACC.ECHNG.EPCTCG.ENMPCG.ENETCG 

END 
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A, 4 THE  DOUBLE  DST  EQUIVALENT  DAY  PROGRAM/  DBLDAY 


79 


C OBLOAY 
C 

C FINDS  THt  DOUBLE  D.S.T.  EQUIVALENT  DAY  FOR  EACH  DAY 
C 

C THE  DOUBLE  D.S.T.  EQUIVALENT  DAY  IS  CHOSE'N  TO  BE  ON  THE 

C SAME  DAY  OF  THE  WFEK  AS  THE  GIVEN  DAY. 

C HOLIDAYS  AND  THE  NIGHTS  BEFORE  HOLIDAYS  ARE  NOT  USED  AS  DAYS 
C OR  AS  EQUIVALENT  DAYS. 

C NO  EQUIVALENT  MORN  I NG/E VEN ING  IS  USED  MORE  THAN  FOUR  TIMES. 

C 

C 

IMPLICIT  INTEGER  (A-7) 

DIMENSION  MTHBGN ( 12 ) • MTHEND ( 12 ) 

DATA  MTHBGN/ 1.32. 61, 92, 122. 153, 183, 214, 245. 275, 

1 306,336/ 

DATA  MTHENL/31 .60,91.121,152,182,213,244,274,305.335. 

1366/ 

DIMENSION  R T A T E ( 5 ) ,HOL (50 ) ,M0BLDY(366) , EDBLDY ( 3661 
DIMENSION  RISE<366).  SET<366>.  RUSAGE<366).  SUSAGE(366> 
LOGICAL  LEAP 

DATA  MDBLDY .EOBLCY.RIIS AGE , SUS AGE/366* 0 , 366*0. 366*0. 366*0/ 
READ  (8.601)  STATE. YEAR, NHOL 

601  FORMAT  ( 4 A4 , A2 , IX , 1 4 . 1 X , I 2 ) 

WKITL(6,691)  STATE, YEAR, NHCL 

691  FORMAT  (1H1,4A4,A2,1X,I4.1X.I2> 

READ  (8.602)  ( HOL ( I ) . I = 1 . NHOL ) 

602  FORMAT  (20<l3.1X)) 

WRITL(6,692)  ( HOL  ( I ) . 1 = 1 ,NH0L  ) 

692  FORM  AT  ( 1 X . 20 ( 1 3 , 1 X ) ) 

00  401  K=1 , 12 

M1  = M THBGN ( M ) 

M2  = MTIIEND  (M  ) 

READ  (8.603)  (RISE(l).  I=M1.M2) 

603  FORMAT  (16(14, IX)) 

WRITE(6,693)  (RISE! I ) , I = M1,M2) 

693  FORMAT  (1611X.I4)) 

401  CONTINUE 

DO  402  M=1 , 12 
M1  = MT  HBGN ( M ) 

M2=MTHEND(M) 

READ  (8,603)  (SET(I).  I=M1,M2) 

WRITE (6,693)  (SET(I),  I=M1,M2) 

402  CONTINUE 
C CHECK  YEAR 

IF  ( YEAR.EG.1972)G0  TO  1972 
IF  ( YEAR. EQ. 1971 >GO  TO  1971 
1972  LE AP= . TRUE . 

DSTBGM=121 
DSTEnD=302 
60  TO  701 

1971  LEAP=.FALSF . 

DSTBGN=llo 
DSTEND=304 
GO  TU  701 

C ADD  ONE  HOUR  FOR  DST  DAYS 
701  DO  499  I=DSTBGN,CSTEND 
RISE ( I )=RISE( I ) + 100 
SET  (I  )=SET(I )+100 
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499  CONTINUE 

IF  (LEAP)  60  TO  703 
C NON-LEAP  YEAR 

C DELETE  FEE  29  * » 

R ISE ( 60 ) =0 
SET ( 60 ) =0 

GO  TO  703  . .. 

C LEAP  YEAR  AND  NON-LEAP  YEAR 
C REMOVE  HOLIDAYS  AND  NIGHT  BEFORE  HOLIOAYS 
703  DO  403  I =1 . NHOL 
HOLDAY  = HOL ( I ) 

RISE  ( HOLDA Y ) =0 

SET ( HCLDA Y ) =0 

IF  (HOLDAY.NE.l ) GOTO  702 

SE  T ( 566 ) =0 

GOTO  403 

702  SET ( HOLDAY-1 ) =0 
403  CONTINUE 

C FIND  DOUBLE  D.S.T.  EQUIVALENT  DAYS 
DO  4 o 9 1 = 1,366 

IF  (.NOT.  LEAP  .AND.  I.EQ.60)  GOTO  409 

MDBLOY ( I )=DBLEQV (RISF  , I ,RUSAGF  ,LEAP,DSTBGN,DSTENO) 

EDBLGY ( I )=CBLE3V( SET. I , SUS AGE . LEAP . OSTRGN . DSTEND > 

409  CONTINUE 
C PRINT  OUTPUT 

WRITE (9. 6 17) YEAR, ST ATE 

617  FORMAT (IX, 14,  ' 00U6LF  D.S.T.  EQUIVALENT  DAYS  FOR  ’.4A4.A2//) 
WRITE (6, 60 7 (YEAR. ST ATE 

607  FORMAT  (1H1.I4,*  DOUBLE  D.S.T.  EQUIVALENT  DAYS  FOR  »,4A4,A2//> 
WRITE (9,604  ) 

WRITE(6,fe04  ) 

604  FORMAT  ( 2X .• DA Y •, 7X ,• MORN ING 2X »• EVEN ING • /12X ,• DBEQDA Y • . 

1 2X, 'DREGDAY' ) 

DO  411  1=1,366 

IT  (I.EQ.60.AND..NOT.LEAP)GOTO  411 
M=MDBLD Y ( I ) 

E = ED3L0  Y ( I ) 

WRITE19, 614)1,  Ml,  E 
WRITE (6, 614)1, M,E 

614  F ORM AT ( 2X , I 3 « •:»,8X,I3,6X,I3) 

411  CONTINUE 

WR ITE ( 6 , 607 ) 

WR I TE ( 6 , 6 04  ) 

DO  410  1=1.366 

IF  (I.EQ.60.AND..NOT.LEAP)G0  TO  410 
J = MDBLD Y ( I ) 

K=ED6L0 Y ( I ) 

Ml=MNTHNO ( I ) 

N1  = D A YNO ( I ) 

M2=MNTHN0(u) 

N2=D A YNO ( J ) 

M3  = M.NTHN0  ( K ) 

N3=DAYNO(K ) 

WRITE ( 6 , 60S ) M1.N1,M2,N2,M3.N3 

605  FORMAT  (1H  , 1 2 , 1H/ , 1 2 , 7X » 12 . 1H/ , 12 , 4X . 12 . 1H/ . 12 ) 

410  CONTINUE 
END 
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C FIND  MORNING/EVENING  DOUBLE  O.S.T.  EQUIVALENT  DAY 

INTEGER  FUNCTION  DELF Q V ( R ISE T , I , US AGE . LE AP , DSTBGN , DSTFNO ) 
IMPLICIT  INTEGER ( A-Z ) 

LOGICAL  LEAP 

DIMENSION  RISE T ( 366 ) .USAGE ( 366 ) 

IF  (RISET(I).EQ.O)  GOTO  706 
DYOFWK=I  - ( 1/7 ) *7 

IF  (.NOT.  LEAP  .AND.  I .GE . 60 ) DYOFWK=UYOFWK-1 
IF  (DYOFwK.LE.O)  DY0FWK=0Y0FWK+7 
ODST I N = R I SF  T ( I )+200 

IF  (I.GE.DSTDGN  .AND.  I.LE.CSTEND)  DOST  I M = R I SET ( I ) +1 0 0 
DRLEQV=0 
M IND I F = 5 

DO  411  J=DY0FWK.366.7 
EQ  = J 

IF  (.NOT. LEAP  .AND.  J . GE . 6 0 ) EQ  = J + l 
IF  (LG.GT .366IG0T0411 
IF  (USAGE ( E Q ) . GT . 3 ) GOTO  411 
DIFF=GETDIF (DDSTIM.RISET (EQ) ) 

IF  (UlFF.GT .MIN DIF)  GO  TO  411 
IF  (OIFF.LT.MNDIF)  GO  TO  705 
IF  (APART (I. EG) . GT . APART ( I . DBLEQV ) JGO  TO  411 

705  DBLEuV=EU 
MINDIE  = D I F F 

411  CONTINUE 

IF  (UELEQV.GT.O)  USAGE < DBLEGV > =USAGE ( DBLEQV ) + 1 

RETURN 

706  D6LEGV=0 
RF  TURN 
END 


C GET  TIME  DIFFERENCE  IN  MINUTES  BETWEEN  TWO  GIVEN  TIMES 
INTEGER  FUNCTION  GETDIF  ( CDSTIM , JT IME ) 

INTEGER  TIME1.TIME2.H0UR1.H0UR2.DDSTIM 
IF (CDSTIM. LT .JTIMEJGO  TO  712 
TIME1=JTIMF 
T I ME2=CDST I M 
GO  TO  711 

712  TIME1=DDSTIM 
TIKE2=JTIME 
711  H0UR1=TIME1 /100 
H0URi;  = TIME2/100 
GETDIF=(HOUR2-HOUR1)*60 
.MINUT1  = TIME1-H0UR1*1()0 
MINUT2=TIMF2-HOUR2*100 
GETD1F=GETDIF+MINUT2-MINUT1 
IF (GETDIF .GT. 720 ) GET0IF=1440-GETDIF 
RETURN 
END 
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C FINO  HOW  FAR  APART  TWO  DAYS  ARE  t RASED  ON  WEATHER) 
INTEGER  FUNCTION  APART  (DAY1.DAY2) 

INTEGER  DA Y 1 1 D A Y2 
DATA  PI/3.14159/ 

DY1=DAY1 

DY2=DAY2 

WETHRl=COS (2.0*PI*(DYl-32.0)/366>.0) 
WETHR2=COS(2.0*PI*(DY2-32.0)/366.0) 

APART  = ABS(WETHR2- WE  THRU *1000.0 
IF  (UAY2.EO. 0 ) APART=5000 
RETURN 
ENO 


C F IND  MONTH  NUMBER  FROM  DAYCODE 
FUNCTION  MNTHNO ( 0 AYCOE ) 

INTEGER  DAYCDE 
IF(DAYCDE.EO.O)MNTHNO=0 

IF(DAYCDE.GT.0.AND.DAYCDE.LE.31 )MNTHNO=l 
IF (DAYCDE.GT.31 . AND . BA YCDE . LE . 60 )MNTHNO=2 
IF (DAYCDE.GT .60.AND.nAYCDE.LE.91 )MNThN0=3 
IF(OrtYCDE.GT. 91. AND. DAYCDE. LE.121)MNTHN0=4 
IF (DAYCDE.GT. 121 .AND. DA YCDE.LE. 152) MNTHN0=5 
IF (DAYCDE.GT. 152. AND. DAYCDE. LE. 182) MNTHN0=6 
IF(DAYCUE .GT.182. AND . C A yCDE . IE . 2 1 3 ) MNTHN0=7 
IF (DAYCDE.GT.213. AN0.DAYCDE.LE.244 )MNTHN0=8 
IF (DAYCDE.GT, 244. AND. DAYCDE. LE. 274) MNTHN0=9 
IF ( DAYCDE.GT. 274. AND. DAYCDE. LE. 305 )MNTHNO=10 
IF (DA YCDE. GT. 305. AND.D A YCDE.LE. 335) MNTHNOrll 
IF (DAYCDE .GT.335 )MNTHNO=12 
RETURN 
END 


C FIND  DAY  NUMBER  FROM  DAYCODE 

INTEGER  FUNCTION  DAYNO ( DAYCDE ) 

INTEGER  DAYCDE 

IF(DAYCDE.LE.31)DAYNO=DAYCDE 

IF (DAYCDE.GT. 31. AND . I)  A YCDE . LE . 60 ) D A YNO=DA YCDE-31 
IF(CAYCDE.GT.60. AND . DAYCDE . L E . 91 ) DA YNO=DAYCDE-60 
IF (CAYCDE.GT. 91. AND. DAYCDE. LE. 121 )DAYNO=0AYCDE-91 
IF (DAYCDE.GT. 121 . ANO.C A YCDE.LE. 152) DAYN0=DAYCDE-121 
IF (DAYCDE.GT. 152. AND . D A YCDE . LE . 182 ) D A YNO=D A YCDE- 1 52 
IF (DAYCDE.GT. 1«2. AND. DAYCDE. LF .213)DAYNO=DAYCDE-182 
IF (DAYCDE.GT .213. AND. DA YCDE.LE. 244 )DAYN0=DAYCDE-213 
IF (DAYCDE.GT. 244. AND . D A YCDE . L E . 274 ) DA YN0=DAYCDF-244 
IF (DAYCDE.GT.274. AND. DA YCDE.LE. 305 > C A YN0=DA YCDE-274 
IF (DAYCDE.GT. 305. ANO . C A YCDE . L E . 335 ) D A YNO=DA YCDE-305 
IF (DAYCDE. GT.335 ) DA YNC=DAYCDE-335 
RETURN 
END 
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A.  5 


THE  DOUBLE  DST  COMPARISON  PROGRAM, 


DBLCOMP 


n n n n n 


C DBLCOMP 
C 

C ROUTINE  TO  COMPARE  THE  ACCIDENT  TOTALS  OF 
C ALL  THE  OATS  OF  THE  YEAR  WITH  THOSE  OF  THEIR 
C DOUBLE-DAYLIGHT-SAVING  TIME  EQUIVALENT  DAYS.  > 

C 

C 

C ACTUAL  EQUIVALENT  DAY  ACCIDENT  TOTALS  ARE  USED. 

C 

C 

IMPLICIT  INTEGER  (A-7) 

REAL  MPCTCG.EPCTCG.YPCTCG 
REAL  MPCTC1 .EPCTC1 
REAL  MPCTC7 .EPCTC2 

REAL  MNMPCG,ENMPCG.MNETCG.ENETCG*YNTPCG 

REAL  MNKPC1  . ENMPC1  * M.NETC1  iENETCI 

REAL  MNMPC2 ,ENMPC2,MNETC2 ,ENETC2 

REAL  RMPCTC  ,1'EPCTC 

REAL  RMNMPC .KENMPC 

REAL  RMPCT1 .PEPCT1 

REAL  RMNMP1 .RENMP1 

REAL  RMPCT2 «FEPCT2 

REAL  RMNMP2  «RENM‘P2 

RFAL  RMNETC.RENETC , R YPCTC . R YNTPC 

REAL  R MNE  T 1 , RENET  1 

REAL  RMNE  T2 , RENE  T2 

REAL  QAVDWh , A VDWK . DMNRM , DE NRM , MNRM , ENRM . Q A VMNH , A VMNH 
LOGICAL  LEAP 
LOGICAL  OST 

DIMENSION  NFATAL<366.4) ,NCASUL (366,41 
DIMENSION  CSTBGN ( 2 ) , DSTEND (2  ) 

DIMENSION  S T A T E ( 5 ) 

DIMENSION  flWKTCT (7,4) ,DWKN0(7) ,QAVDWK(7,4) ,AVDWK(7) 

DIMENSION  MNHTOT (12,4) ,MNHN0(12) , Q A VMNH ( 1 2 , 4 ) , A VMNH (1? I 
DIMENSION  M NHALL ( 12  > 

DIMENSION  DYDFM2 ( 13 ) 

DIMENSION  EQLAY (366,3) 

DATA  LEAP/. FALSE./ 

DATA  SM ACC. SCMACC.SE ACC, SC EACC/0, 0,0.0/ 

DATA  SMAC1.SDMAC1 .SEA Cl. SDEAC 1/0 .0,0.0/ 

DATA  SMAC2.SDMAC2.SEAC2.SDEAC2/0.0.0.0/ 

DATA  DSTBGN/llb.121/ 

DATA  CSTLND/304,302/ 

DATA  SMNON.SENON.SMDST .SEDST .SNOON.SMIDN/O. 0.0. 0.0,0/ 

DATA  D WK TOT. DWK NO. MNHTOT, M|V UNO. MNH ALL /28* 0.7* 0,48*0,1 2* 0,12*0/ 
DATA  DYEFM2/0, 0.31.6D, 91, 121, 152 ,162,213,244,274,305,335/ 

DATA  NERROR/O/ 

OAT A SMD YNM,. SED YNM. SMEQNM, SEE ONM/O, 0,0,0/ 

DATA  SMDYN1 . SEDYN1 , SMEQN1 , SEEQN1/0 , 0 , 0 , 0/ 

DATA  SMDYN?,SEOYN2,SMEQN2,SEEQN2/0, 0,0,0/ 

DATA  NMEODS.NEEGDS/O.O/ 

DATA  NMEQD1 ,NEEGCl/0 .0/ 

DATA  NMEODP ,NEECC2/0 . 0/ 


GET  STATE  ANO  YEAR 
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REAO(10,2630)YFAR, STATE 

2630  FORMAT ( IX . I 4 , 35X . 4 A 4 . A 2//// ) 

IF  (YtAR.EG.1972)  LE  AP  = .TRUE . 

C * ' 

C READ  QUARTER-DAY  ACCIDENT  TOTALS 
C 

READ  (8 .2601 ) ( (NFATALI I . J) . J=1.4) .1  = 1.366) 

2601  FORMAT (/////////// (7x.4( lx. 19) ) ) 

READ (8 ,2641 )NFACC .NFTLTY 

2641  FORMAT (//28X.I7/28X.I7) 

READ (8 ,26  01  ) ( (NCASUL ( I , J ) . J = 1 , 4 ) , I = 1 . 366 ) 

READ(8 .2642 ) NIN J AC .NTNJRD 

2642  FORMAT (//39X, I7/39X, 17 ) 

IF  (LEAP)  GOTO  2720 

NFATAL(59,4)=NFATAL(S9,4)+NF ATAL(60.4» 

NCASUL  (59.4  )=NCASUL  (59 .4  ) + NIC  A SUL  (60,4) 

2720  WRITE (6,2631 ) YEAR .STATE 

2631  FORMAT! •1,,10X,I4.1X.4a4.a2) 

WRITE (6, 2610 ) ( I , ( NF  AT  AL ( I » J ) . J=1 ,4) , 1 = 1.366) 

2610  FORMAT ( ' NUMBER  OF  ACCIDENTS  WITH  F AT AL I T IES • /// ( IX , 1 3 . 

1 3X  < 4 ( 1 x , 1 9 ) ) ) 

WRITE (6 .2643 )NF ACC .NFTLTY 

2643  FORMAT!//'  NUMBER  OF  FATAL  ACCIDENTS  ='.17/ 

1 • NUMBER  OF  FATALITIES  =',I7) 

WRITE (6 ,2631 (YEAR .STATE 

WRITE  (6,2611 )( I , (NCASUL ( I , J) , J=1 ,4 ) ,1=1,366) 

2611  FORMAT! ' NUMBER  OF  ACCIDENTS  WITH  ANY  INJURY’/// 

1 ( IX. 13 , 3X .4 ( lx . I 9 ) ) ) 

WRITE (6,2644)NINJAC,NINJRD 

2644  FORMAT ( //  ' NUMBER  OF  I N JUR Y-C AUS ING  ACCIDENTS  =',17/ 

1 • NUMBER  OF  INJURED  (INCL.  FATALITIES)  =*,I7) 

C 

C GET  DAY-OF-WEEK  AND  MONTHLY  TOTALS  AND  AVERAGES 
C 

DO  2420  1=1,366 

IF  (.NOT.LFAP  .AND.  I.EQ.60)GOTO  2420 
DYOFnK=WEEKDY (YEAR, I ) 

MHOFYR=MNTHNO ( I ) 

DO  2421  J = 1 , 4 

DWKTOT(DYOFWK,J)=DWKTOT(DYOFWK,J)  + NCASUL(I.J) 

MNHTOT( MHOFYR. J ) =MNHTOT ( MHOF YR , J)+NCASUL( I. J) 

MNHAlL (MHOFYR ) =MNH ALL ( MHOF YR ) +NC ASUL ( I , J) 

2421  CONTINUE 

DWKN<J(DYOFWK>=DWKNO(DYOFWK)  + l 
MNHNO ( MHOFYR  >=MNHNO ( MHOFYR  > + 1 
2420  CONTINUE 

WRITE (6,2631 ) YEAR, ST ATE 

WRITE (6,2620 ) (D, (DwKTOTID. J) . J=1 .4 ) . DWKNO ( D ) ,D=1.7) 

2620  FORMAT ( ’ ACCIDENT  TOTALS  FOR  EACH  DAY  OF  THE  WEEK’// 

1 (3X,I2, • : ' ,4 ( IX, 19) ,3X. 'NO.DAYS=' ,19)  ) 

WRITE (6 ,2631 (YEAR .STATE 

WRITE16.2621)  (M,  (MNHTCT  (M.J)  ,~J=1.4)  . MNhALL  ( M ).  MNHNO  ( M ).  M = 1 , 12  ) 

2621  FORMAT! ’ ACCIDENT  TOTALS  FOR  EACH  MONTH’// 

1 (3x.I2.*: ’.4(lX»I9 ) ,3X. ’TOT=’ .I9.3X. ’NO.DAYS=* .19  ) ) 

DO  2426  D = 1 , 7 
DO  2428  J=1 ,4 

QAVDwK (D,J)=(l.*DWKTOT (D,J) ) /DWKNO (D) 

2428  CONTINUE 
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SIGCWK=DWKTOT(D. ll+DWKTOT (D,2)+DWKT0T(D,3)+DWKT0T(D.41 
AVDWK(C)=(l.*SIGCWK)/DWKNO!D) 

2426  CONTINUE 

DO  2427  14  = 1,12 
DO  2429  J=1 ,4 

QA VMNH  (Ml.  J)=<l.*MNHTOT(M,J)  )/MNhNO(M> 

2429  CONTINUE 

SIGMNH=MNHT0T<M.1)+MI\IHT0T(M.2)  + MNHT0T(M.3)+MNHT0T<M.4) 

AVMNHt  M ) = ( 1 . + SIGMNH ) /MNHNO  t M ) 

2427  CONTINUE 

WRITE (6, 2631 ) YEAR. STATE 

WRITE(6,2622)(0,(QAvnwK(0.U).J=l,4>,AVDWK(D).D=1.7> 

2622  FORMAT ( 1H  ,•  AVERAGE  ACCIDENTS  FOR  EACH  DAT  OF  THE  WFEK • // 

1 (3X.I2,':',4(lX,F9.1>.9X.'AV='.F9.1>) 

WRITE (6. 2631) YEAR, STATE 

WRITE (6,2623)  (Mi  (OA VMNH (M , J) , J = 1 ,4 ) , A VMNH ( M ) ,M  = 1 ,12  ) 

2623  FORMAT ( 1 H ,»  AVERAGE  ACCICENTS  FOR  EACH  MONTH*// 

1 (3X ,12,  • I • ,4 (1X.F9.1)  ,9X,  'AV=-  .F9.1  ) ) 

C 

C 

C READ  EQUIVALENT  DAYS 
C 

1 = 0 

2703  1=1+1 

READ (10, 2605, END=2704 ) (EQOAY ( I , J ) , J=1 .3) 

2605  FORMAT  < 2X , 13, 9X , 13. 6X , 13) 

GOTO  2703 

2704  ND A YS= I -1 
WRITE(6,2632)YEAR, STATE 

2632  FORMAT ( • 1 ' ,14,  ' DOUBLE  D.S.T.  EQUIVALENT  DAYS  FOR  • ,4A4.A2// 

1 2X, 'DAY', 7X.’M0RNING'.2X, 'EVENING'/ 

2 12X , 'CBEQDAY • ,2X , *CBEQDAY • ) 

DO  2433  1=1 , NDAYS 
WRITE(6,2633)  (EGDAY(I.J) ,J=1,3) 

2633  FORMAT (2X , 1 3 , • : • , 8X , I 3 , 6X  , 1 3 ) 

2433  CONTINUE 
C 

C GET  MORNING/ EVFNING  ACCIDENT  TOTALS  FOR  D.S.T.  AND  NON-O.S.T. 

C 

DO  2434  1=1,366 

IF  (.NOT. LEAP  .AND.  T.EQ.60)  GOTO  2434 

IF  (l.GE.OSTBGN (YEAR-1970)  .AND.  I .LE .DSTEND ( YEAR-1970  ) )GOTO  2721 
SMNON  = SMNON+NCASUL ( I .1  ) 

SNOON=SNOON+NCASUL ( 1.2) 

SENONrSENON+MCASULt I .3) 

SMIDN=SMIDN+NCASUL( 1.4 ) 

GOTO  2434 

2721  SMDST=SMOST+NCASUL ( I « 1 ) 

SNOON=SnOON+NCASUL( I « 2 ) 

SEDST=SEDST+NCASUL(I«3) 

SMIDN=SMIDN+NCASUL (1.4) 

2434  CONTINUE 
C 

C GET  DAY'S  ACCIDENTS  VS.  FQUIVAl  T-DAY'S  ACCIDENTS 
C 

WRITl(6,2606)YEAR, STATE 

2606  FORMAT ( '1 ' ,13X, 'COMPARISION  OF  ',14,'  ACTUAL  AND 
1 'DOUBLE  D.S.T.  EQIJI V AL ENT -D A Y ' , 
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2 • INJURY-CAUSING  ACCIDENTS  FOR  ’,4a4,A2) 

WR ITE ( 6 , 2652 ) 

2652  FORMAT (///IX, ’DAY  DATE  *************  ACTUAL  ACCIDENTS’. 

1 »*************., 

2 10X. • ' ***********DOUfciLE  D.S.T.  EQUIVALENT  DAY 

3 • ACCILENTS************ ’ ) 

WRITE  ( 6 . 2b(j8  ) 

2608  FORMAT ( IX  .1 IX • ’MORNING*  «3X. ’EVENING A .7X. ’NIGHT+’ »2X. ’MIDDAY’ . 

1 1 3 X , ' — LQUIV  MORNING--’ .3X. ’--EQUIV  EVENING — ’,9X. 

2 * N I GHT  + MIDDAY’/ 

3 1 X » 9X  * * (4AM-10AM)  • .IX.  * (4PM-10PM) • »5X. 

4 'MIDDAY  +NIGHT ’ . 13X . • DAY  DATE  ACCDNTS’.3X. 

5 ’DAY  CATE  ACCDinITS  ’.  9X .'  MIDDAY  +NIGHT'///) 

00  2404  I=1,NDAYS 

DA Y=EGDAY ( I . 1 ) 

MEQD AY=EQDA Y ( I ,2) 

EEGDAY=EQDAY ( 1,3) 

MHDAY=MNTHNO (DAY ) 

DYDA  Y = DAY-DY8FNi2  ( MHDAY  + 1 ) 

DWKCAY=WEEKDY( YEAR, DAY) 

MHMEQ=MNTHNO (MECCA Y) 

DYMEQ=MEU0AY-DYBFM2 (MHMEQ+1 ) 

DWKMLG=WEEKDY ( YEAR .MEQDAY ) 

MHEEU=MNTHNO(CEOCAY) 

DYEEQ=EEQDAY-DYBFM2 (MHEEQ+l ) 

DWKEEG=WEEKOY ( YEAR.EEQDAY) 

IF  (UAY.LT.l  .OR.  DAY. GT. 366 

1 .OR.  MEQDAY. LT. (I  .OR.  MEQD  A Y . GT  . 366 

2 .OR.  FEQDAY.LT. 0 .OP.  EEQDAY.GT.366)  GOTO  2709 
IF  (.NOT. LEAP  .AND.  (DAY.EQ.60  .OR.  MEQDAY. EO. 60 

1 .OR.  EEQDAY.EQ.60 ) ) GOTO  2709 

DST=.  FALSE. 

IF  (UAY.GE.DSTBGN(YEAR-3970)  .AND.  D A Y . LE  . DSTEND ( YE AR -1 970 ) ) 

1 JST=.TRUE. 

OYBEPR=DAY-l 
IF  (DAY.EQ.l )DYBEFR=366 

IF  (DAY.EQ.61  .AND.  .NOT. LEAP)  DYREFR=59 

MBEFUR-MEQDAY-l 

IF  (MEQDAY. EQ. 1 )MBEFOR=366 

IF  (MEQDAY. EG  .61  .AND.  .NOT. LEAP)  M6EF0R  = 59 
C 

DMACC=NCASOL(DAY ,1) 

DEACC=NCASUL(DAY,3) 

C 

IF  (MEQDAY. EG. OIGOTO  2705 
C MORNING  DOUbLE  D.S.T.  EQUIVALENT  DAY 
NMEGUS=NMEGOS+l 
IF(DST)  NMFGD1=NMEGD1 +1 
IF(.NOT.DST)  NMEGD2=NMEQD2+1 
MACC^NC ASUl (MEQDAY, 1 ) 

SMACC=SMACC+MACC 

SDMACC=SDMACC+DMACC 

IF(DbT)  SMAC1=SMAC1+MACC 

IF(DST)  SDMAC1=SCMAC1+DMACC 

I F ( . NOT  . DST  ) SMAC2=.SMAC2  + MACC 

IF  ( .NOT  . DST)  SCMAC2  = SDMAC2  + DMACC 

MEGNM.D  = NCASUL(MREF0R,4  J+NCASUL  (MEQDAY  ,2) 

MDYNhD=NCASUL(DYEEFR,4)+NCASUL(DAY,2) 
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GOTO  2706 

C NO  MORNING  DOUBLE  D.S.T.  EQUIVALENT  DAY 

2705  M ACC  = 0 
MDYNMD=0 

MEQNMU=0  ' 

DMACC=-DmACC 

C 

2706  IF  (EEqDAY.E(j.O)GOTO  2707 

C DOUBLE  D.S.T.  EUUIVALFNT  EVENING 
NEEGDS=NEEGDS+1 
IF  (UST)  NE EGDl=NEEOnl+l 
IF  (.NOT.  DST  ) NEEQFj2  = NEEQD2  + 1 
EACC^NC ASUl  ( EEqU A Y , 3 ) 

SEACC=SEACC+EACC 

SDEACC  = SGEACC+I)EACC 

IF  (DST ) SEAC1=SEAC1+FACC 

IF  (DST )S0EAC1=SDEAC1+CEACC 

IF  ( .NOT  .DST )SEAC2=SE AC2+EACC 

IF  ( .NOT. UST  )SDEAC2  = SI)EAC2  + DEACC 

EOYN«U=NCASUL ( DAY, 2 ) +NCASUL (DAY. 4 > 

EEGNMDrNCASUL < EE GDA Y . 2 ) +NC ASUL ( EEQD A Y . 4 ) 

GOTO  2706 

C NO  DOUBLE  D.S.T.  EQUIVALENT  EVENING 

2707  EACC=0 
EDYNMD=0 
EEQNND=0 
UEACC=-DE ACC 

C TOTALS 

2708  SMDYUM=SMDYNM+MDYNMD 
SNEQl\lM  = SMECNN  + MEGNMO 
SEDYNNi  = SEDYNN  + EDYNMD 
SEEUI\iK  = SEEGNN  + EEGNMO 
IF  (DST)  GOTO  2722 
SMDYN2=SMDYN2+MDYNMD 
SME0N2=SMEGN2+MEQNMD 
SEUYN2=SEDYN2+EDYNMD 
SEEGN2  = SEE(jN2  + EEGNMD 
GOTO  2723 

2722  SMDYN1=SMDYN1+M0YNMD 
SMEQN1=SMEGN1+MEQNM0 
SEDYN1=SE0YN1+EDYNMD 
SEEQl\il  = SEEONl  + EEGNMD 

2723  WRITE  ( 6 , 26 0 2 ) Di/JhD A Y , MHD A Y , DYD AY  .DM ACC  ,DE ACC  .MDYNMD .EDYNMO  • 

1 DWKMECJ  .NHMEG  ,DYMEQ  ,MACC  , 

2 OWKEEu  .NHEEG ,OYEEQ ,EACC .MEQNND.EEQNMD 

2602  FORM  ATI  IX  ,1X  ,11 ,1X,I?,*/,.I2,1H:,I7.3X,I7,8X.,<',I6,I6,2X, 

1 •>,i12X,1H(»I1»1X»I2,,/,,I2,2H:),I7»3X»1H(»I1,1X»I2,,/,» 

2 I2,•:)•,I7,ex.,<•.I6,I6,2X.•>,) 

GOTO  2404 

2709  NERROR  = NERROF'  + l 

2404  CONTINUE 

WRITE (6 ,2604  ) 

2604  FORMAT  (///IX,  • - — , 

1 « • 

2 • •) 

WRITE  (6 ,2683)  SDMAC1  ,S0EAC1  .SMDYNl.SEDYNl  .SN.AC1.SE  AC  1, 

1 SMEQN1 , SEEGN 1 

2653  FORMAT (//IX  , ’DST  PER  1 CD « , 1 6 , 3X , 1 7 , 8X , • < • , 1 6 , 1 6 , 2X , * > • , 12X , 
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1 10X.l7.3X,10X.l7.ex. •<* , 16, I6.2X. •>• ) 

WRITE  (6 ,2654)S0MiAC2,SDEAC2,SM0YN2  , SEDYN2 , SMAC2 , SE AC2  , 

1  SMEQN2.SEEGN2 

2654  FORMAT ( / IX  « 'REST  OF  YR  • , 16 , 3 X . 17  * 8X . • < ' . J6 , 16 . 2X . • > • , 12X . 

1 10X.I7,3X,10X,17.8X.,<'.I6.I6,2X,*>') 

WRITM6.2603)  SQMACC ,SCE ACC .SMDYNM. SEDYNM. SM ACC. SEA CC. 

1 SMEGNM .SEEGNM 

26  0 3 FORMAT (//IX , "TOTAL • , 4X , 1 7 , 3X , 17 , BX , • < • , 16 , 16 , 2X , • > • , 1?X , 10X , 17 , 3X , 
1 10X.I7.8X. '<’.I6.I6.2X. •>’) 

WRITE ( 6 .264  6 HJMEGDS .MMEGD1 .NMEQD2 .ME EGOS .NEEQOl .NEEQD? 

2646  FORMAT (///////////////1X.22X. 'NO.  OF  MORMIMG  EQUIVALENT  DAYS  = '. 

1 17. 5X. MOST  PERIOD  =',I5.'.  REST  OF  YEAR  =»,I5,*)* 

2 //1X.2PX. 'MO.  OF  EVENING  EQUIVALENT  DAYS  ='.I7, 

3 bX.MOST  PERIOD  =',I5,',  REST  OF  YEAR  =',I5,*)'> 

MCHNG=SMACC-SDMACC 

ECHNG=SEACC-SDEACC 

MNMCHG=SMEQNM-SMCYNM 

ENMCHG=SEEQNM-SECYNM 

KCHNG1=SMAC1-SDMC1 

ECHNG1=SEAC1-SUEAC1 

hNMCGl  = SMEQ!\ll-SMDYNl 

ENMCG1  = SEE(JN1-SECYN1 

MCHNt>2=SMAC2-SDNAC2 

ECHNG2=SEAC2-SDEAC2 

MNMCG2  = SME(>N2-SMDYN2 

ENMCG2=SEEqN2-SEDYN2 

WRITE (6 .2650 )MCHNG,MCHNG1.MCHNG2.ECHNG,ECHNG1.ECHNG2, 

1 MN'MCHG  .MNMCGl  .MNMCG2  .ENMCHG  .ENMCGl  .ENMCG2 

2650  FORMAT (////IX .22X, 'CHANGE  II  MORNING  ACCIDENTS  =',I7, 

1 5X.MDS1  PERIOD  =',I5.*.  REST  OF  YEAR  ='.  15.*)*/ 

2 IX. 22X. 'CHANGE  IN  EVENING  ACCIDENTS  =*.I7. 

3 5x . ' < UST  PERIOD  ='.I5.'.  REST  OF  YEAR  = ',  15.')'// 

4 IX, 22X. 'CHANGE  IN  NIGHT-MIDOAY  ACCIDENTS  =’.I7. 

5 5X.MUST  PERIOD  ='.I5.*.  REST  OF  YEAR  ='.  15,  • ) */ 

6 IX  ,22X CHANGE  IN  MIDDAY-NIGHT  ACCIOENTS  =',I7, 

7 5X.MDST  PERIOD  =',I5.'.  REST  OF  YEAR  ='.  15,*)') 

WRITE (6, 2638)YEAR, STATE 

2636  FORMAT ( 'l'/////40x,  • CHANGE  IN  THE  LEVEL  OF  ACCIDENTS*/ 

1 4 0 X , ' UNDER  THE  EFFECT  OF'/ 

2 40X. 'DOUBLE  DAYLIGHT-SAVING  ALL  YEAR'/ 

3 4 0 X , • FOR  ',14,'  FOR  *,4A4,A2///) 

MPCTCG=(MCHNG*100. J/SCMACC 
EPCTCG=(ECFNG*100. J/SCEACC 

RMPC IC=MPCTCG+0 . 0 
REPCTC=EPCTCG+0.0 
WRITE (6, 2637) RMPC TC.REPCTC 

2637  FORMAT (3X.40X. ' MORNING  ACCIDENTS:  '.F6.1.'**// 

1 3X.40X,  * EVENING  ACCIDENTS:  *,F6.1.'X') 

MNKPCG=(MNMCHG*100. ) /SMDYNM 
ENMPCG=(ENMCHG*100. ) /SEDYNM 
RMNMPC=MNMPC6+0 . 0 
RENMPC=ENMPCG+0.0 

MNETCG=< (100.*SMACC)/SOMACC)/( ( 1 . *SMEQNM ) /SMDYNM ) -1 0 0 . 

Ef:ETCG=  ( <100.*SEACC)/SDEACC )/< ( 1 . *SEE QNM ) /SEDYNM ) -1 0 0 . 

RMNETC  = M!NETCG  + 0.0 
RENETC=ENETCG+0.0 
WRITE (6, 2655) YEAR,  ST ATE 

2655  FORMAT ( 'l'////40X, • NET  CHANGE  IN  ',14,'  ACCIDENTS  FOR  '.4A4.A2 
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1 /4 0 X « • UNDER  DOUBLE-DAYLIGHT-SAVING  TIME • ) 

WRITt <6,2649)RMPCTC.KMNMPC.RMNETC.REPCTC,RENMPC.RENETC 

2649  FORMAT  ( ///// 1 X , 4 0 X , F6 . 1 , • 26  CHANGE  IN  MORNING  ACCIDENTS*/ 

2 1X.40X.F6.1  . '%  CHANGE  IN  NIGHT  + M, TODAY  ACCIDENTS*// 

3 1X,40X,F6.1,*%  NET  CHANGE  IN  MOKNTNG'  ACCIDENTS*//// 

4 1X«40X»F6.1, ’%  CHANGE  IN  EVENING  ACCIDENTS’/ 

5 IX, 40X.P6.lt •*  CHANGE  IN  MIDDAY  + NIGHT  ACCIDENTS’// 

6 IX ,40X ,f 6.1  , *%  NET  CHANGE  IN  EVENING  ACCIDENTS*) 
MPCTC1=(MCHNG1*100. )/SDMACl 
EPCTC1=(ECHNG1*100. )/SDEACl 

RMPCT1=MPCTC1+U. 0 
REPCTl=EPCTCl+0.0 
MNMPC1=(MNMCG 1*100. ) /SMDYN1 
ENMPC1= ( ENMCG1*100.>/SEOYN1 
RMNMHl=MNMPCl+0.0 
RE NMP1=ENMPC 1+0.0 

MNETL1=  ( (10  0 .*SMAC1 1/SDMAC1  )/(  { 1 . *SMELiNl ) /SMD  YN1  ) -1 00  . 

ENET  Cl= ( <1G0.*SEAC1)/SDEAC1 )/( ( 1 . *SE EQN1 ) /SED YN1 ) -1 0 0 . 

RMNET 1=MNETC1+0.0 
RENE  I l = ENETCl  + 0.0 
WRITE (6, 2655) YEAR, ST ATE 
WRITE (6,2656 ) 

2656  FORMAT (40X  , ’ FOR  THE  D.S.T.  PERIOD  OF  THE  YEAR*) 

WRITE (6.2649)RMPCT1.RMNMP1,RMNET1.REPCT1,RENMP1,RENET1 

MPCTC2=(MCHNG2*100.)/SDMAC2 

EPCTC2=(ECHNG2*100. ) /SDEAC2 

RMPCT2  = M:PCTC2  + 0.0 

REPcT2=EPCTC2+0.0 

MNMPC2= ( MNMCG2+100. ) /SMDYN2 

ENMPC2=(ENMCG2*100. J/SEDYN2 

RMIJMP2  = Mi\lMPC2  + 0 . 0 

RENKH2=ENMPC2+0.0 

MNETC2=((100.*SMAC2)/SDMAC2)/( ( 1 . *SME0N2 > /SMD YN2 > -10 0 . 
ENETC2=((100.*SEAC2)/SDEAC2)/( < 1 . *SEEGN2 ) /SE0YN2 ) -1 0 0 . 

RMNEI 2=MNETC2+0.0 
RENET2=EUETC2+0.0 
WRITL(6,2655)YEAR, STATE 
WRITE(6,2657) 

2657  F0RmAT(40X,*  FOR  THE  NON-D.S.T.  PERIOD  OF  THE  YEAR’) 
WRITE(6,2649)RmiPCT2,RMNMP2.RMNET2,REPCT2,RENMP2,RENET2 
WRITE (6,2639) YEAR, ST ATE 

2639  FORMAT ( • 1 ' /////l X , 30 X , • COMPARISON  OF  • ,14,  • ACCIDENT  TOTALS* , 

1 • FOR  * . 4 A4 . A2 ) 

WRI  I L ( 6,2651 ) 

2651  FORMAT ( 2 7 X « • USING  APPARENT  PERCENT  YEARLY  ACCIDENT  CHANGES* ) 
WRITE(6,2648) 

2648  FORMAT (///50X, ’6-MONTH  D.S.T.  FULL-YEAR  D.D.S.T. •/) 

Mi6MITH=SMIN0N  + SMDST 
E6M.TH=SEN0N  + SEDST 

MFYEAK=(MPCTCG/100.  + 1.)  * ( SMNON  + SMDST  ) 

EFYEAR=(EPCTCG/100.  + 1.)  * (SENON  + SEDST) 

T0T6Mi  = M6MiTH  + SN00N  + E6MTH  + SMICN 
TOTFYR=MFYFAk+SNOON+FFYEAR+SMIDN 

WRITE (6 ,2640  > M6MT H , MF Y E AR , SNOON , SNOON , E6MTH , EF YE AR , SMIDN , 

1 SMIDN, TOT 6M , TOTFYR 

2640  FORMAT (30X,*  MORNING  ACC  I CEN TS I * , 4 X , 1 7 , 12X  , 1 7/ 

1 30  X . * MID-DAY  ACCICENTS: * . 4 a , 17. 12X, 17/ 

2 30  X , • EVENING  ACCIDENTS: • ,4X , 17 ,12X. 17/ 
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3 30X  « ' NIGHT  ACCIDENTS;  • .4X.I7.12X,  17/// 

4 3 0 X « • TOTAL* , 17 x , 1 7 , 1 2X . 1 7 ) 

YCHNGE=TOTF YR-T0T6M 

YPCTCG= (10G.*YCHNGE ) /T0T6M  . , 

RYPCTC=YHCTCC+U.OOO 
WRITE<6,2645) YCHNGL.RYPCTC 

2645  FORMaT(////30X, • CHANGE  IN  YEARLY  ACCIDENT  TOTAL  =*.I7/// 

i 3 u x « ' percent  change  in  yearly  accidents  =’,f7.2,’x’> 

WRITE (6. 2b39) YEAR, ST ATE 
WRITE (6,2647) 

2647  F ORRAT ( 30X , * USING  NET  PERCENT  YEARLY  ACCIDFNT  CHANGES’) 
WRITE (6 ,2648 ) 

MNTF YR= (»Nf TCG/100 . + 1.)  * (SMNON  4 SKDST) 

ENTFYR= ( E N I TCG/100.  4 1.)  * (SENON  4 SFDST  ) 

TNTFYR  = V.NTf  YR  4SNCON4F  NTFYR4SMIDN 

WP ITL (6 ,2640  )MbKTH .MNTFYR , SNOON .SNOGN , F6MTH , ENTF YR  , SMI ON , 

1 SMIDN.TOTGM.TNTFYH 
YNTCHG  = TI\lTF  YK-T0T6M 
YNTPCG=(  luu  .»  YNTCHG  ) /T0T6M 
RYNTPC  = YNTI  CG40.000 
WRITE ( 6 , 264  5 ) YNTCHG, KY NTPC 
WRITE (6, 2639) YEAR, ST ATE 
WRITE (6 ,2b61  ) 

2661  FORRAT ( 20X  . • USING  APPARENT  PERCENT  ACCIDENT  CHANGES*. 

1 * FOR  D.S.T.  AND  NON-D.S.T.  PERIODS’) 

WRITE(6,2648) 

MFYEAK=(MPCTC2/100.  4 1.)  * SMNON  4 (NPCTC1/100.  4 1.)  * SMDST 
EFYE AR= ( EPCTC2/100 . 4 1.)  * SFNON  4 (EPCTC1/100.  4 1.)  * SEDST 
TOTFYR  = f'FYFAR4SNOON4FFYEAR4SM10N 

WRITE  (6, 264  0 )M6MTH.MF  YEAR,  SNOON,  SNOON,  F6MTH.  EFYE  AR,  SMI  ON. 

1 SNIDN.T0T6N.T0TFYR 

YCHNGE=TOTF YR-TOT6M 
YPCTCG=(100.*YCHNGE)/TOT6|w 
RYPCTC=YPCTCG4 0.000 
WRITE(6,2645)YCHNGE,RYPCTC 
WRITE (6, 2639) YEAR, ST ATE 
WRITE (6 ,2667 ) 

2667  FORMAT (2UX, ' USING  NET  PERCENT  YEARLY  ACCIDENT  CHANGES’, 

1 ’ FOR  D.S.T.  AND  NON-D.S.T.  PERIODS’) 

WR I TE ( 6 , 2648 ) 

MNTFYR=(MNFTC2/100.  4 l.)*SNNON  + (MNETC1/100.  4 l.)*SMDST 
ENTFTR=(ENETC2/100.  4 1. ) *SE  NGN  4 (ENETC1/100.  4 l.)*SFDST 
TNTFYRrNlvTF  YR4SNOON4F  NTFYR4SMION 

WRITE (6,2640 ) M6MTH.NNTFYR , SNOON , SNOON . F 6MTH , ENTF YR , SM I ON . 

1 SMIDN.T0T6N.TNTFYR 
YNTCHG=TNTF YR-TCT6M 
YNTPCG= ( 100 .* YNTCHG ) /TOT 6M 
R YNTPC  = Y(mTPCO40 . 000 
WRITE (6,2645)YnTCHG,RYNTPC 
END 
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c 

c 

c 

C GET  MONTHNO 

C ' 1 

INTEGER  FUNCTION  MNTHNO(D) 

IMPLICIT  INTEGER  <A-7) 

DIMENSION  D YrFMH ( 12 ) 

DATA  DYBFMH/O, 31. 60, 91, 121. 152. 182, 213, 244, 274, 305. 335/ 
IF  (u.EQ.O)GOTO  2902 
DO  2901  MNTHNO=X.ll 
IF  (U.LE.DYBFMH(MNTHN0+1 ) ) RETURN 

2901  CONTINUE 
MNTHNO=12 
RETURN 

2902  MNTHNG=0 
RETURN 
END 


C 

C 

C GET  DAY  OF  WEEK 

C . , 

INTEGER  FUNCTION  WEEKDY ( Y ,D ) 

IMPLICIT  INTEGER  (A-7) 

LOGICAL  LEAP 
DIMENSION  UAN7 ( 2 ) 

DATA  JAN7/4 , 5/ 

DATA  LEAH/.rALSE./ 

IF  (U.LT.l)  GOTO  220? 

IF  (Y.E0.1972)  LE AP= .TRUE  . 

WEEKDY  = D - ( ( D/ 7 ) * 7 ) + JAN7<Y-1970> 

IF  (.NOT. LEAP  .AND.  n.GT.60t  WEEKDY=WEEKDY-1 

IF  (WEEKUY.GT.7)  WEEKCY=WEEKDY-7 

RETUhN 

2202  WEEKDY=0 
RETURN 
END 


94 


APPENDIX  B 

B.  1 


- OUTPUT 


SAMPLE  OUTPUT  OF  COMPAR 


1972  DISTRICT  COLUMBIA 
NUMBER  OF  ACCIDENTS  WITH  FATALITIES 


0000*^003000 


OOOOOOOOOOOOOOOOOOOOOOOOOOOOr* 


ooooooooooooooooooooooo 


I 


0 0 0003  330303 


0030003330030 


oooo^aooooor-»oooooooo*-«ooooooooooooooo 


003  3 0303333303 


0333  3 3 33  0 0 300033*H 


03*-t0000  3 003000  3000300000000030 


O 


300000000 


oocooooooooooooooooooooooooooooooooooooo 


ooo  ooooooooo 


o o o o o o 


00000003000000«-l0000000000000000000030000000000000r*00 


ooo*ho»ho«hooooo*hoooooooo«hoo 


0000*-I**000tH  3000000000000000000000  3 


0000000003303033300033300  3 O O O 3 rH  O 33 


330  33003000300  3 00330000* 


lOOOOOOOOOOOOOOOOOOOOOOOOf 


OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 


*Hc\ifOj-iDvDr^<£)(7'o»Hc\jfOj-ir)LDr^-cDCT'  orHcvjfo^-insXjr^co  o'  o h w to  j in  vor~co  <r<  o .-4  <\j 


in  NO  N O O'  Oi-tC\jK)^-lT)vDr-©(> 

i’j^^i-^inininmintninininin 
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000000000000000000000«-»0000000000«-»00000000000000»-lOOOOOOQOOOOOOOO 


ooooaooooooo  3 0030ooooooooor-ioooooooooooooo.-iooo*-*ooooooooooo.-toooi-»o 


OOOOOOOOOOOOOOOOOOOOOOOOOOOOOftOOO  30  3 00  3000  300000«-t00.-10.-«000»H00003 


000000C0000000000 


ooooooooooooooooooooooooooooooooooooooooooooo 


>£pw3)a\o^HOjro^invor'®cr'OHc\jro3-ir)\cr*coa'o*-<o<roa-mvX)r,*<Da'o«Hc\jK)^-insX)h»a)<ro»-ic\jroa'tn'flf^®<y'o,^f^»o^’ir»'i)r*-aD 

3 03  3 O'  (7 OOOOOOOOOOrtHHHWHHHrtrtOJW^(\(OJ(M(M(\)(\|(\jflrOK)(OK)K)KllOlOK)JJjjj  J- 

^^^rHrHr^^rHr^^^r^^^(\jcvooc\jc\jc\jr^c\jrvjrvjc\icgrgcvfMC\jc\jfV(^c\j<\jc\jCVC\irgcvcvjcNic\)(\jrvjf\jfVrgc\jojojfVfNjojCvjc\jcgc\j(Ni(\j(>j(>joj 


OOOO^OOOOOCOOOOOO«HOOOOOOOOOOOOOOOOO«HOOOOO 


000000000000003000001 


ooooooooooooooooooooooooo«hoooooooooooooooooo»h30ooo30oooooooooo« 


000«H00030030000000000000000r^OO«H0000000000  0 300000 


O O O O rH  O 3030  3 3 0 0* 


OOOOOOOOOOOOOOOOOOfHOOOOOOOOOOOOOOOOOOOOOTHOOOOOOOOOOOOOOOOOOOOO 


f0®’lTJ3^3CT'O*Hf\JP03-lDv0^-3(T'  O«H(\JK}3-in3r-C0<r'OfH<\jr0^invGr-®a'O»HC\lf0.3-  in\DMCff'OH(\i(OJiOyDN®J'Or<(V"''J  * 

c'J(^<^<\joj<\iOjK)fOrOK)K)K)K)rOK>fO^-^^^a-cj-^^’^^-tOtniOir)«r»inmiOir)m333333v03vCvDr^r^r,“r>‘r^r,'r-r*-r*-r^.o3-v'a''1x't 
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ooor-toooooooooooooo»-tOGOoo«-iooGO»-iooooooooooc«HOoooooooooooo 


OOO33O0O0O300OOOO3  0O0OO3OOO0aOO0OOOOOOO0OO0OO.-t0OO0OOOT-IO 


O'  O 


3 3 3 r-l  3 3033  303330  33333333330033  33330  33  3033303333303000030 


II  II 

V) 


Id 

o 

t/] 

o u 

U *-H 
<1  h- 

ooooooooooooo*-iocooo»hooooo«hooooooooooocooooooooooooooooo  *-H 

_J  _J 

<i  <x 


(NjrOj-irivflf^coa'  Ot-i 

rO(OK)rOK)iOfOK)(OK) 


c\jrocfiT>vDr~cDcr'  o th 

C\j<\JC\l<\jC\J-\J<\IC\irOfO 
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1972  DISTRICT  COUiRBIA 
NUMBER  OF  ACCIOENTS  WITH  ANT  INJURY 
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1972  DISTRICT  COLUMBIA 
ACCIDENT  TOTaLS  FOR  EACH  DAY  OF  THE  WEEK 


1: 

139 

264 

329 

126 

NO. DAYSr 

52 

2: 

138 

253 

371 

129 

NO. DAYSr 

52 

3 1 

168 

218 

325 

130 

N'O.DAYS  = 

52 

4 : 

155 

276 

408 

144 

NO .DAYSr 

52 

5: 

146 

318 

483 

329 

NO.DAYS= 

52 

6 l 

139 

321 

422 

328 

NO. DAYSr 

53 

7 : 

110 

222 

288 

132 

NO. DAYSr 

53 

1972  DISTRICT  COLUMBIA 

ACCIDENT 

TOTALS  FOR 

EACH  MONTH 

I! 

85 

170 

198 

111 

TOTr 

564 

NO. DAYSr 

31 

2: 

91 

129 

178 

87 

TOTr 

485 

NO. DAYSr 

29 

3: 

70 

133 

180 

96 

TOTr 

479 

NO. DAYSr 

31 

41 

64 

163 

217 

133 

TOTr 

577 

NO. DAYSr 

30 

5: 

80 

172 

256 

112 

TOTr 

620 

NO. DAYSr 

31 

6: 

73 

150 

242 

115 

TOTr 

580 

NO. DAYSr 

30 

7: 

78 

151 

230 

109 

TOTr 

568 

NO. DAYSr 

31 

a: 

78 

142 

228 

93 

TOTr 

541 

NO. DAYSr 

31 

9: 

85 

189 

228 

121 

TOTr 

623 

NO. DAYSr 

30 

10: 

103 

177 

229 

119 

TOTr 

628 

NO. DAYSr 

31 

11: 

106 

154 

236 

103 

TOT  = 

599 

NO.DAYS= 

30 

12: 

82 

142 

204 

119 

TOT  = 

547 

NO. DA YS= 

31 

1972  DISTRICT  COLUMBIA 

AVERAGE 

ACCIDENTS 

FOR  EACH  DAY 

OF  THE 

WEFK 

1 : 

2.7 

5.1 

6.3 

2.4 

AV  = 

16.5 

2 t 

2.7 

4.9 

7.1 

2.5 

A V = 

17.1 

3 • 

3.2 

4.2 

6.3 

2.5 

A V = 

lb. 2 

4 : 

3.0 

5.3 

7.8 

2.8 

Ayr 

18.9 

5 : 

2.8 

6.1 

9.3 

6.3 

A V = 

24 . 5 

6 I 

2.6 

6.1 

8.0 

6.2 

A V = 

22.8 

7 ! 

2 . 1 

4.2 

5.4 

2.5 

A V = 

14.2 

average 

1: 

1972  OISTRICT  COLUMBIA 
ACCIDENTS  FOR  EACH  MONTH 

2.7  5.5  6.4 

3.6 

AVr 

18.2 

2: 

3.1 

4.4 

6.1 

3.0 

AVr 

16.7 

3: 

2.3 

4.3 

5.8 

3.1 

AVr 

15.5 

4: 

2.1 

5.4 

7.2 

4.4 

AVr 

19.2 

5: 

2.6 

5.5 

8.3 

3.6 

AVr 

20.0 

61- 

2,4 

5.0 

s„i 

3.8 

AVr 

19.3 

7: 

2.5 

4.9 

7.4 

3.5 

AVr 

18.3 

e: 

2.5 

4.6 

7.4 

3.0 

AVr 

17.5 

9: 

2.8 

6.3 

7.6 

4.0 

AVr 

20.8 

10: 

3.3 

5.7 

7.4 

3.8 

AVr 

20.3 

11: 

3.5 

5.1 

7.9 

3.4 

AVr 

20.0 

-121  - 

2.6 

4.6 

6.6 

3.8 

AVr 

17,6 
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1972  CONNECTICUT 

NUFBER  OF  ACC  I OENT S WITH  ANY  INJURY 
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1972  CONNECTICUT 


ACCIDENT 

TOTALS  FOK 

EACH  DAY 

UE  THE 

WEEK 

i : 

345 

978 

1025 

363 

NO.DAYS= 

5? 

2 i 

392 

943 

1 072 

417 

NO.DAYS= 

52 

4 : 

433 

1016 

1152 

477 

NO.DAYS= 

5? 

4 : 

455 

1004 

1242 

619 

NO.DAYS= 

52 

s: 

4 u4 

1169 

1727 

1263 

r 0 .DAYS= 

52 

6: 

314 

1275 

1528 

1199 

kc.days= 

53 

7: 

258 

894 

1142 

405 

no.cays= 

53 

1972  COT  NECT ICUT 
ACCIUENT  TOTALS  E OK  LACK  FONTH 


i : 

270 

518 

564 

318 

TOT  = 

1690 

NO  .DAYS  = 

31 

2 : 

20? 

528 

599 

311 

TOT  = 

1647 

NO.DAYS= 

29 

3: 

18  2 

472 

605 

343 

70T  = 

1602 

NO.DAYS= 

31 

4 : 

180 

622 

678 

313 

TOT  = 

1793 

NO .UAYS= 

30 

5 l 

216 

627 

833 

367 

TOT  = 

2043 

NO . 0 A YS  = 

31 

6 I 

203 

66  0 

871 

474 

TOT  = 

2208 

NO • D A YS  = 

30 

7 : 

189 

664 

752 

421 

TOT  = 

2026 

NC.DA YS= 

31 

a : 

179 

620 

735 

426 

TOT  = 

1960 

MO.DAYS= 

31 

9 : 

221 

626 

809 

435 

TO  T = 

2091 

NC.DA YS= 

30 

10: 

271 

631 

778 

458 

T G T- 

2138 

NC.DA YS= 

31 

n : 

262 

581 

736 

377 

TGT  = 

1956 

NO.DAYS= 

30 

12: 

219 

730 

908 

500 

TOT  = 

2357 

NO.DAYS= 

31 

1972  CONNECTICUT 

AVERAGE 

ACCIDENTS 

FOR  EACH 

DAY  OF  THE 

WFFK 

1 : 

6.6 

18.8 

19.7 

7.0 

A V = 

52.1 

2: 

7.5 

18.1 

20.6 

8.0 

AV  = 

54.3 

3: 

8.3 

19.5 

22.2 

9.2 

A V = 

59.2 

4 : 

e.e 

19.3 

23.9 

11.9 

A V = 

63.8 

5: 

7.8 

22.5 

33.2 

24.3 

AV  = 

87.8 

b l 

5.9 

24.1 

28.8 

22.6 

A V = 

81.4 

7 : 

9.9 

16.9 

21.5 

7.6 

A V = 

50.9 

1972  CONNECTICUT 

AVERAGE 

ACCIOENTS 

FOR  EACH 

MONTH 

1 : 

8.7 

16.7 

18.8 

10.3 

AV  = 

54.5 

2: 

7.2 

18.2 

20.7 

10.7 

AV  = 

56.8 

3: 

5.9 

15.2 

19.5 

11.1 

AV  = 

51.7 

4 : 

6.0 

20.7 

22.6 

10.4 

A V = 

59.8 

5: 

7.0 

20.2 

26.9 

1 1.8 

AV  = 

65.9 

6 : 

6.8 

22.0 

29.0 

15.8 

AV  = 

73.6 

7 : 

6.1 

21.4 

24.3 

13.6 

A V = 

65.4 

a: 

5.8 

20.0 

23.7 

13.7 

A V = 

63.2 

9: 

7.4 

20.9 

27.0 

1 4.5 

Av  = 

69.7 

10: 

8.7 

20.4 

25.1 

1 4 .8 

AV  = 

69.0 

11: 

8.7 

19.4 

24.5 

12.6 

A V = 

65.2 

12: 

7.1 

23.5 

29.3 

16.1 

A V = 

76.0 
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MORNING  EVENING  NIGHT  + MIDDAY  --EQUIV  MORNING--  --FQUIV  EVENING--  NIGHT  + MIDDAY 

UAM-10AM)  (RPM-IOPM)  MIDDAY  +NIGHT  DAY  DATE  ACCDNTS  DAY  DATE  ACCDNTS  MIDDAY  *N I GHT 


AAAAAAAAAAAA/NAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 


ov£v£*-*(\ij-o{\3-v£r.r~iDfOfM»a;o»oir)or^o^ovDPOcD(t)r>-o(7'vD(rcc»noK)(r'ocr'©oor-(occw^iv£vcircao4 

CgCyCVKJCwd-sOCMCVjrHr^CMd-kDrH  H(VlOK)U1Cg^Hlf)HJinrMOVrMir-rU3-tflH  CM  CM  v£  K)  CM  PC  lT>  PO 


oooooooooooooooooooooooooooooooooooooooooqoqooooogooo 


vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvv 


omK,)io*Hcg^N*-»(rvflvoa-(\jfViinof-»aNO(\iocu^««Hcoina'Of-ij-ir)©c\jr>-»o*-io^-^-oor-ooK)*Ho©^ooK) 

C\,H(\i3-K)C\J(MOJWlOfOfO<\iC\ir*K50J(\jK)(\l  *-i  & CM  CM  PC  CM  CM  CM 


o^r-®(7'<o«HC\K)^sflN£h-©Ki^-o(T'Or-i(\iNLr^ffiO'K)3-invDOricvjror‘©<roctfnN£ooc\.oo(ro*-i^NCc;cj' 

*-HrH*-lCMCMCMCMCMCMCMCMCMCMCM  •-<  CM  PO  PO  CM  CM 


afs'rHC\K3ctir)v£rs-»H(\K)j-in'£r^oc\jroj-ir)\i]r^*-ir\jKijir)vflr-rHr\,fO^-ir.  \£r~cc\K)aa»£ooovK)j-ir.  >£oarg 


ooooooooooooooooooooooooooooooooooooooooooooooooooooo 


000000003 

ooooooooo 


ooooooooo 

ooooooooo 


oo  oooooo 

\\\\\\\\ 

oooooooo 


ooooooooo 

\\\\\\\\\ 

ooooooooo 


ooooooooo 

VVVN.X'S.V.V'x 

ooooooooo 


ooooooooo 

ooooooooo 


00000000000000000000000000000000000000000000000300030 


AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 


oo'oo3'©«-tror- 

:mpocm\DCM.*-po.-i 


ajajojojtorDoo 

H rt  (M  (\|  IO  C\J  rt  Oo 


o t>J  H o KI  Ml  O \D 


*-*r-ina'(\j®f\jir>o 

C\lJinrM(\jHK)fOIO 


MDO^^OCNO 

-o  vD  *— » •— * cm 


00l0\£0"000^ 

(\  -H  (\  J IO  CM 


OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 


vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvv 


p^oJ'nDj-  oiPro^iPj'sr-jforYij'nin  -oxcr'o  dt'lti  m d j-  Mr-*-i:oo.or^r'»  ? h jjbh  m ,0  -m  *o  r 

«— I r— 1 rt(\J(\/(M(NlCMHHHH(MCVHHHHfOC\jH(\J  «-*•—«•-*  O'  C\J 
| | I I I I I I I 


v£r-(T'K^c\jc\ia-ju)®ir)(7'vflNccjj-vDrvinr^®®(j'j-j'Dvi)vDvOir-rvirioo'oina30'?vDvi;ccr^irj'3  j r-  - oo 

III  l(\|<-t*HI||IIOIIIIIIII*HIIIIII<HIIIII(\JI  I I I I I I • I • • I I • I • ^ 

till  I II  II 


H(\jfOJiT)vDr-cDiT'OHt\jK)j-inv£r^coa'OH(\rojinaj(^fflC'o^ 


h cm  mi  ? in 

\ \ \ \ \ 

CM  CM  CM  CM  CM 


vO  r»  <o  ex'  o 


cm  ro 


\\\\\\\\ 

CMCMCMCMCMCMCMCM 


kPsDr'®a'o^-«x 

\\\\\\\\\ 
CM  CM  CM  CM  CM  CM  CM  CM  i\ 


& O' 


CM  fO  3-  lT>  vD 


CM  rO  d"  vO  f*** 


CM  PO  cJ-  iT»  vD 


CM 


J IT)  i)  S 


cm  j in  j) 


(\|  Ml  J iT  v3  P' 


CM 


at 


M 


119 


/2  3 


AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

oinMrvflinooH(MncK)jnMO(M^inN(\jvflHO(DvC\i)(\iHicio(\ivo 

(\jcvifOjHc\jiO(\jfOd-inH(Vj(vroicio^(VHK)K:(\j^irw(\j(NifOK53-irfOvfl 

OOOOOOOOOOOOOOOOOGGOOOOOGOOOOOOOOO 

vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvv 

j-in^)^Ha'ir)^-j-ir)j-»Hv£iaoiri(jNCT'inc\j^)j-Lr)K)rocrc\j3-orsjfo^c\jr-r-ir) 

CMCM<MCMK)CMCM(\.HCMCMCM(M*HCM«-«CMiM<v<MCMCMfO^CMCMCMCMCM<Min;j-«Hro 


i 

t 


aaaaaaaaaaaaaaaaaaaaaaaaaaaaa; 

O (\1  O a J ^OJv£(V®a'tMriO(\iCa(MJ'«J)(rOOOOOOOOO 

po  k)  a-cMPOCMiocMinincMioroa-.a-vo^cM 

I 


OOOOOOOGOOOOO3-r^^O0NCDO<\JOin*-tf*“OvX3r^J- 

K)<\jinro*-«voa-»ooic\i3-^ir)^ro<M 


vvvvvvvvvvvvvvvvvvvvvvvvvvvvvj 

t 

( 

I 

i 

O'lOOOiniOJliniOHfflvflff'HOJff'Nrtl/IrtOOOOOOOOO 
cm  cm  cmcmcmcmcmiocmcmcmcmcmcmioidcmcm 


K?  ^lTvCOn(MK)hCCff'Or( 
C\J  CM(MCMCMCMCMCM«HtH«HK) 

00  CO<Oa3CDCO®CDCO<X'<O^J-m 


C\jK)ctinv£r‘C00'^-*H<\K1 
OJ  H n H 

irminioioioioioioioioio 


^ir^jeff'oofKitfoo 

coinininmininr^h-iDinoo 


CM  PC  IT  vC  ® O'  O H O O it 
C\J  CM  CM  H H H (\J  (M 

\\\\\\\\\\\\ 
^^t£\filOvOv£vfiv£h‘V£vfl 


CMPOcftTOOOOOOOOO 
CM  CM  CM  CM 

wssswwww 

vO  vO  \0  n£>  O O O O O O O O O 


K,).3-mvor^t-i(MPO3-ir)v£f'>r-icMK5.d-ir>Xf'-pHCMK)3-irj\0r^fHCMfO;fir)vor-*-tcMiiooovcr,»«HCMiO3'mvor,-»HCMfO^'invj3r-ooQoooooo 


OOOOOOOOOOOOOOOOOOOOOOOOOO 


o o o o o o 


vJ  U U 


ooooooooooooooooo3ooooooooooGoooooooooooooooooof0^-in'i)a'P0r-^'inK)CJ>®^-ir)^'O 

H CM  H CM  CM  H H f-l«H  1—i 

03000000000000000000003000000000000000000000000«H«H»H«H*HCMCM'»“1»-«CMOCMOOCMOj 


:OOOOOOOOOOOOOOOOOOOOOOOOOOOOOGOOOOOOOGOOOOOfHCMr0^-ir)s0h*»-tCMf03-insX)f^»HCM, 


A AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA  AAAAAA  AAAAAAAAAAAAAA 


r-a'cDCDK)cr'cr>vor'TH<xia)CM3-^-ir)cr'Oco^-cM*-t»-io3-ir)r*)cofOr^m«-«r^^O'CMoo^'roinr^cMCM3-(7'ro(T'(7'ooir)r“^-ooooooooo 

HlO(OJ'^C\,!MCM(MinjH(M<MfMHK)jHriJCMCMinj<-iHH(M(MlOK)H»HrtCM  ^(MCMvCCMCMIOC\jCMHC\jCMK)rOfOCM 


• GOOOOOOOOOOOOOOO 


o o o o o o o 


CMCMCMCMCMK>3-CMHf03-f0 


vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvv 


h i*J  D X)  D J''  ■— < 'r) 
(M*lC\lN)rtrt(MC\) 


iTJ'  ^ H O B D O 
(M  JK)H(MHHh 


I .0  iD  3 il  N H 
H CM  H H (O  H CM 


3-  DJ'^O.DOd-r^ 
(MC\,HrtHHHH 


vfl\3-0'j-DHin 


h-j'X)^-cT'X)cr'Ln 

HHhCVJHCMHH 


ON  J 1/10'  D vfl  K) 
CMHH(\|K)HCM(M 


pO  H if)  H •— I (JN  r-l 

CMCMCMJK5CMH(\J 

I I I I I I I I 


rHCT'r^o'j-<MCMH^ovor^r^3-»-tf^>OvDiiocrvccMr^cD^rr«3-r^ioiDcrvoocor^irmcM*^j<Mvi3J'»^ooir)rs*oir>ir)CT'CMco^-ir)(MfOvoco^‘crr^ 

I I I I I I I I I I I *H  r-l  I H r-i 

I I I I I II  II 


"vjp'jj-invur^cDcro^icMPO^-mvor^cccror-tcMro^-Lnvijr'ajcr'OHrHCMKicriD'Xip^cocjNo 

CM  HHHrlHHHHHHCMCM(M(MCM(M(M(M(M(MlOlO 


J in  vO  N D 

CM  CM  CM  CM  CM 

\ \ \ \ \ 

CM  CM  CM  CM  CM  CM 


CMfOJlTvflNCOa'OHCMrOJin 

CM  HHHrlHHHHHHCMtMCMCMCMCMCMCMCMCVtOrO  HHHHHHr.HHHAJ(M<MCM(\)(\l 

cMro^ir)'jr^«H\j'Od-j3vijr*^«HCM'O^ru3v£r**-rHCM'O.^’Lns0r^«HCM,o^-ir)\X)f*~t-i(Mro^-ir)vi)r^fHCMfOcrir)v0fs*»HCMro^-inv0fs'»HCM 


120 


AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA/NAAAAAAA 

OOOOOOOCOOOOOOOOOOOOOOOOOOOOOOaOOOOOOOOOOOOOOOOGOOOCiOOOOOOOOCOG 


p^osi3oaDinr-f^mccor-iric\jacccr'«Hco(T'cDir)Hd-focDfHoa'^-(T'in»ooHvX)^-v£Ln3-ars~vD<%jfO®c\jcj-r-oojfO(ocg3-^-vi}c\(»oiD^-o^ 

CMCMt0O;fCM*-<CMfHinincMiHCMr0CMCM's03-CMr-irocMCM3’t0CMCMCMCMiM:3-3-  rocMCM»omfOCM«HCMro.3-^-®CM»HCMtO;*-^-(D<M»-i<Mto^-tntocM*c 


vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvv 


ooaoooooi 


lOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO* 


oooooooooooooooooooooo o 


aooccoocaoooaoaoaoooooooaoacoooooooooooooooGOcoocooooooGoocoGGG 

ooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo 

ooacooooaoooooaoooooooooaoooooooooooooooooooooooooooooooooooooo 


to  incMcr'^vDr»r*»'<x)tnin®a'into«-if^r*'CM3-3-in®o 

H H H H H W H H H 


LT  J I'd  I 


CM  t O CM  CM  (M  CM  rH  <\ 


v0in®CM*-«CM®(7'®T'in® 
H CM 

\\\\S\\S\N\\ 


r-r-»c\jTHCM^3crooin^or^ro 

CM  H H «-t 

CMT'C7'POPO(r>(7'(T'acr'C7'Cr'rO 


rO  to  K)  cm  ro  rO 

\\\\\\\  \ \ \ \ \ 
lO'OOOO'CDO'O'OlOOCD 


tO  CM  *0  »0  -O 

O O'  O'  O K)  O (MO  K1  >0  O 


10  ^-iTJvDr-r-iCMfOJ-invXir^ 


<MK)j-jnvX5h-f-«cMf0^ir)vi3r^»HCMio^'ir)v£ir'OCMro^-in43r«»HCMfio^-ir.  vor^«HCMro^inv£r*“^-»CMfO^-»rvijr^«-HCM 


AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 


ooaoaoooaaooaoaoooaoooooooooooooooooooooooooooooooooooooooaooo- 


•H^-co«-tocM^’0<r'\DvOsDavDCMcocJ'3-cDr*-r'fONCororoin»HCMr^K)o»-»oir)^DK)»Hvom3-ioa-for>sD*-<^-jfOj-r03-u305->oir>®®^^'^r 
H(M(M'OJOM(\JiO,OCMiOK)(^iOCMHiOJjCM(\j(M<MlOvfliOH(MrOfO(\liniO  K)J'0r0JKlr0C\j,\|f0(M'DinCMCM'0''MK)inj^K)K)j'j.'r''.i.\ 


VVVVVVVVVVVM«'VVVVVV\/V\/V\/VVVVVVVVVVVVVV\/VVVVVVVVVVVV\/VVVVVVVVVVVV 


'D^fvoosinj'T'inAi 

HHCMCMCVJCMHCMHCMCM 

I I I I I I I I I I I 


sD-J-lDlOrO^-  "M  CM  O CM 
CMCMCMCVJfOJ-CMCMCMCM 

I I I I I I I I I I 


cm*h-mp^ct ' n o'  "O  r>  & 

CM  CM  ID  CM  H CM  CM  CM  fO 

i i i i i i i i i i 


T*  r*~  o M1  b h - 1 ■o  vo  in 

(M(M(MHrOHJK)J(MCM 

I I I I I I I I I I I 


OHiOOU'OtlHO  3 
CMrOCMCMCMCMfHCMfOvfl 

I I I I I I I I I I 


CM-3CM<MCMrOiPCMC\CM  — 

I I I I I I I I I I I 


(7'irc7'm^'CM^’Croa'Or-ivi)r-ir^N£OvDincMo^-®incMvflJ’oroin(7'invot|oroino®crcr\Da'vD^(7'0'vDoina-a'r^p**»or^ 

H H H I H H 


\Dr*cD(j'OH(Mio  3-  in  ® ® cr  a^-iCMmcrin®r'®cr'0^cMPOj-invor-coa'OHT-icMrojin®r- 

CM  CM  CM  CM  fO  HHHHHHHHHHCMCM<M(M(\J(\l(M(MCMCMfOlO 

4-3'3-d>3-if)ininininininininininininininif)inininininininininininininininiOviiii)vDvovOvi) 


\ \ \ 
\Q  \0  vO 


nCMfO  j 


\ \ \ \ 
vX)  vD  v£)  si) 


in  vd 


r*-  ® O'  o *■«  cm  »o  y »'  >4?  p*» 

H H fd  CM  CM  i\  |\  >\  >\ 

xxxv'xv.'svxx  X 

vO  vO  ^ vfl  ^ O s*.'  \4 


to  3-  iT»  ® r- 


cm  to  3-  in  ® r*» 


cm  to  3-  in  ® r> 


CM  to  3-  in  ® |X. 


cm  ro  3-  in  ® rx 


CM  IO  3-  in  \D  P*> 


cm  to  3-  in  sfl  r- 


cm  •o  j in  \0 


CM  *>*■«-»  4) 


121 


ae/ 


I 


AAAAAAAA 


AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA, 


v£^'DlTlOHOCOd-0'in 

cm  cm  to  in  3-  cm  cm  cm  cm  3- 


vvvvvvvvvvv 


rooovor^j-coiHCjNco 

3-cmcmit)cmcm3-3-cm«h 


vvvvvvvvvv 


otoo'(Mcoinf\jinHvDa' 

CMCMCMlOCMr-lCMCMr-lCMCC 


vvvvvvvvvvv 


(J'®C0f'-rO3-(7'CCrOr^ 


VVVVVVVVVV 


^oK>cor-r-o^ojh-ff> 

(MCMromrorocMCMCMr'-.s- 


vvvvvvvvvvv 


rto^mmcDcvirtio^, 

3‘CMincMroi'~rorocMc-(, 

! 

l 

vvvvvvvvvv' 


OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOQOOOOOOO<MK)C0i 

Od  CM  *H  ; 

1 


cr  oooooooooo 
ooooooooooo 


cooocooooo 

oooooooooo 


c ooccoocco 

\\\\\\\\\\ 

oooooooooo 


ooooooooooo 

\\\\\\\\\\\ 

ooooooooooo 


oooooooooo 

\\\\\\\\\\ 

oooooooooo 


OOOOOOOOCMvTh* 
CM  CM 

OOOOOOOOf'-vflvD 


OOOOOGOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOr'**H<\J 


r-\DOvDvO(\jo^<Hd-K5\£ir)K)a]rr)or^co3-j-Hor^r'ir)C7'iDt*'j-r'Cof^inc7'r-jK)Hr^ffl\£itoind]vfir^irir^r^ffiod'i/iin3‘ir'cEiin(03‘3- 


>Dr«-ro^-or-«orHnjco(T'Ovr'V)f03-Lnvor-<D(7''0»Hc\jrooinvDr*-h*a30r-id-inA/in3‘in«of^®a'o^»<\jfO^-in>Dh‘^)(7'o»-((\(fo^-if)c\j<\jr>‘vX), 
«-<  i-H  ft  HHHHHHHH(\|(\i(M(\|C\l(M(MWHHfO  H H f\J  CM  CM  Od  CM  CM  CM  rO 

CT'(T'»OfO(7\<T'0»H:0(r*(7'Cr*rHCr>(T'(X'(7'Cr‘J,*C7'(7'«CM(r»<7'(r>CM(r*0'*(J'CMCMCT'OrMCM«-<0«-t«Hr-l«Hr-l»H«HCMCMC\J<MCMCMCMC\jCMCMCM<MCMCMCM«HO»H  cm) 


vGo-.-icMr03-insors-f-«cMr03'inv£r^«HCMro.3-inN£r-r-icMro3inM)r-r-tCM 


i 


aaaaaaaaaaaaaaaaaaaaaaa 


A A A A A 


OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOa'rs'C\l 

c\i  cm  ro 


c\jrofoir'cooof\jfoo^r~»Hc\jc\jrooj-®crvD»HCT'ocrr^h'Oir)0,'«Hvcr'®j-o^Ho©(t)inoc\jr^cor-co\£)r*tn®»H3'd-a'3-®ror-j-vflr*-o(j' 
3"  K)  3"  3"  3"  3-  vDrO(M'OjrOC\jHd)fOfOfO(\J,OH^K1fOlOrO(\)C\)(\J£\jPOrOd)J-K)'OlOrOKlif)'OOJ>OfOK)K)(MK)POfOK)iniO(\J(\j(MH'OrOJ'OC\J 


VVVVVVVVVVVVVVVVVVVVVVVVVVVVVVVVVVVVVVVVVVVVVVVVVVVVVVVVVVVVVVV 


T>r»^-rr'vj'Ov07'oj- 

C\J(\jC\jC\jC\jC\j(\jrHf\jC\j 

I I I I I I I I I I 


o r-H  — < O'OT'j-ojcoin 
lO(\l(\|(\|K}rt(\J(\;HaH 

I I I I I I I I I I I 


n-o^in-ooinT'j'o 

I I I I I I I I I I 


»o  n in  o'  -o  don  j o in 

I I I I I I I I I I I 


ojiniH-or-N/ij 

C\JCVJPOC\.C\J'NJC\)tH»HOJ 

I I I I I I I I I I 


h in  ^-3-iT)sO*H(7'r^,o 

(\JK)C\lC\J(\jC\j(\JC\jHHfM 

I I I I I I I I 


3-  invcr~v£CMcMf'~J'vOr03-f~ococoroco3-cocoo3-3-ror*“inf*-a3*-<3-3-cr'rororor<'r^vUvijO\£)  vcoi 

H I H H H H H 


l(\l3ini!'3<MvC3f^ 


O'  o h (\|  ro  3 
cm  ro 

\ \ \ \ \ \ 

mj  vd  r-  r>- 


in  dj  ® O'  o 

\ \ \ \ \ \ 
r-  r^-  r-  r- 


iH  Od 


ro  J in  vO  N ® O' 


\ \ \ \ \ \ \ 
r*-  r»  r*-  r-  r-  r*- 


o H (M  ro 
CM  CM  CM  cm 
\ \ \ \ 
r-  r-  r- 


3-  in  vD  r- 
CM  CM  Od  CM 
\ \ \ \ 
r-  i- 


<x>  cr>  o *-• 
c\j  cm  ro  ro 
\ \ \ \ 
r-  r — r-  r- 


rH  CM  ro  3- 

\ \ \ \ 
CD  CO  CO  CD 


in  \o  n ® 

\ \ \ \ 
® CO  CO  CO 


a'OTHCM'03-ir'vor^-®CT'OHCMro3-invor-®cr 
hhhhhhhhhh(MC\(M(MCMCMCM(MCMCM 
\ \\\\\\\\\ \\\\\\\\\\\ 
©Q3COCOCO<X)COCO(OCO®COCOGOCOCOCOCO®COCO 


ro  3 in  d)  N 


CM  ro  3"  in  >0  P'- 


oj  fo  3-  m ® r- 


C\l  rO  3 in  \fl  I** 


CM  rO  3 in 


vC  r*-  H CM  ro 


m vD  n 


CM  ro  3 n vll  N 


(M  rO  3 in  d)  S 


CM  ro  3 in  \0  f- 


122 


AAAAAAAAAAAAAA/\AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA/\A/\AAA 

w<i»flNoo(j'NcooinRicyiNH(r'Hoininhff'MMnoirooinMno>f*(M»H>fivflHoo(Mi)inoioff'^^ttia)NoooioiroNNooc 
3 to  *)  to  <xir)*^a-ropo»-*io<v<\i3'ir<v»Hio>o»03-in«-tc\iK)oj<\)KJd-cg<vK)cs(K)invi)  cMCMCM.3-in.-icM»Hin*Hcrin  nru  ro  10  in 


®o^cgoor*ooooooooooooooooooooooooooooooooooooooooooooooooooGocooG 

eg  « ^ ® eg  eg 


vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvv 


CPo^fOooina‘iPOj'OP“in(PCMiOd‘CMvOa-inr^<p^«-4®tn;i-3-in®«H<pp*>o«Ho®i‘^ooij-:j-cMP'ro»-»^-in®incpooo«-icpocMCMooo 

C\  CM  CM  CM  HIOe*IO(V«HHCM(\|J(VOj(Ni(\rtrlO<(V<NJ(V<>j(V<M(NJC\|Hrt(\J(\|llO(\J3-  riH(\|C\M0(\j(\|HK)K)l0»0  C\  H (\  J K1 


« 

vC 

3 

m 

a 

o 

® 

r- 

® 

CP 

o 

CM 

r- 

K) 

CM 

to 

® 

CP 

vC 

3 

® 

CP 

O 

oH 

CM 

90 

3 

in 

® 

O' 

o 

«H 

CM 

o 

o 

in 

v£ 

® 

CP 

o 

CM 

m 

3 

in 

c 

o 

CP 

o 

CM 

to 

O 

o 

c 

CM 

•H 

to 

CM 

CM 

CM 

to 

CM 

oH 

IO 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

fO 

to 

1— 1 

*H 

«— i 

rH 

t-i 

t— i 

«H 

CM 

CM 

CM 

CM 

V 

X 

N 

X 

V 

X 

X 

X 

x 

X 

\ 

X 

% 

x 

X 

X 

X 

X 

\ 

\ 

x 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

® 

r- 

o 

o 

in 

r- 

P~ 

p» 

r- 

in 

® 

m 

® 

in 

in 

m 

m 

in 

® 

in 

® 

® 

3- 

® 

® 

® 

CO 

® 

® 

® 

CO 

® 

® 

® 

CP  <p 

o 

o cp  (P 

CP 

CP 

CP 

CP 

CP 

(P 

CP 

CP 

CP 

CP  O 

o 

CP 

CP 

CP 

CP 

CP 

o 

o 

o 

to 

3 

m 

vfi 

o 

o 

CM 

KJ 

m 

r~ 

- 

CM 

rO 

in 

® 

p* 

- 

CM 

to 

3 

in 

vC 

H 

CM 

to 

CP 

in 

1^ 

- 

CM 

d- 

in 

vC 

o 

o 

CM 

to 

m 

p- 

r-i 

CM 

K) 

in 

o 

o 

CM 

to 

3 

m 

v£ 

o 

o 

G 

CMCPP^^-OIOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 


r^vOP-Gf-io^oooooooaoooooaoooooaoooooooooooooooooooooooooooooooooooG 
cm  eg  cm  fO 

CMrHOOCMOtHOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 
H H H 

l03‘in'£P»OCMOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 


AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 


^%ovo»Hoo(n®»noioa'vD(7'®®r**ovovoroir)ai,*»r-®or-r'3-cr'\i)'X)CMinoo^-imoo^-®r»»-*vo»ocMOvooifJCMOoo\o*-*p-»ooo-> 

eg  in  ® CMc\i(\i3,inHC\jHiriHjin(\i(\iHOjio»oinH(\jH(\ic\iNN(\iaj(\i(Mio\flffl  <\i  c\j  »o  ,p  s 


OCM'l0®OO*-»OOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOOO 

J J 03  rO 


V VN/VVVVVVVVVVVVVVVVVVVVVVVVVVVVVVVVVV  vvvvvvvvvvvvvvv  vvvvvvvvvvv 


Hoos(\jHj-^-ojr*-,OHj-irit)in^or\jHHU'o'\i(Min^  o\ir'-cpo:j--*o:\jor-HCPr'~:r’,o  ohsoiptoo  o 

(NJfOCM^CVJtNjHHOJtMinCNJfVJHiOfOKJfOMfVlCNJHCgd-fOCMHHCXKMin^CXKMWCMHiniOWCMfNjCNJCMJCMW  *-«<MCMu''3-r-«*-i*-i»M*\»'\.-’  — — — 
II  II  I I I I I 


<P  ® ® \fl  m 


inp‘®;*-»o®CMio®fOor^®:rci-p-or,-r-»CMCPinr-vfiin®f^vOP-®^cM<PfOv£>®t^vi)®r*-<P3-roo©®in.»iocp.x'  <p  p <x  v 

rH  I I I I I *H  I I «H  «“ I III  iH  | f— 4 I I I I I r-i  I lllr-lrHlll«-lll«HIIIIII«Hr^<Mlllllllll  1 1 

III  II  I II  II  III 


OrtHOjKij-m'Dr'Cocr'OHCMfOj^vfliv-tDir'OHAjro3-ino)rs-(DO'OHC\iiojin^Divcotr'OH(\i'Ojinvi)r^a?iT'o*-(V'1'  s • 

Ki  m H H H H H H H H H CM  (M  (M  (\  «\  .V  .\  >\  A \ 

CC  cD^CMJ'ff'flNJ'gN|MJ'ff'(MM'<MM^lMMMMM3'CMMMMMMJ'^OOOOOOOOOOOOOOOOOOOC  OO  o o o O v'  - * ' - 


i^-in«i)h'r-icMto^-ir)®rs~»H'AJfOJ-  in  ® < 


ioj',OJinvDr'-rHC\jr'0:t‘ir>®r-~'-*e\j>o.yj'>®r^— «>\io  s x "•  — \ 


123 


AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 


r-if)(0^-cooor*^ovOvDvc*H<\j<Dinr*‘^-<\j(j'coininroa'(T'o<\iH(rroa'a'c\j<\ja-a'»iO(7'(T'cyjcsjHmK)^ir)oo»Hr*-oo»rooinr»-o 
h oy  j in  h cvcjirr*‘C\jc\j<vro^-K)»oojcv<\j  P'Ococ\j.HC\jc\jc\jvOCMOj*H(\jcs.c\jioc\jCMcHCMCMOjv0ccc\j«-rojojK)vDio  oj  <m  po  po  k) 


ooooooooooooooooooooooooooocooooooooooooooooooooooooooooooooo 


vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvv 


vDd-c\jcoinooin(\j»HoinK)ror-iina'vD'i3or»>oo(r'f^ir)oh*^-ff'cva'ir)oh'^-o'cou'ioor^^-ONc\j^inorH3'vDcvivfloo^HK)if)a'v£)o 

h h aj  r\j  (\  cvc\j-^-(\jc\c^c\jooro^-if\jc\j»-i  m<)rt(\iHHrtinHH(\iH(\jH(Mrtn(\jH(\irtinrOrt(v(MnrtirH  cm  cm  c \i  to  h 


\£MDff'JOOMDr<OHhCD0'OJif;v£(\|HCO^MrOH(\Ulvl)(C(rOrtir;>£N(80'OHir\£yflMJ'Ol0(VilfUflinoolflJI0Jlf1O 
C\J  OJ  CM  CM  C\J  (\J(M(\)IO  HHHWHHH  H «H  «-»  cH  r-l  »H  »H  H H <rH  *H  tH  H H H 

\ \ \ \ \ \ \ \ \ 


roj-in'Xip^ooK)3-ir)v£r^*-ic\jro^-ir)sCfs‘r-i<\jooir)s£ir^r-iojf03-if)vDr' 


c\ifOcrinvcr><^ovjr03-irivxir^*HC\jro^-invX)ooc\jK^^-inv£io 


oooooooo  oooooooooooooooooooooot 


OOOOOOOOOOOOOOOOOOOOOOOOl 


o o o o o 


oooooooooo 

\\N\\\\  \ \ \ 
0030000000 


3000030300 

\\\\\\\\\\ 

oooooooooo 


3000000000 

\\\\\\\\\\ 

oooooooooo 


ooooooooooo 

\\\\\\\\\\\ 

000300000  oo 


oooooooooo 

\\\\\\\\\\ 

oooooooooo 


oooooooooo 

\\\\\\\\\\ 

oooooooooo 


o OOOOOOO  00000000000000000000000003003000000000000000000000000 


AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 


«-»r-t®(s'3-OOCDHia' 

CM  O-  ^ CM  fO  (M  ^ 


ir>a'a'vjjHfnro(\jojH 


CDOOOHO(MvJ3(MvOvO 

vJ3d-POOO»HfOC\j3- 


®ocr'infO(T‘o®ro® 

j-fOKiHCMHOjiMH 


inotMiBNOoajNy 


IO(\|OOCM(M©(OvflO 
CD  ID  CM  CM  CM  CM  vD 


0000000000000003000000000000000000000000000000000000000000000 


V vvvvvvvvvv  vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvv 


■c  oo-HOio^-N-gN-^^r-cr^-'O  \inH\jg'noHinjjro?  0'Hir>(s~'\j-omfs>r«:}-Ln:r'.n:\jo,'Moor'-oj7'inrs>rs-j':3-Ajinvi)«-*ir>:om-o 

r-i»03'*HrHcr<\,c\)rocMC\j»HrMro»H^HrHCM«-iCMCMforocMrHaj»HC\jcMroc\j«-*<\jfOC\jCM^-»o«Hr^^-icMr^  vOK)(\jHCM(\jHfyjinn 

II  It  l«l 


iPCT'f*)cr3-iDoir)cDC7'C'Osjr^oro^oir)ccr-tnc-inror-ir-ccvnrOLn'XvDc-Lr!r-p-cD03-iniDj-incoinc\rH(rrooj-roinir)^-crr-o>K)jcj- 

l«Hr-ll  I I H H I I I I I H H(M  H I I «— i • — 1 0 I I I I I I H I I I I I I I H H I I I I I I I I I I I H OJ  H I I I I — J I I I I 
II  II  I I I I | | | II  III  I 


Hc\jro.3-ir>vi)fs'-cocr‘0r-ic\jf0.d-in^or-a) 

«H  H H H H rl  rl  rl  H 

\\\\\\\  \ \ 


CT'o»Hc\iro^-ir)vDC'® 

«H<\J(\jC\jC\JC\JC\JC\IC\JC\J 

\\\\\\\\\\ 


cr'Or-ic\jiocrinvors-cocr'  oh  oj 
CM  fO  H H H 

\\\\\\\\\\\\\\ 


ro^-invor^®a'Or-icMro3-iDvoh'®a'o*H 

i*H«HtHC\JC\JC\IOJ<\l(\JCVJC\JCVJC\JfOfO 

\\\\\\\\\\\\\\\\\\\ 

(\JC\JC\JC\JC\JC\JOJC\JC\JC\JC\JC\JCMC\JOJ(\JC\JC\JCM 


^ jm  non 


i\j  fO  j in  vQ  Is- 


OJ  ro  IT)  vO  Is- 


f\j  j-  in  v0  Is- 


(\ifOjiOvi)isH(MrojinvDrs 


<\j  ro  j in  © n 


OJ  ro  J in  vO  N 


cm  ro  in  ® r- 


124 


UST  PERIOD  818  1527  < 4405  2069  > 952  1410  < 4196  1928 


IT 

vO 


<T» 

\C 


O 

o 


co  in 
in  io 
ro 


Ll  Ll 
O O 


co  co 

U.  Lj 

cr  cr 


q o 
o o 


tr  cr 
u.  u. 
Q.  a 


co  co 

O JD 


•H  on 
O'  O' 


a:  cr 
< < 
LJ  Lj 


Ll  li. 
O O 


co  co 

Lj  L_ 

cr  cr 


O'  »H 
o j- 

cw 


o a 
o o 

cr  cr 

Lj  U/ 
0-  Q_ 


CO  CO 
O 3 


CO  CO 
rl  in 


II  II 

co  co 

h-  h- 

2 2 
UJ  Ll. 

a o 


O C H 
Ll  O 

O h 15  U 
0 2 2 
j y-i  z 
uJ  u.  o o 

> U.  4 u 

Lj  lj  co 
_j  i r 
uJH  o 
u I I u. 
X •-  O 
I—  H 

cr  _j  r~ 

2 cJ  >-  O' 

a < n 

2 CJ 

uj  x cr 
O Lj  O 
2 _J  U. 


2 2 
UJ  UJ 
O O 

LJ  O 
CJ  CJ 
<.  <a 

CD  CD 
2 2 


cr  lj 
o > 

X Uj 


x a 

CD  Q 


CD  ID  ID  ID 


<3  d 
X X 
CJ  CJ 


X X 
CJ  CJ 


125 


MORMNG  ACCIDENTS:  20.3* 


20.3*  CHANGE  IN  MORNING  ACCIDENTS 

-O.A*  CHANGE  IN  NIGHT+MIDDAY  ACCIDENTS  A5.7*  CHANGE  IN  MORNING  ACCIDENTS 
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-7.7*  CHANGE  IN  El/FNING  ACCIDENTS 
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